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spectroscopiè measurements. — 9. Other evaluations and reviews. — 
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Nomenclature 


Definitions of the principal symbols are given below. Other symbols are defined 
in the appropriate sections. 


a = 2re*/hc Fine structure constant. 

B Magnetic induction (gauss). 

e Velocity of light (em/s unless otherwise stated). 

e Electronic charge (esu). 

ele Electronic charge (emu). 

= Neje Faraday constant. 

Vp = 4rp,/h Gyromagnetic ratio of proton. 

vy Apparent gyromagnetic ratio of proton in oil. 

H Magnetic field strength (oersted). 

in Planck’s constant (ergs). 

h = h/2x 

k Boltzmann constant. 

A Wave-length (cm). 

Agl As Conversion factor from Siegbahn scale to absolute units, i.e., 
1 XU = (4;/4;) cm. 

x Wave-lengths in terms of Siegbahn scale, i.e., A’= Agm/(Ag/As)- 

M Atomic or molecular mass (amu).. 

Mj= Nm, Atomic mass of proton. 

m Mass of electron (g). 

My : Mass of proton (g). 

Ms Magnetic moment of electron (emu). 

Un = eh/4am,c Nuclear magneton (emu). 

Ho = eh/4ame Bohr magneton (emu). 

Ho Magnetic moment of proton (emu). 

My Apparent magnetic moment of proton in oil. 

N Avogadro’s number. 

v= w/20 Frequency (8-1). 

Ves Vos Va See wi, We, @, below. 

Vu Hyperfine splitting of hydrogen. 

R= Nk Gas constant (erg/mole °C). 

Rao = 2n*met/ch® Rydberg constant for infinite mass. 


Ry = Ralm,/(m,+m,] Rydberg constant for hydrogen. 

V Voltage (volt). 

o, = He/m,¢ Cyclotron resonance frequency of proton (57). 
wo, = He/me Cyclotron resonance frequency of electron (8-1). 
oO, = jy; Nuclear resonance frequency of proton (8). 
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A) Object. 
1. — Reasons for present study. 


From time to time evaluations of the atomic constants have been made 
by various reviewers. Two of the most recent examples are those of DUMoND 
and CoHEN [1] and BEARDEN, EARLE, MINKOWSKI and THOMSEN [2]. One 
area in which their results caused great concern was X-ray crystallography. 
In the past, X-ray wave-lengths were commonly expressed on the Siegbahn 
scale, i.e., in XU; one XU was approximately, but not exactly, one milli- 
angstrom (10-!! em). After long consultation with the principal workers in 
this field a committee headed by Sir LAWRENCE BRAGG[3] recommended 
1.00202 as the best value for the conversion factor from XU to absolute units, 
based on many direct measurements with ruled gratings by various inves- 
tigators. This figure was thought to be correct to within .003%. On the other 
hand, both the DuMonDp and CoHEN and the BEMT [2] evaluations gave ap- 
proximately 1.002 06. 

This discrepancy has led to requests of one of the authors (J. A. B.) for a 
new direct determination of /,/A,. A preliminary study led to the conclusion 
that there were no new techniques available which would assure any major 
improvement in such a measurement and that it would be extremely difficult 
to obtain a significant increase in accuracy. At the same time a study was 
made of previous direct experimental determinations; no indications were 
found of any important systematic errors or of any uncertainties in the theory 
of ruled grating measurements as applied to X-rays. Thus the one remaining 
hope of resolving the dilemma seemed to lie in a new survey of the best values 
of the constants, which was therefore begun. 

The A,/A, question is actually only one aspect of a problem raised by the 
two above-mentioned evaluations. In both cases it was found that the errors 
computed by external consistency were about 2.5 times those obtained by 
internal consistency. This can only mean that the accuracy of some or all 
experiments has been seriously overestimated. If this occurred more or less 
uniformly, no harm has been done, since the stated errors were based on ex- 
ternal consistency. However, such a viewpoint is certainly an optimistic one. 
If the discrepancy is due primarily to a particular experiment or group of 
experiments, the results of the evaluations may be significantly in error. 
(Another aspect of the discrepancy concerns the best value of Avogadro’s 
number N; this has been discussed at length in an exchange of letters between 
STRAUMANIS [4] and DuMonp and CoHEN [5]). 
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A further reason for a new study of the constants is the experimental data 
reported since the DuMoNnpD and CoHEN evaluation was completed in 1952. 
For example, there have been several determinations of the velocity of light, 
an X-ray experiment on calcite, and a measurement of the anomalous mag- 
netic moment of the electron. While none of these results by itself implied 
any major change, together they seemed to warrant a new survey of the 
situation. 

A final motive of the present report has been to attempt to summarize 
in one place the principal experiments which have been employed in calculating 
atomic constants. This should be of particular interest to textbook writers 
and physics teachers in showing just how the traditionally quoted measure- 
ments fit into the present picture. 


2. — Previous evaluations. 


In the early work in atomic physics there was a tendency to select one best expe- 
riment as fixing a particular constant or ratio. For example, Miilikan’s oil drop expe- 
riment was long regarded as giving the most reliable value of e. This approach gra- 
dually changed to a more critical review of all pertinent data. The International Cri- 
tical Tables [6] contained a list of the principal physical constants compiled in 1923; 
a few years later an article by HENNING and JAEGER in Handbuch der Physik [7] gave 
a more extended discussion of the subject. There was no adequate treatment of the 
errors in either case. 

The most important single advance in the field came in 1929 when BirGE [8] published 
a comprehensive evaluation of the physical constants with detailed discussion of the 
errors involved. In 1932 he made two additional contributions [9], an excellent summary 
of the theory of errors and a resurvey of the value of e, h, e/m, and x. The former 
provided the foundation for much of the subsequent work; it included discussions of 
the desirability of using probable error as a universal criterion, the distinction between 
internal and external consistency, and the use of the method of least squares. The 
latter had been previously employed by Bonp [10] for simultaneous determination of 
e and h in a study of atomic constants. Birge’s later work included a second com- 
prehensive survey [11] published in 1942 and a new evaluation of e and N in 1945 [12], 
as well as numerous other contributions on more specialized topics. 

Birge’s method was usually one of successive approximation. He began with the 
most accurately known constants, which were independent or only slightly dependent 
on the others. After fixing these, he used them consistently in subsequent calculations 
and worked his way through the less precise data step by step to determine all the 
required constants. The least squares method was sometimes used to simultaneously 
evaluate two constants, but not as a generalized procedure. While the procedure of 
successive approximations is open to some theoretical objections, it worked well in 
practice and could be understood relatively easily. The same approach was also used 
in studies of atomic constants by von FRImsEN [13] and STILLE [14]. 

A new step was taken by DUNNINGTON [15] in 1939 by making several simultaneous 
least squares calculations of e, m, and h. His primary aim was to find the cause of 
the inconsistency of the data then available by omitting one questionable item in 
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each calculation. The generalized least square method was used in three studies by 
DuMonp and Conen [1, 16]; their chief purpose was to determine the best values of the 
constants, rather than search for discrepancies. The values given in the last report 
are now in general use. One of the major contributions of DuMonp and CoHEN was 
to emphasize the importance of correlations between the results of various experiments 
when the same factor contributes a significant part to the error of each one. If such 
correlations exist, it is necessary to eliminate them by reformulating the problem or 
to treat them by a special method developed by COHEN. 

In 1951 BEARDEN and Warts [17] published an evaluation of the atomic constants, 
using the generalized least squares method. This work was continued by BEARDEN, 
EARLE, MINKOWSKI and THOMSEN [2]. The results were presented only in condensed 
form, since they appeared almost simultaneously with the last DUMoNp and CoHEN 
report and involved no significant differences in the values of the constants. 


8. — General considerations. 


There are at least four primary benefits of these evaluations: (1) Perhaps 
the most obvious one is to introduce a single set of values which can be used 
by all workers in the field until it becomes obsolete and is superseded by a 
new study. (2) A reviewer presumably can take sufficient time to study care- 
fully all items in the literature which are likely to influence the best values 
of the constants. Even so, many reviews have omitted or misinterpreted some 
significant data; there is certainly no guarantee that the present study will 
be an exception. (3) A more subtle reason lies in the need for obtaining a self- 
consistent set of results. For example, a certain X-ray experiment performed 
several years ago gives h/e as the final result; but the value of c? is needed 
in the calculations leading to this. If one blindly uses such a figure for h/e 
together with a recent determination of ¢ to obtain e/he, it is equivalent to 
assuming two contradictory values of ¢ in different parts of the calculation. 
On the other hand, a careful evaluation of a set of physical constants is un- 
likely to contain such obvious inconsistencies. (4) A study of the constants 
may point out certain discrepancies, which can lead to important new exper- 
imental and theoretical investigations. 

Out of the various surveys the group referred to as the atomic constants 
has crystallized in relatively definite form. It is usually understood to include 
e, h, m, c, and N. Of these, ¢ has sometimes been treated as a variable to be 
adjusted and sometimes as an auxiliary constant, whose value is known with 
much greater accuracy than the others. The same has held for the Faraday 
constant / = Ne. In recent years, as a result of the possible correlation of 
errors pointed out by DUMoND and CoHEN [1], it has also become necessary 
to include among the atomic constants /,/A,, the conversion factor between 
absolute wavelengths and X-ray wave-lengths measured on the Siegbahn scale. 
This factor hardly belongs to the group in its own right, but must be treated 
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with them as long as X-ray experiments play an important part in the results. 
Hence the group may now be said to include six constants: e, h, m, c, N, 
and /,/A,. 

At this point a few words should perhaps be included regarding the basis 
on which errors are stated. In the present report, all uncertainties will be 
given as probable errors (P.E.), which means that there is a 50-50 chance that 
the true value lies within the stated limits. This was the criterion always 
used by BIRGE and employed by DuMoND and CoHEN in their first two reviews. 
In their most recent work they use standard deviation (1.48 P.E.); this dif- 
ference is primarily one of personal preference. 

Probable error is rigorously defined only for statistical errors, e.g., one 
resulting from the use of a Geiger counter for a finite counting time. When 
systematic errors are important, as is usually the case, many workers prefer 
to give a «limit of error». There are a few instances in which this term may 
be strictly applicable. For example, the width of a gap which is measured 
by means of very accurate feeler gauges is known to lie somewhere between 
two values; no further information is known on error distribution. However, 
in most cases of systematic error, there is a large subjective factor in the esti- 
mate of either P.E. or limit of error. There is ample evidence that the latter 
term as used in practice gives no ironclad guarantee that the true value must 
lie within the stated limits. Hence it seems quite as justifiable to use P.E. 
even in the case of systematic error. It may be defined as an estimate such 
that the experimenter feels it is equally probable that the actual systematic 
error is greater or is less than the P.E. (Theoretically a standard experimenter 
might be found by choosing one whose estimates on many experiments have 
proven high and low equally often.) The P.E. is frequently obtained from a 
stated limit of error by arbitrarily dividing by three. While this rule has not 
been followed blindly in the present report, it is approximated in many cases. 


B) Experiments. 
1. — Introduction. 


A brief description will now be given of some of the various types of exper- 
iments used to determine the atomic constants. Only a few of these are suf- 
ficiently accurate to play an important part in fixing the precise numerical 
values of the constants today. However, the others are still of value in con- 
firming the results of the more precise experiments or in pointing out possible 
sources of systematic error in them. The odds are about 20 to 1 that the 
result of an individual experiment should not differ from the true value by 
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more than three times its probable error. When such a disagreement exists 
and can not readily be explained, there is cause for concern, regardless of the 
actual magnitude of the probable error. Thus many older experiments are 
of more than historical interest. While the list given here is certainly not 
complete, it is believed to include the principal methods of determining the 
various constants over the past thirty years or longer. 

For convenience in discussion, the experiments may be divided into three 
groups. The first consists of a number of experiments, usually simple in prin- 
ciple, which can be explained essentially in terms of classical theory. The 
second involves experiments which require some of the elementary ideas of 
quantum theory. The third class involves more complicated problems in 
quantum mechanics as well as the use of quantum electrodynamics. 

As an aid to the reader who may wish to omit the earlier work, the exper- 
iments actually used in the present evaluation are starred (*). Most of these 
involve errors of the order of a few parts per million (ppm); the errors in the 
other experiments may be judged accordingly. 


2. — Classical experiments. 


2°1.* Velocity of light. — The velocity of light ¢ is a most important physical 
quantity and has a proved a difficult one to evaluate properly. Excellent dis- 
cussions of this topic include BirGE’s 1941 report [11] on the best values of 
the constants, a careful review by DORSEY [18] in 1944, and more recent survey 
articles by EssEN [19] and MULLIGAN [20]. Various reviewers have determined 
what they considered the most probable values of e; the results are shown 
in Table I. It is clear that these values have fluctuated rather drastically, 
often by several probable errors. 


TABLE I. 
3 Value Error 
Reviewer Date (km/s) (km/s) 
En heels i 1929 299 796 4 
TEA RI es ESA RO a 1941 776 4 
MORSE LS ee e ot oa a, ys 1944 7713: 10 
VONSERIESEN [PS] see n 1937 780 20 
SOCIA OO ta 1947 778 20 
DuMonp and CoHen [21]. . . . . . 1950 789.3 0.9 
BEARDEN and WarTs [17]... . ... 1950 790.0 0.7 
HOSSEIN O e e ints che iret 1951 790 — 
DuMonp and CoHEN [21] . . . . . 1952 79219 0.6 


The principal experimental values of the velocity of light are given in 
Table II. 
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TABLE II. 
Experimenter Date () Yate @ BARGES) Method 
(km/s) (km/s) 
Rosa and Dorsey (1) 1907 299 784 10 Ratio of units 
MERCIER.(C) ane 1923 782 30 Standing waves 
MICHELSON (3). . . . 1927 798 15 Rotating mirror 
MICHELSON, PEASE 

and PEARSON (4) . 1935 774 4 Rotating mirror 
MITTELSTAEDT (5) . . 1929 786 10 Kerr cell 
PERITI) hee Bae 1940 apni 10 Kerr cell 
ANDERSON (7) SHE 1937 "LITE 10 Kerr cell 
ANDERSON (8) . . . . 1941 776 6 Kerr cell 
HousTOUN (°) Say 1950 775 9 Crystal grating 
BERGSTRAND (19)... 1951 (9) 793.1 2 Kerr cell | 
RANK, RutH and VAN- | 

DER SSLUIS.(41) 1952 776 15 (©) Infra-red bands | 
RANK, SHEARER and 

WIGGINS (12) G2 Sy 1954 789.8 3 Infra-red bands 
PLYLER, BLAINE and | 

CONNOR (13)2 5.9 san. 1955 792 | 6 Infra-red bands 
Essen and Gorpon- 

RIT Eas) Rea 1948 792 3 MW cavity 
EISSENE(2°) TR, re 1950 792.5 1 MW variable cavity 
HANSEN and Bot (16) 1952 789.5 23 MW cavity 
EF ROOME (22) 25S te acc 1952 792.6 sl MW interferometer 
BROOME: (18) tl 1954 793.0 i MW interferometer | 
ELORMANK( o) ote at ead 1955 795.1 Be RF interferometer 
PASTAS ONM (2) osama ie 1949 792.4 2.4 Shoran (radar) . 
ASLAKSON (21)... . 1951 794.2 1.4 Shoran (radar) | 


(a 

(?) Values in first group (8 experiments) are obtained from BIRGE (*); in some cases cor- 

| rections to original values have been made. 

| (°) Errors in first group are probable errors as given by BIRGE; those in second are those 
given by experimenters. (Probable errors unless otherwise stated in text). 

) Mean value, reported by BERGSTRAND in 1951, of all his experiments to date. 

©) Original error revised by RANK in later paper (2). 


| 
) Date of publication. . 
| 


B. Rosa and N. E. DoRsEY: Bull. U.S. Bur. Stand., 3, 433 (1907). 

MERCIER: Ann. Phys. (Paris), 19, 248 (1923); 20, 5 (1923). 

A. MICHELSON: Astrophys. Journ., 54, 1 (1927). 

A. MICHELSON, F. G. PrasE and F. PEARSON: Astrophys. Journ., 82, 26 (1935). 

MITTELSTAEDT: Ann. Phys., 2, 285 (1929). 

HUTTEL: Ann. Phys., 37, 365 (1940). | 

.C. ANDERSON: Rev. Sci. Instr., 8, 239 (1937). 

.C. ANDERSON: Journ. Opt. Soc. Amer., 31, 187 (1941). 

R. A. Houstroun: Nature, 164, 1004 (1949); Proc. Roy. Soc. Edin., A 63 (Pt. I), 95 (1950). 

L. E. BERGSTRAND: Nature, 165, 405 (1950); Ark. f. Fys., 2, 119 (1950); 3, 479 (1951). 

D. H. RANK, R. P. RurTH and K. L. VAN DER SLUIS: Phys. Rev., 86, 799 (1952); 

Journ. Opt. Soc. Amer., 42, 693 (1952). 

(1?) D. H. RANK, J. N. SHEARER and T. A. WIGGINS: Phys. Rev., 94, 575 (1954). 

(3) E. K. PLYLER, L. R. BLAINE and W. S. Connor: Journ. Opt. Soc. Am., 45, 102 (1955). 

(1') L. Essen and A. C. GorDoN-SMITH: Proc. Roy. Soc., A 194, 348 (1948). 

(15) L. ESSEN: Proc. Roy. Bore A 204, 260 (1950). 

(@°) K. Bou: Phys. Rev., 80, 298 (1950); final corrected value of c was given in paper at 
on Frequency Measurements Conference, Washington, D.C., January 14, 1953. 

. FROOME: Proc. Roy. Soc., A 218, 123 (1952). 

. FROOME: Proc. Roy. Soc., A 223, 195 (1954). 

FLORMAN: Journ. Research NBS, 54, 335 (1955). 

ASLAKSON: Trasn. Amer. Geophys. Un., 30, 475 (1949); Nature, 164, 711 (1949). 
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The experiments at optical frequencies are all basically similar to Fizeau’s. 
toothed wheel experiment. (Brief descriptions are given by JENKINS and 
WHITE [22]). Basically this method employed a rapidly rotating toothed wheel 
which sends out pulses of light 7 seconds apart, each one containing many 
waves. The light covered a round trip over a known distance and was de- 
tected by the observer if it covered this path in some multiple of 7 and so. 
passed through an opening in the wheel on its return. In later experiments. 
the wheel was replaced by a many-sided rotating mirror or by a Kerr cell, 
which can be used as an electro-optical shutter at frequencies of the order: 
of 10 me. All of these experiments measure group velocity by determining 
the time to cover a given distance. Rosa and DORSEY determined c by measur- 
ing capacitance in electromagnetic units and calculating it in electrostatic 
units from geometrical dimensions, the ratio being c?. MERCIER used a 
rf. oscillator with a measured frequency of the order of 100 me and determined 
the wave length of these oscillations by measuring standing waves on wires. 
When this value was corrected to free space wave-length, it could be mul- 
tiplied by the frequency to obtain e; this is thus a direct measurement of 
phase velocity. 

The second group of experiments will now be considered in somewhat 
more detail: 


a) Houstoun performed an optical experiment using a piezo-electric 
crystal as a shutter; an rf field was applied which made it act as a diffraction 
grating during part of each cycle. The stated error of 9 km/s is purely statis- 
tical; no estimate of systematic errors is given. The method seems capable 
of considerable refinement. 


b) BERGSTRAND’s experiment employed the geodimeter, an instrument 
which he developed for measurement of length in surveying applications. It 
is based on a modulated Kerr cell; the photomultiplier tube used as de- 
tector was modulated by the same source (8.33 MHz). Elaborate precau- 
tions were used to eliminate or average out small errors, such as slight dif- 
ference in transmission of various parts of the Kerr cell at a given instant. 
A major improvement lay in using a much longer base line (about 7 km) than 
in previous Kerr cell work; this greatly reduced the effect of any small un- 
corrected time lags in the cell and detector. BERGSTRAND gave .2 km/s as 
the standard deviation of the result. This appears to be stated less conser- 
vatively than some other (e.g., ESSEN’s); but his result certainly seems to be 
by far the most accurate optical one. 


c) The next series of experiments to be considered involves a combi- 
nation of infra-red and microwave spectroscopy. For each vibrational state 
of a molecule, the rotational levels can be represented to high accuracy by 
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means of the rotational quantum number and two constants. If a large number 
of lines are measured for an infra-red vibration-rotation band (involving a 
transition between two vibrational level), these constants may be determined 
by a least squares adjustment; they are thus given in terms of wave numbers. 
These constants may also be found in frequency units by a microwave measu- 
rement of a transition between two rotational levels in the same vibrational 
state. Thus the same quantity is known in both frequency units and wave 
numbers; the ratio of these two values must be c. 


RANK, RUTH and VANDER SLUIS carried out absorption measurement on 
an HCN band in the near infra-red (about .8 um). Their error was originally 
given as 6 km/s; but in a subsequent paper [23] it is stated that the grating 
used may have introduced a systematic error of 14 km/s. In the RANK, SHEARER 
and WIGGINS experiment, involving an HCN band at about 1.6 um, wave- 
lengths were measured with a Fabry-Perot interferometer; by small adjust- 
ments of air pressure in the interferometer it was possible to obtain an exact 
order of interference for the line to be measured. The error of 3 km/s appear 
to be the sum of four errors; if the probable error is calculated by the usual 
method, it will be less than 2 km/s. PLYLER, BLAINE and CONNOR measured 
a large number of lines in the CO band at 4.67 um; they used both absorption 
and emission lines. The error given is standard deviation. The last two 
measurements both achieved higher accuracy than Rank’s original experiment 
and both agree satisfactorily with other recent values of c. In all cases, the 
errors in microwave determinations of the band constants were much smaller 
than the spectroscopic errors. 


d) Some of the most precise recent measurements have been carried 
out using evacuated resonant cavities with microwaves. In some respects, 
the resonant cavity may be thought of as the microwave analogue of the Fabry- 
Perot interferometer; a good general discussion of this type of experiment is 
given by MULLIGAN [20]. Briefly, a cylindrical metal cavity is employed; in 
order to obtain resonance the length must be some integral number of half 
wave-lengths. However, this wave-length is not the free-space wave-length 
although it has the same order of magnitude. It may be calculated in terms 
of free-space wave-length by solving Maxwell’s equations; the relationship 
depends on which of many possible modes of oscillation is taking place. A 
small correction is applied for the finite Q of the cavity. Thus the free- 
space wave-length is found in terms of geometrical measurements, length and 
diameter; when the resonant frequency is determined, ¢ may be computed. 
(If two resonant frequencies are measured for different modes, it is possible 
to eliminate one geometrical measurement, e.g., diameter.) This method es- 
sentially measures phase velocity; it eliminates the correction for index of 
refraction of air since the experiment is performed in vacuum. 
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Essen and GORDON-SMITH used a fixed cavity with the TM,, and TM,,, 
modes of oscillation. A further refinement was introduced in the subsequent 
EssEN experiment by having an adjustable piston at one end of the cylinder. 
The modes used were TE,,, and TE,,,, where p took various small integral 
values. In order to determine wave-length it was primarily important to know 
the difference in length between two resonant positions of the piston; the 
absolute length was not critical. This was an important feature since the 
results indicated that the electrical diameter might be about 2 um greater 
than the geometrical value. ESSEN used two methods to correct for surface 
films and took the mean; there is a difference of about 1 km/s between these 
values. Values of c determined at three different frequencies checked to within 
.4km/s. The probable error of 1km/s appears to be conservative. 

HANSEN and BoL performed an experiment similar to that of Essen and 
GORDON-SMITH. Despite the great care and ingenuity they used, there seem 
to be at least two possible sources of serious error. First the upper plate of 
the cavity had a mass of about 28 kg and produced considerable compression 
in the spacer rods as well as dents in the lower plate, which required important 
corrections. Secondly there may have been errors due to surface effects. To 
check this point, BOL repeated the experiment, taking special precautions to 
remove impurities, and found no appreciable change; this does not preclude 
a surface layer of silver of high resistivity, as found by Essen. Essen mi- 
nimized this factor by using a variable piston and used a counterweight with 
the piston to avoid appreciable strain in his measuring gauge. Consequently 
there seems good reason to prefer ESSEN’s result and to conclude that Bow’s 
probable error should be raised considerably. 


e) Several recent experiments are based on microwave interferometers. 
FROOME’s first experiment used essentially a Michelson interferometer. The 
method is simple in principle and capable of high accuracy. However, it is 
complicated by diffraction effects, since the dimensions of the components 
are only a few wave-lengths. It is designed to employ spherical waves, but 
some corrections are necessary and the calculations are rather involved. 


FROOME’s second experiment, which employed a more symmetric arran- 
gement of components, is shown in highly simplified form in Fig. 1. The car- 
riage was set so as to obtain an interference minimum and then moved through 
about 81 wave lengths (about 1m) and carefully set for a minimum. The 
distance was measured in terms of very accurately known length gauges, a 
micrometer being used for a small final adjustment. This was done for a series 
of ten different transmitter-receiver distances in order to obtain the diffraction 
correction. To further investigate diffraction effects, three series of measure- 
ments were made, one with convex plastic lenses in the horns, one with half 
lenses, and one without lens. Slightly different methods of making the dif- 
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fraction corrections were tried. Although these corrections amounted to as. 

much as 100 km/s, the mean of different series taken with various lenses and 

various methods of diffraction correction 

peo checked within about 1km/s. The final 

spari value given is based solely on the mea- 

surements using convex lenses. The stated 

error of .3 km/s is apparently based on 

standard deviations; this would imply a 
probable error of only .2 km/s. 

The FLORMAN experiment used the 

interferometer principle at a much lower 

Fig. 1. frequency, 173 MHz. Since the wave- 

length is thus about 1.73 cm, it was de- 

cided to perform the experiment on a flat dry lake bed. Two receivers were 


ATTENUATOR 


TRANSMITTING HORNS 
RECEIVING HORN. 
DETECTOR 


1.5 mT07.5m 
MOVABLE CARRIAG 


placed 1.5km apart, the distance being measured to about 1ppm. A 


portable transmitter was then moved in a path encircling both and the max- 
imum phase difference was noted. The 1.5 km distance was then divided by 
half of this difference (in seconds) to obtain the velocity of the waves. Cor- 
rections were made for ground effects and for index of refraction of air. The 
error was given as 3.1 km/s, including a systematic error of .7 km/s and a 95%. 
confidence interval for the statistical measurements. This implies a probable 
error of about 1 km/s. 


f) ASLAKSON has made several determinations of c by means of « Shoran » 
(radar) measurements of a number of base lines ranging up to 235 miles. This 
is apparently the only direct radio measurement of group velocity. The errors 
reported are statistical mean deviations. The principal error in the 1949 measure- 
ment was a «signal intensity correction », which was substantially eliminated 
in the 1951 work. The probable error of the second experiment appears to 
be about 2 km/s. 


It is evident that all recent values lie in the vicinity of 299793 km/s. The 


only glaring exception is the HANSEN and BoL result, which differs from this 
figure by about ten times its probable error. Reasons for this discrepancy 
have already been discussed and appear to provide a plausible explanation. 


On the other hand, the older experiments, which are mostly optical ones, yielded 
Birge’s average of 299776 km/s. Among the optical experiments, only the Mittel- 
staedt and 1926 Michelson values may be considered in reasonable agreement with 
recent results. 

Of course, the older experiments were much less accurate; for example, in the 
Michelson, Pease, and. Pearson experiment, individual determinations were scattered 
over a range of 100 km/s. Systematic errors may have been considerably larger than 
estimated. BERGSTRAND [24] has recently studied the older optical work and reevaluated 
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the errors. For example, in the Michelson experiments he considers the variation in 
local atmosphere due to the expanding air blast from the turbine which drove the 
rotating mirror; an error of 8 km/s is ascribed to this in the evacuated pipe deter- 
mination. However, it should be noted that Birge obtained an error of only 1.4 km/s 
by external consistency; this error is primarily determined by the agreement between 
different experiments. Even if all six older optical experiments are arbitrarily assigned 
equal weights, the resulting value of ¢ is raised to only (299779 + 3) km/s, the error 
being determined by external consistency. Thus, while it is relatively easy to find 
additional possible sources of error in certain experiments, it is extremely difficult to 
explain why the statistical average was so far in error. 

So far as is known to the writers, the only partial explanation put forward as to 
the direction of the errors in the older experiments was given by Daynorr [25]. He 
has pointed out that Anderson used the group index of refraction for 5500 A to correct 
results to vacuum although the photocell and light source employed probably had 
their maxima at a considerably lower wave-length. Correcting for this might increase 
the Anderson result in the order of 1-5 km/s. This point is well taken, but can at 
most account for a small part of the discrepancy. 

Attempts to reconcile optical and microwave values of c thus depend chiefly on 
accepting Bergstrand’s result and rejecting most other optical experiments. (There 
is some question as to whether the values obtained from infra-red bands should be 
classed as optical or microwave ones; they would appear to be essentially optical ones 
if the relation H = hy can be regarded as certain). Hence there would seem some small 
room for doubt as to possible dispersion of light in vacuum. 

Under the circumstances, it seems worthwhile to briefly mention some indirect 
evidence against such a variation. a) Measurements of velocity of y-rays by LUCKEY 
and WerIL [26] and CLELAND and JASTRAWM [27], with errors of the order of 1%, show 
no evidence of a difference in c at these extremely short wavelengths. d) Observations 
‘on eclipsing binary stars [28] fail to show any conclusive evidence of dispersion in the 
visible region. e) The Rydberg constant, which varies inversely as c, is used in both 
the optical and X-ray regions; this provides a rather sensitive test. 

Thus, while the direct evidence for any variation in ¢ seems inconclusive, the 
indirect evidence against it appears rather strong. The question then arises as to whether 
the older optical values should be considered at all. The only ones in this group which 
are comparable in accuracy with any recent results are the Michelson, Pease and Pearson 
and the Anderson 1941 values. In both cases the errors were revised by Birge; in fact, 
the MPP paper [see (4) in Table II] did not contain even an estimate of systematic 
error. Dorsey took a more pessimistic view than Birge in both cases; he estimated that 
the MPP systematic error might be as high as 10 km/s. In light of our present know- 
ledge, it would seem that Dorsey’s view is fully justified. Thus, in the present eva- 
luation of €, it seems reasonable to include the infrared experiments (with errors of 3 
and 4 km/s), but to reject all optical values except Bergstrand’s. 


If the experiments were weighted with the probable errors assigned by the 
experimenters, ¢ would primarily be determined by the BERGSTRAND, HANSEN 
and Bor, and FROOME results. However, it seems preferable to adopt a system 
of somewhat arbitrary weights based on two factors: a) The stated probable 
errors with some modifications due to the reasons discussed above, particularly 
in the case of HANSEN and Bot, and 6) the difference in methods used. It 
seems desirable to give comparable weight to several different methods (Kerr 
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cell, microwave cavity, microwave interferometer, radar, infra-red bands) even 
though not strictly justified on the basis of probable errors; this procedure 
gives some insurance against completely unsuspected systematic errors. The 
final evaluation is shown in Table III. In the case of RANK et al., ESSEN, 
FROOME, and ASLAKSON, all but their last results are omitted and the entire 
weight given to the most recent values. Probable errors, computed from 
statements of the experimenters, are included only for comparative purposes. 


TABLE III. 
È t Value c; ht pine 
Experimenter (km/s) P.E. Weight w,; | Deviation d, 
BERGSTLERAND hea I e 299793.1 .14 5 + .3 
RANK, SHEARER and WIGGINS 789.8 3 1 =:3.0 
PLYLER, BLAINE and CONNOR . 792 4 1 eee 
ESSN oth A aaah EA ee hone 792.5 1 3 — 3 
HANSEN, and. (BOL! Giese. te 789.5 ne 1 — 3.3 
PROOME pets i IE 793.0 Ye 5 + .2 
WEORNLANT PI ovens teed ee tear bie 795.1 ] 2 + 2.3 
A SUAISON Gs oc! coat sth IE 794.2 1 +1.4 
Mean = Ywy,c;/X wi = 299 792.82 km/s 
(P.E. ext = -6745V D wyd?/(8 —1) Sw; = .34 


If the HANSEN and BoL value is entirely omitted the result will be raised 
to 299793.0 km/s. However, it will be retained in the present report and the 
result given as 


e = (299792.8 + .4)km/s (1 ppm). 


If the above analysis is correct, c is probably the most accurately known 
of the principal atomic constants. It is therefore convenient to treat it as an 
auxiliary constant (i.e. consider it as exactly known); it is then assumed that 
experimental errors are uncorrelated through the use of c. In discussing the 
various experiments in the following sections, it will always be made clear 
how ¢ is involved, if at all. However, the above numerical values of e will 
be used in stating the final results in cases where this will put the results more 
nearly in the form given by the experimenters. 

It should be emphasized once more that the direct evidence against a varia- 
tion of c between optical and microwave values rests primarily on BERGSTRAND’S 
experiment. Since the microwave values are in such good agreement between 
themselves, further microwave work will probably not throw any additional 
light on this question. On the other hand, one more direct optical measure- 
ment confirming BERGSTRAND would do much to settle any lingering doubts. 
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2:2. Electrical standards. — The older studies of atomic constants usually ineluded 
the conversion factors between international electrical units and the corresponding 
absolute quantities. On January 1, 1948 the old international system was abolished 
and all electrical measurements since then have been in terms of absolute units. At 
first glance it might seem that the experiments giving the old conversion ratios need 
no longer be considered in evaluating atomie constants. However, while the change 
has undoubtedly been a desirable step and reduced the multiplicity of electrical units, 
our knowledge of the constants still depends on the same experiments, regardless of 
this change in units. There are two reasons why they must be considered. 

In the first place these results still contribute to the error in other experiments 
involving electrical measurements. The absolute electrical units are defined in terms 
of the mks mechanical units with the added provision that the permeability of free 
space is defined as 47-10-’. Precise determination of absolute units in terms of these 
definitions is a very painstaking process and, in general, involves errors greater than 
those encountered in the intercomparison of electrical units. Thus the errors involved 
in precision electrical measurements hinge on the same experiments as before. 

The second factor, which does not appear to have been considered seriously hereto- 
fore, is the possible correlation of various experiments through the common use of 
electrical standards. This would not occur with the present system if each experimenter 
made his own determination of absolute electrical units as needed. However, in practice, 
all workers in the United States depend on the Bureau of Standards for calibration 
of electrical standards. An error in these standards may thus be common to all such 
experiments (DuMoNp [29] has pointed out the dangers of such correlations in 
the case of /,/A,). With the precision now attainable in some experiments, the pos- 
sible effect of this correlation should not be overlooked. 

The various absolute measurements are summarized in reviews by Curtis [30] and 
SILSBEE [31]. The determination of the ohm consists essentially of two steps. First 
a self-inductance or mutual inductance is very carefully constructed and its dimensions 
accurately measured. Since the permeability of free space is known by definition, the 
value of the inductance can be calculated from its geometry. Hence the reactance at 
any desired frequency may be computed. The second step consists in comparing this 
reactance with a standard resistance by means of a bridge or other suitable circuit; 
a number of different schemes have been used for this purpose. A recent measurement 
has been made by THomas, PETERSON, CooTER and Korrer [32] and was in good 
agreement with the previously accepted values. In his review Curtis quotes a weighted 
average error of 14 ppm for 8 experiments, which would imply an error of about 4 ppm 
in the mean value. This error appears to be a reasonable one. 

The absolute determination of the ampere depends essentially on measuring the 
force between two conductors carrying the same current. The usual procedure has 
been to use three coaxial coils in series with the middle one movable and attached to 
the arm of a balance; frequently two such sets were employed, one at each end of the 
arm. A recent experiment by DRIscoLL [33] uses a Pellat balance, which uses a large 
stationary horizontal coil and a small movable vertical coil near its center; the movable 
coil is subjected to a torque rather than a force. The measured current is passed through 
a standard resistor and used to calibrate the standard cell. However, the determination 
of other currents in terms of these standards does not depend on the knowledge of 
the ohm in absolute units, but only requires that the standard resistor and cell do not 
drift with time. The three most accurate measurements quoted by Curtis in his review 
yield an error of the mean of only 2 ppm, but some of the values he rejected were off 
in the order of 100 ppm. It should be noted that there was a discrepancy of 3 ppm 
between this average and the conversion factor finally adopted in converting standards: 
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to absolute units. The recent Driscoll determination differed from the accepted stan- 
dards by 7 ppm with the experimental error estimated at 8 ppm. (The discrepancy 
quoted in the Driscoll abstract was 19 ppm, but this was revised to include a correction 
for current distribution in the smaller coil due to a change of resistance with mecha- 
nical stress). The probable error in the electrical measurements, including small drifts 
in the average voltage of the standard cells, would appear to be about 5 ppm. To this 
there must be added the error in the acceleration of gravity, since it is necessary to 
know g to measure force in terms of known masses. The comparison of the various 
determinations of g given by DRYDEN [34] would seem to indicate a probable error of 
about 3 ppm. The total error then becomes 6 ppm. 

Thus the uncertainties in the various electrical units may be estimated as follows: 


ohm 4 ppm 
ampere 6 ppm 
volt 7 ppm 


‘The error in the volt is, of course, found by combining the errors of the first two mea- 
‘surements. These errors represent the highest accuracy obtainable in electrical mea- 
surements assuming perfect calibration of experimental measuring devices in terms 
of standards. In view of these errors, the correlation problem is not yet too serious 
a one, as will be evident from the following sections. However, it may soon be of 
greater importance and should be carefully kept in mind in future evaluations of the 
constants. 


2°3.* Electrochemical measurement of the faraday. — The faraday may be 
determined by measuring the amount of a known electrolyte which is depo- 
sited when a measured quantity of electric charge has flowed through a coulo- 
meter (electrolytic cell). Early experimental work on silver and iodine was 
carried out by Rosa, BATES, and VINAL [35, 36]. At the time, the ampere was 
legally defined in terms of the electrolysis of silver, and there was strong em- 
phasis on obtaining excellent reproducibility. This was indeed attained but 
there was a serious disagreement between the silver and iodine values of the 
faraday. There are at least four suspected sources of error: (1) the atomic 
weights used in the calculations may be incorrect, (2) the material deposited 
on the cathode may have inclusions of small pockets of liquid, (3) some of 
the material deposited on the cathode may redissolve in the solution, and 
(4) there may be partial isotopic separation. 

The value of the faraday in the old international units depended only 
on the atomic weight of silver and the definition of the ampere. Hence in 
principle the experimental part consisted of calibrating standard cells in terms 
of an ideal silver coulometer and holding these standards constant through 
the years. Taking the atomic weight of silver [37] as 107.89 and using the 
conversion factor of .999835 abs. ampere/NBS int. ampere [38] then gives 
F = 9651.3 emu/eq. (The error in atomic weight may be of the order of 
-01°% due to uncertainty as to the relative abundance of the two isotopes of 
silver.) Making similar corrections to the iodine values of the faraday [35] 
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using recent mass spectroscopic data for the atomic weight [39] yields 
F = 9651.8 emu/eq. Iodine has only one isotope and the error in its atomic 
weight appears negligible here. Also, in this experiment, a careful effort was 
made to eliminate errors of types (2) and (3) by checking the loss of iodine 
at the anode against the amount gained at the cathode. The limit of error 
for the chemical measurements appears to be about 50 ppm, which implies 
a probable error of roughly 20 ppm. When this is combined with uncertainties 
in electrical measurements, the resulting probable error should be in the order 
of 30 ppm. 

CRAIG and HOFFMAN [40] have recently made a new determination of the 
faraday, using the electrolysis of sodium oxalate (Na.0,0,). H, gas was 
given off at the cathode, there being no deposit in the usual sense. They weighed 
the total amount dissolved and determined by titration the amount remaining 
after a run. This procedure obviously eliminates errors of types (2) and (3); 
since titration depends on the number of moles but not the molecular weight, 
it also removes errors of type (4). It is, of course, necessary to know the exact 
purity of the compound dissolved. CRAIG and HOFFMAN used samples with 
purity certified as 99.96%. This involves an uncertainty of about 50 ppm. 
The impurities are analyzed and found to be inert in the reaction. At the 
time of the experiment there was a serious discrepancy between the chemical 
atomic weight of sodium and the mass spectroscopic value; this has now been 
removed by a revision of the chemical value [41]. Hence the Craig and Hoff- 
man value becomes 9651.1 emu/eq. The present situation is summarized by 
them [42] as follows: (1) The purity of the sodium oxalate is thought to be 
slightly lower than the certified value. (2) There is some variation in F with 
current in the oxalate coulometer (due perhaps to a competing reaction) which 
suggests the need for extrapolation to infinite current. (3) These corrections 
are in opposite directions and it is impossible at the present time to estimate 
the correction of the combined effect. 

From the above discussion, it seems evident that all three experiments are 
in reasonable agreement. The precise errors are difficult to estimate. If the three 
are averaged (giving half weight to the silver value because of the greater un- 
certainty in atomic weight), the resulting value is (9651.42 + .15) emu/eq, 
where the error is computed by external consistency. Since the weighting is 
somewhat arbitrary and since all the experiments have some errors in common, 
for example the accuracy of electrical standards in absolute units, the error 
will be doubled and the result given as 

Ne 


we (9651.4 + .3) emu/eq (30 ppm). 


24. Millikan oil drop. — This well-known experiment consisted in measuring the 
charge on a single oil drop by means of the force exerted on it by an electrostatic field. 
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The rate of fall of the drop through a viscous fluid (air) was then given by Stokes’ Law 
plus small experimentally determined corrections. The charge so determined was 
always approximately an integral multiple of an elementary charge, i.e., the charge 
of the electron. The experiment is described in many general physics texts, e.g., SEARS 
and ZEMANSKY [43]. For some time newer values derived from X-ray experiments of 
BEARDEN and others were in serious conflict with the oil drop results; the problem 
was apparently resolved by new determinations of the viscosity of air. The situation 
as of 1945 is well summarized by BIrGE [12]; he obtained e= (4.8130+.0075)-10—° esu 
as an average of oil drop determinations by MILLIKAN [44], HoPPER and Lapy [45], 
BAcCKLIN and FLEMBERG [46], and IsHipa, FuKUSHIMA and SUETSUGU [47]. 


Subsequently, according to a footnote in a later paper by Hopper [48], IsHIDA 
discovered an error of about 1% in voltage measurement. This would lower the Birge 
value to about (4.804 + .008)-10-1° esu, in excellent agreement with other values of e. 
Unfortunately in this same report Hopper discusses a correction to Stokes’ Law due 
to the proximity of the condenser plates and, using this correction in a new experiment, 
obtains e= 4.8259-10-! esu (with the Birge figure for viscosity of air). While he claims- 
greater internal consistency with this correction, it leaves the various experiments in 
serious disagreement. He believes the results may imply a lower value of the viscosity 
of air than now accepted. Another possible source of error pointed out by IsHIDA [49] 
lies in the facts that many drops are oblate spheroids rather than spheres. 


If Hopper’s corrections are accepted, the resulting mean appears to be close to 
4.82 esu. On the basis of internal consistency as well as comparison with other methods. 
of finding e, there seems to be reason to suspect some systematic error in his work and 
to reject his result. Since the oil drop method has comparatively low accuracy, it 
does not seem necessary in the present report to pursue the question further and at- 
tempt a more critical analysis of the conflicting data. 


2:5. e/me by deflection methods. — The so-called e/m experiments may be divided 
into three classes, as has been pointed out by DuMonp and Corner [16]. One class gives. 
a direct determination of e/me (e/m in emu); the other two yield e/me only with the 
use of the faraday. The measurement of the faraday is subject to numerous errors, 
as previously discussed; hence it can not be treated as an accurately known auxiliary 
constant. Thus it is important to separate these three categories. The deflection exper- 
iments now to be considered form part of the first class; the other data will be dis- 
cussed in the proper section of Part B, 3. The following data include corrections to 
absolute electrical units; most of the corrected figures are taken from DuMonp and. 
CoHEN [16]. 

Early measurements were carried out in the well-known work of Sir J. J. THom- 
son [50]. The principle of his method is described in many texts, e.g., SEARS and ZE- 
MANSKY [43]. It consists essentially in measuring the deflection of a beam of electrons: 
in a known magnetic field and then finding the electrostatic field which will just cancel 
this deflection. 


To secure greater accuracy the obvious modification in Thomson’s method is to 
observe the deflection of electrons through a circle, or through many revolutions, rather 
than a small are. This is facilitated by the well-known fact that an electron in a 
uniform magnetic field revolves in a circle with uniform angular velocity w = He/mce,. 
independent of its velocity. (Since this is the principle of cyclotron operation, this. 
expression is often referred to as the cyclotron frequency of the electron). DUNNING- 
TON [51] devised a very careful experiment in which an electron beam traveled through 
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a narrow circular ring of 340° in a uniform magnetic field. The electrons were acce- 
lerated by a high frequency voltage at the beginning of the trajectory and either acce- 
lerated or decelerated by the same voltage source at the end; this gave a resonance dip 
in the current when the time of flight was exactly one cycle. The field was adjusted 
for resonance and measured by means of the current in the solenoid used; w was deter- 
mined from the known frequency and angle; thus e/me was calculated. The principal 
difficulty lay in the formation of surface charges and corresponding spurious potentials 
(sometimes 10 to 20 V) on the walls of the evacuated chamber. DUNNINGTON attempted 
to allow for this by extrapolating to infinite energy; this correction was of the order 
of 5 parts in 104. Saw [52] has criticized this on the basis of his own work with 
surface charges and questioned the assumptions on which the extrapolation was based. 
Dunnington’s result was (1.7598 + .0004)-107 emu/g; if Shaw’s criticism is correct, 
the error should be increased to perhaps .0010. 

Suaw [52] measured e/m in a rather complicated method employing a cylindrical 
condenser with a radial electrostatic field and an axial magnetic field. Under these 
conditions he showed that a beam of electrons would be focused by the combined fields 
if it traveled through a circular are of about 127°. In general, electrons struck the 
collector for two different velocities or for none; however, with the proper combination 
of electric and magnetic fields, they struck the collector for a unique velocity. (The 
velocity of the electrons was varied by changing the accelerating voltage). Under these 
conditions it was possible to calculate e/me from the measurements of the fields and 
the radius of the orbit. Shaw’s value was (1.7582 + .0013)-107 emu/g. 

Perry and CHAFFEE [53] measured e/me by a direct velocity determination. An 
electron beam was accelerated by a high voltage £, from which the energy of each 
electron could be calculated in terms of e. It was then passed between two pairs of 
electrodes a known distance apart, to which a radio frequency voltage was applied. 
Electrons came through to a fluorescent screen without deflection only if they covered 
this distance in exactly one cycle of the r.f. Since the frequency could be accurately 
measured, the electron velocity was thus obtained and e/me computed from the relation 
He/e = mv?/2 (with proper relativistic correction), where £ is in statvolts. The result 
was (1.7611 + .0010)-107 emu/g, which seems rather high. There may have been 
some error due to an axial magnetic field needed to focus the beam and secure suf- 
ficient intensity ; they describe this as an undesirable but necessary feature. KIRCHNER [54] 
performed a similar experiment with a resulting value of (1.7590 + .0009)-107 emu/g. 

The Busch method makes use of the focusing property of a magnetic field parallel 
to the direction of a beam of electrons. The electrons are accelerated through a known 
voltage and then emerge through a small aperture in the anode as a slightly divergent 
beam. Assume that the velocity is principally directed along the z axis with small 
random components in the zy plane. If a magnetic field is applied along the 2 axis, 
the electrons will revolve in circular orbits in the zy plane with angular velocity 
w = He/me. Since they all have nearly the same velocity v along the z axis, they will 
return to a focus on the axis after traveling a distance | = 22v/w. v is found from 
the accelerating voltage (as in the previous experiment) and e/me calculated. The 
method was first used by BuscH [55]. More accurate determinations by Wour [56] and 
GEODICKE [57] gave (1.7590 + .0018)-107 emu/g and (1.7587 + .0008)-10? emu/g re- 
spectively. 

If Dunnington’s error is increased as previously explained, all results except Wolf's 
have roughly the same error, viz., .0010-107. When all are weighted equally the re- 
sulting mean is 1.7593-107. The error is about .0004-10’ by internal consistency 
and slightly less by external consistency. Inclusion or omission of the Wolf result 
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makes no significant change. The result of all the deflection measurements may there- 
fore be taken as 


e/me = (1.7593 + .0004)-10? emu/g (230 ppm). 


2°6.* Avogadro’s number (X-ray crystal density method). — Avogadro’s 
number N is defined as the number of molecules per mole (gram molecular 
weight). Thus, for a substance with molecular weight M, the mass of each 
molecule is M/N. If this substance has a perfect cubic crystal structure with 
unit cell of length d containing f atoms its density o is fM/Nd?. For a crystal 
such as calcite, where the unit cell is a rhombohedron, d* must be multiplied 
by g, the volume of a rhombohedron with sides of unit length. (This can be 
calculated in terms of the dihedral angle x between two faces of the unit cell.) 
Hence N is given by 


N={fM[0dp. 


Experimentally the problem of using’ this relation to obtain N reduces to a 
measurement of d and @ (i.e. x) by means of X-rays and an accurate deter- 
mination of density and molecular weight. 

This work has been reviewed at some length by BEARDEN[58] and BIRGE [12]. 
Several different crystals have been employed. However, since the calcite 
determinations appear to be the best, only calcite will be used here. Accurate 
determinations have been made by BEARDEN [59] and BROGREN [60]. BEARDEN 
determined the density of calcite by weighing crystals in kerosene and using 
Archimedes’ Principle; the same crystals were later checked by the Bureau 
of Standards. The mean of these measurements together with one by DEFOE 
and CoMPTON gives the value 0 = 2.71029 g/cm’, which should be accurate 
to within 10 ppm. « is determined by reflecting X-rays from two faces of the 
crystal and g is then computed. BEARDEN found g = 1.09594 + .00001 at 
20 °C. BROGREN obtained a value 20 ppm higher. (He also calculated a value 
40 ppm higher based on a mean of five different measurements, including his 
own; but this appears to have given equal weight to less accurate determi- 
nations.) In the present calculation y will be raised 20 ppm to 1.09595 (a 
mean of the BEARDEN [59], BROGREN [60] and Tu [61] values) and assigned an 
error of 15 ppm. The molecular weight of calcite, when computed directly 
on the physical scale from recent mass spectroscopic date (see BROGREN [60] 
for details), is 100.119; the error is about 20 ppm due chiefly to uncertainty 
in the abundance ratio of calcium isotopes. However, the effective molecular 
weight is slightly higher due to impurities. BIRGE calculated from BEARDEN’S 
data that this correction should be 31 ppm. BROGREN used an average based 
on both BEARDEN and STRAUMANIS to obtain 50 ppm. In the present work 
a correction of 40 ppm will be made and assigned an error of 20 ppm. This 
gives M,,, = 100.123 + .003 (30 ppm). 
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The distance d is measured by observing X-ray diffraction from calcite, 
using the Bragg Law. This involves a rather complicated three-dimensional 
phenomenon and is well described by CoMpTON and ALLISON [62]. Briefly the 
results show that whenever a diffraction maximum occurs, it may be con- 
sidered as arising by direct reflection from 
parallel planes passed through the atom in some 
simple manner. These are the so-called Bragg 
planes; their intercept on each crystallographic 
axis is the reciprocal of some integer (usually 
small in practice) times the length of unit cell 
along that axis. Waves scattered from all atoms 
in such a plane are in phase since they are pro- 
duced by direct reflection; for a maximum it 
is also necessary that waves from each of the parallel planes reinforce one 
another. Fig. 2 shows that the necessary condition for this is 


nA = 2d sin 6, 


where » is an integer, d is the separation between planes, 4 is the wave-length, 
and @ is the angle between the plane and the incident beam. 

In early X-ray work it was soon apparent that the diffraction angle 0 for 
X-ray wave-lengths could be measured much more precisely than the value 
of d was known. An arbitrary value for the grating space of calcite, d' (Sieg- 
bahn) = 3029.51 XU at 20 °C, has therefore been assumed for X-ray wave- 
length measurements and crystallographic work; the XU is approximately 
10-!!' em. (This Siegbahn scale was originally based on rock salt with 
d'= 2814.00 XU at 18 °C, but it was soon found that calcite was a more 
suitable as a standard.) Since the value of d’ is assumed, it has no probable 
error, except possibly 5 ppm due to temperature and index of refraction cor- 
rections and to variations between samples of calcite crystals. This grating 
constant d’ may be regarded as the length of a cell containing one half of a 
calcite molecule (i.e. f = 4). It must still be multiplied by a conversion factor 
Ag], to convert it to absolute units. Since /,/A, is also needed in other X-ray 
experiments, its determination should be treated as a separate experiment to 
avoid correlation between the various data; this point has been emphasized 
by DuMonp and CoHEN [29]. The conversion factor A,/A, will be treated in 
the next section. The equation for N may be rewritten as 


N'= N(A,/A;)* = {Med , 


where d’ is the grating space in terms of the Siegbahn scale; in other words, 
N' would be the value of N if 1 XU were exactly 10-!! em. When the various 
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values given above are inserted, the result is 
N'= (6.06145 + .00022)-10?3 mole! (38 ppm). 


A rhombohedral crystal may also be referred to a hexagonal set of axes. 
The volume d*y must then be replaced by 3a?e 3/2, where c is measured along 
the optic axis and a lies in a plane perpendicular to it. This makes it pos- 
sible to compute N (or N’) without measuring the angle x. BROGREN used 
a «unit cell» containing 18 molecules and found a’= 4997.97 XU and ¢’= 
= 17.02599 XU (corrected to 20 °C), both having errors of 6 ppm. These 
results yield N’ = (6.06140 +.00021)-1023. VAN BERGEN [63] also measured 
erating constants in the hexagonal system. Using a «unit cell» of 36 mole- 
cules, he then obtained, at 20 °C, a'= (9959.93+ .05) XU and e’ = (8513.06+ 
+ .06) XU. The corresponding value of N’ is thus found to be (6.06136 + 
+ .000 21) - 102°. 

The agreement between results in the rhombohedral and hexagonal systems 
gives a good check on the various X-ray measurements involved. However, 
the major error, i.e., that of the molecular weight of calcite, is not eliminated. 
The mean value of N’ will be taken as 


N'= (6.0614 + .0002)-1023 mole (33 ppm) , 


with the error substantially the same as for the individual values above. 

This value may be compared with a value of NV’= (6.06178 + .00023)-1023 
obtained by BIRGE [12], using data for five different crystal—calcite, diamond, 
NaCl, LiF, and KCl. Since the value given above for calcite differs from 
Birge’s figure by almost twice his probable error, it may appear to be rather 
questionable to discard data for the other four crystals. However, a consi- 
derable part of the discrepancy is due to the rock-salt (NaCl) data, based on 
the work of Tu [61]. The available data on NaCl is notin good agreement. 
An early measurement of Siegbahn gave a value of d 126 ppm higher than 
Tu while STRAUMANIS obtained one 20 ppm lower. A more recent determi- 
nation by VAN BERGEN [63] yielded a value 78 ppm higher with the error 
quoted as onlu 11 ppm. These wide variations may be due to a mosaic structure 
of NaCl as described by RENNINGER [64]. If VAN BERGEN’s result were used 
instead of Tu’s, the N’ obtained from NaCl would be the lowest of the five 
crystals instead of the highest and the Birge average would drop to 6.06150. 
On the other hand, if all NaCl data are omitted his average would be 6.06164. 
Thus the discrepancy is reduced to one probable error or less. When it is also 
recalled that the present analysis for calcite includes the BROGREN and VAN 
BERGEN work as well as recent mass spectroscopic data on atomic weights 
and abundances, none of which were available to BIRGE, it seems quite justifi- 
able to adopt the present calcite value. 
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2°7.* Ruled grating measurement of X-ray wave-lengths. — The conversion 
factor 4,/4s between XU and absolute units (milliangstroms) is found directly 
from measurement of X-ray wave-lengths by ruled gratings. At X-ray wave- 
lengths the index of refraction is less than unity and therefore total reflection 
occurs for X-rays incident on a grating from air or vacuum. (The simple 
theory must be slightly modified if absorption is appreciable.) It is only 
with total reflection that a grating is practical for use with X-rays; this 
means that the incident beam must strike the grating at a small glancing 
angle of the order of 1°, which corresponds to 89° angle of incidence in optical 
terminology. The wave-length is calculated from the>usual optical formula 
NA.= dsin 0. 

EcKART [65] has investigated the possibility of systematic errors (due, for 
example, to refraction, multiple scattering, and surface waves) when the above 
equation is used for an ideal grating at X-ray wave-lengths. He concluded 
that no appreciable error arises from these causes. 

Precision grating measurements suitable for calculating /,/A, have been 
‘carried out by BEARDEN (5 results). BACKLIN, SODERMAN, and TYREN (4 re- 
sults). Most of these experiments have been discussed by BIRGE [12] and in 
greater detail by BEARDEN [58]; however, BIRGE gives only one of TYREN’s 
values, while BEARDEN’s paper omits these entierly. 

The various X-ray lines used in the above experiments were Cu-K,, Cu-K,, 
Cr-K,, Cr-K,, Si-K,, Al-K,, Mg-K,, and Na-K,, the wave-length region covering 
roughly 1 to 12 A. BEARDEN and BACKLIN used plane gratings and determined 
absolute wave-length by directly measuring the grating space d. SODERMAN 
and TYREN, on the other hand, employed concave gratings and compared 
higher orders of X-ray spectra with first order spark spectra of various ele- 
ments (e.g. the Lyman series of hydrogen-like ions, which could be calculated 
theoretically and required only a small experimental correction). The various 
results cover a spread of 130 ppm and show good agreement in view of the 
errors given, usually about 50 ppm. 

Using his own data plus BACKLIN’s and SODERMAN’s, BEARDEN [58] obtained 
a mean of ,/Z,=1.00203 + .00001 (10 ppm), while BIRGE [12] later gave 
1.00203 +.00002 (20 ppm) with the inclusion of TyREN’s Al-K, measurement. 
It should be noted that the first four BEARDEN experiments have some syste- 
matie errors in common (estimated at 40 ppm) and are thus somewhat 
overweighted in the above averages since they are considered independent 
observations. In the conclusion of his paper, TYREN [66] recommended 
1.00201 as the best average of all experiments and SIEGBAHN [67] supported 
this view. After prolonged study the X-ray Analysis Group of the Institute 
of Physics adopted 1.00202. In reporting this decision, BRAGG [3] states 
that the value is believed to be correct to within .003%; this is probably 
intended as a limit of error. In the present report, this last average will 
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be used with the probable error given by BIRGE, i.e. 
AglAs = 1.002 02 + .000 02 (20 ppm) . 


As BIRGE pointed out, the external consistency of the data suggests that this 
error is quite conservative; this remains true when the TYREN experiments 
are used. 


2°8.* Atomic mass of the proton and other atomic masses. — The atomic 
mass of the proton enters into the calculation of atomic constants at several 
points. It has been determined by two principal methods. The first is by 
means of mass spectrometric measurements of « doublets », i.e. two ions which 
have almost the same e/M value. Thus the corresponding mass spectrographic 
lines will lie close together; in practice, this means that the small difference 
in mass can be measured to high accuracy. A commonly used set of three 
such doublets which yields the mass of the proton in terms of the standard 
*O is (*H),—*D, ?D)s—™C, #C(*H),—*0. 

The second method forms an interesting contrast to the first; it makes. 
use of energy differences determined from nuclear reactions to compute mass. 
differences. While such measurements are not highly precise, they are only 
used to determine a small correction to the mass, i.e., the amount by which 
16Q/M,, differs from the integral value 16. Thus, when a large number of such 
reactions furnishing various paths from 1H to 1°O are used, a least squares solution 
can yield an accuracy comparable to that of mass spectroscopic determinations. 

Recent data on the atomic mass of hydrogen have been summarized in 
a review by DUCKwoRTH, HoGG, and PENNINGTON [39] and are shown in 
Table IV together with a new value computed by DRUMMOND [68]. 


TABLE IV. . 
; ; | Mz: 108 

Experimenter Date | Method (physical scale) 

BAINBRIDGE Wh.) «lag. e e 1947 Prewar mass spec- |1008128.3 + 2.8 
troscopic data 

WALD E RE OO 1951 Mass spectrometer 41 +2 
CoLLINS, NIER and JOHNSON 1952 | Mass spectrometer 46 +3 
OGATA and MATSUDA . . . 1953 Mass spectrometer 45 +2 
MaTTAUCH and BIERI . . . 1954 Mass spectrometer 45.9 + .5 
AIA INTA ARR ee Ce RC 1951 | Nuclear reactions A ORE 
DRUMMOND E eee uae 1955 Weighted average 44.0 + .7 


A good discussion of the work in mass spectroscopy and errors in the older 
measurements is found in the proceedings of the Symposium on Mass Spectro- 
scopy [69] at the Bureau of Standards in 1951. 
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It is evident that the older data used by BAINBRIDGE probably contained 
serious systematic errors while the last six values concur reasonably well. 
The good agreement between the mass spectroscopic work and the nuclear 
data is very encouraging. The DRUMMOND average is obtained from a least 
Squares evaluation of atomic masses, based on both nuclear and mass spectro- 
scopic data, including all the above work except that of MATTAUCH and BIERI. 

The present report is intended primarily as a study of atomic constants, 
rather than nuclear ones. Consequently the following scheme will be adopted 
+ without a more detailed study. The DRUMMOND value will be given a weight 
of three and the MATTAUCH and BIERI value a weight of unity. (This assigns a 
considerably higher weight to the DRUMMOND result than indicated by the 
relative errors, yet it appears justified on the ground that his study includes. 
a considerable number of experiments employing two quite different princi- 
ples.) The error will be conservatively set as .0000010, slightly higher than 
either the DRUMMOND or MATTAUCH and BIERI error. This accuracy of 1 ppm 
is still ample for the purpose at hand. Thus the result becomes 


M,, = (1.0081445 + 0000010) amu (physical scale) (1 ppm). 


The atomic mass of the proton is now readily found by subtracting the atomic 
mass of the electron, Nm. Since the latter is only a small correction, it need 
not be known to high accuracy. Thus the source of this datum is not critical 
and no appreciable error is introduced by it. The result is: 


M, = (1.0075957 + .0000010) amu (physical scale) (1 ppm). 


Other atomic masses used (except silver) have been taken from or checked 
against the DUCKWORTH, HoGG, and PENNINGTON data. While several values 
are often given by the latter, the differences are generally not important in 
the present applications. One point should perhaps be mentioned in connection 
with the atomic mass of oxygen and the conversion factor r from chemical 
to physical state. The chemical scale is based on taking the atomic mass of 
oxygen as exactly 16, but does not specify its source. According to NIER [70] 
the chemical mass unit will have a mass of 1.0002783 amu. on the phy- 
sical scale if atmospheric oxygen or oxygen from limestone is used as standard,. 
but only about 1.000268 if the source is water or iron ore. Thus the con- 
version factor is no longer a matter of experiment, but rather of convention; 
according to WICHERS [71], it seems likely that r = 1.000275 will soon be 
adopted as a conventional value. This figure will be employed in the present 
report. The uncertainty in abundance of oxygen isotopes affects the CRAIG 
and HOFFMAN measurement of the faraday and crystal measurements of Avo- 
gadro’s number, but is not critical in either case. 


292 J. A. BEARDEN and J. S. THOMSEN 


2°9.* Ideal gas constant. — While not directly needed in the evaluation 
of most of the atomic constants considered here, the ideal gas constant is 
needed in obtaining one important derived quantity, i.e., the Boltzmann 
constant. The value of pV for an ideal gas at 0 °C and zero pressure is given 
as (2271.16 + .04) J/mole (chemical scale) by the Subcommittee on Funda- 
mental Constants of the Committee on Physical Chemistry of the National 
Research Council [72], the error being the standard deviation. This is deter- 
mined by precision measurements on almost ideal gases to determine the virial 
coefficients and permit extrapolating to zero pressure. Converting to the phy- 
sical scale and changing to probable error gives 


(pV)2=°,., = (2271.78 + .03)-107 erg/mole . 


The absolute temperature at 0 °C is now known by definition. The Tenth 
General Conference on Weights and Measures [73] adopted a proposal origi- 
nated by KELVIN in 1854 and renewed by GIAQUE in 1949 to define the absolute 
temperature scale by a single fixed point. This was taken as the triple point 
of water and defined as 273.16°K. The zero of the Centigrade scale had 
previously been defined as .01 °C below the triple point [74]; hence 0 °C is 
now exactly 273.15 °K. Thus the gas constant R = pV/T becomes 


R = (8.3170 + .0001)-107 erg/deg mole (12 ppm). 


It should be noted that neither RT nor kT is really known any more ac- 
curately at 0°C as a result of the new definition of the temperature scale. 
Under the old system the product RT at 0 °C was known far more accurately 
than either R or 7. Since in practice correlation coefficients are not usually 
employed in calculating error, the error in RT at or near 0 °C appeared to 
be far greater than it actually was. With the new definition a good estimate 
of the error can be obtained without resort to correlation coefficients. 


3. — Quantum experiments. 


3:1. Photoelectric effect. — The first precision experiment to give h or h/e made 
use of the photoelectric phenomenon. The theory of this effect was first given by 
EINSTEIN. (See, for example, SEARS and ZEMANSKY [43] or MILLIKAN [44]). It is summa- 
rized by the equation 


MU ax e 
hope 108V 7) 


in which » is the frequency of the incident light, g is the work function of the material, 
and mv 1/2 represents the maximum kinetic energy of an emitted electron. This energy 
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is in turn equal to 108Ve/c, where V is the minimum retarding potential necessary to 
stop the photoelectric current. Since v = ¢/A it is only necessary to plot V as a fune- 
tion of 1/A; the slope of the resulting straight line is he?/e. In practice, the situation 
is complicated by varying amounts of impurities on the photo-sensitive surface and 
by photoelectric emission from other surfaces due to scattered light. 

MILLIKAN [75] performed a very careful experiment designed to eliminate these 
difficulties. He constructed an ingenious experimental setup in which a new surface 
could be machined on the photo-sensitive material while it remained in a vacuum 
tube. The materials used were sodium and lithium, which are sensitive in the visible 
region, while the anode was copper, which has a photoelectric threshold well in the 
ultra-violet; thus over a large range of frequencies scattered right had little effect. 
MILLIKAN obtained 6.57-10-?” as the value of h, using e = 4.774-10-1° and 300 as 
the conversion factor between statvolt and volt. The value of ¢ used in calculating fre- 
quencies is not given; it may perhaps be assumed that it was taken as 3-101 cm/s 
as in the voltage conversion. If the results are then corrected to absolute electrical 
units and the present value of e inserted, Millikan’s work gives 


h/e = (1.379 + .007)-10-1" erg s/esu. 


More recently Lukirsky and PRILEZEAV [76] performed an experiment with photo- 
electrons from Zn, Al, and Ni; they found 6.543-10-?? erg s as the value of h with an 
error of .1% to .2%. The «reverse current » due to emission from the anode was not 
eliminated, as in Millikan’s work, but experimental corrections were carefully made. 
Another photoelectric measurement was made by OLPIN [77]; he used sodium surface 
sensitized by a mixture of sulfur vapor and air, which was sensitive throughout the 
‘entire visible spectrum. His published value [78] of h was 6.541-10-27 « correct to at 
least three significant figures ». However, a similar run with pure sodium gave 6.6: 10-27; 
while this is admittedly less accurate due to smaller frequency range, it is not clear 
why it was neglected entirely. 

Since neither OLPIN nor LUKIRSKY and PRILEZEAV give the values of e and € 
which they assumed, it is impossible to determine precisely their experimental values 
of hje. In any event their results are roughly .4% lower than Millikan’s. While both 
results are given with lower errors than Millikan’s, it is not evident that there is any 
essential improvement to justify the higher weight. In fact, with the corrections required 
in the LuxKtrsky and PrILEZEAV work, the contrary might be expected. It should also 
be noted that both MirLigAN and Orrin obtained higher values of h with sodium, 
which they rejected. DUBRIDGE [79] and RornrR [80] have discussed the effect of thermal 
energies on the results. They point out that Einstein’s equation, holds strictly only 
at absolute zero and can be in error by several hundredths of a volt at room tem- 
perature; however, the error in the slope dV/dv may be considerably less. In light 
of ail these considerations Millikan’s value will be lowered by .2% and his error retained. 
Hence the result of the photoelectric measurements is 


h/e = (1.376 + .007)-10-1” erg s/esu (510 ppm). 


3:2. X-Ray photoelectric effect. — Photoelectric experiments may also be carried 
out with X-rays. In this case the electrons are ejected from the lower energy levels 
of the atom rather than the conduction band; consequently the equivalent « work 
function » is much greater. The energies of the emitted electrons are high enough so 
that they may be conveniently measured by the radius of curvature r in a uniform 
magnetic field. The velocity is then given by Hre/me (neglecting a small relativistic 
correction which is required); the corresponding energy is then (Hre)?/2mc?, The in- 
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cident frequency v is given by (¢/A’)(A,/A,) and the «work function » is approximately 
(he/A,)( A,/A4g), where A, is the absorption edge for the ejection of electrons from the 
particular level oli and the primes denote measurement on the Siegbahn scale.. 


The resulting equation is 
268 (1 1 
amp ld) = (Hrye\a ai)? 


2, may be determined from measurements on the absorption edge as done by KRET- 
SCHMAR [81] or eliminated by taking measurements at two values of 4’ as done by 
ROBINSON [82]. These experiments are discussed at greater length by DUNNINGTON [15],. 
who concludes that the latter method is more accurate. Robinson’s results, when cor- 
rected for the present value of e, yield 


Ss a, /A,) = (3.828 + -002)-10%4 (esu)?/g ergs (520 ppm) . 
mh 


Ropinson did not give an error for the final average, but quoted .075% as a limit of 
error in an earlier report. It seems safe to reduce this at least to .05% as given above; 
this is still about five times the statistical error of the result. 


3:3. Ionization and excitation potentials. — The so-called h/e measurements by 
means of measurements of ionization or excitation potentials are quite similar in prin- 
ciple to the photoelectric determinations. The energy necessary to ionize an atom or 
to excite a given level is determined in terms of wave-length 4 by spectroscopy mea- 
surements; this same energy is then measured by determining the voltage V necessary 
to accelerate an electron so that it will just ionize or excite the atom in question. It 
follows that he/A = 108Ve/e or h/e = 10°VA/c?. The measurement of ionization potentials. 
has been carried out by LAWRENCE [83] on Hg, and by VAN ATTA [84] and WHIDDINGTON 
and WooDpRoore [85] on He, Ne and A. DuNnNINGTON [15] has analyzed this data and 
obtained an average h/e of (1.3743 + .0035)-10-1’ erg s/esu (the value of e assumed 
is somewhat different from the present best value; but this is unimportant in view 
of the error). 

A much more accurate measurement, involving the excitation potential of He,. 
has been carried out by DunNINGTON, HEMENwAY, and RoucH [86]. The appratus. 

was quite similar to that used by DUNNING- 

bv TON [87] in his e/me experiment and is shown 
schematically in Fig. 3. A magnetic field was- 

FIRST SLIT applied perpendicular to the paper and was ad- 
Fore justed so that electrons emerging from the elec- 
tron gun had the proper radius of curvature to; 
GSELECTRON just pass through the first slit system. A small 
is amount of He was introduced in the gap bet- 
ween the slit systems. When the voltage bet- 
ween them was approximately zero (actually a 
Fig. 3. small voltage V, was necessary to correct for 

contact emf and similar effects), electrons which 

did not collide with the atoms passed through 

the second system and reached the collector; thus a current maximum was obtained’ 
for V,. When the voltage was raised by the amount of the excitation potential to 
a new value V, a second maximum occurred since electrons which lost enough energy 


SECOND SLIT 
SYSTEM 
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to excite one He atom by collision regained an equal amount in the gap and passed 
through the second slit system. The excitation potential measured was between the 
1s21S, (ground state) and the 2p!1P, levels. Dunnington’s results gave 


h/e = (1.3790 + .0002)-10-1 erg s/esu . 


‘The error given is principally a statistical standard deviation, but also includes an 
allowance for systematic errors in the electrical measurements. 

Additional sources of error should perhaps be considered. Dunnington has noted 
three possible sources of error in voltage measurements: contact potential, space charge, 
and surface charge. Since only the difference between V, and JV, is required, he eliminated 
any effect of contact potential as well as any surface or space charge phenomena which 
were independent of voltage. However, there is reason to think that the last two 
effects may have been voltage dependent. At voltage V,, most of the electrons reached 
the collector; the few which had one or more collisions were stopped at slit S. At vol- 
tage V,, those which had only one collision got through; all of the others were stopped 
by S. It seems quite plausible that this difference may have had some effect on space 
and surface charges in the second slit system. Of course, this effect may have been 
at least partially corrected for by an extrapolation to zero beam current which Dun- 
nington made. 

Dunnington also investigated other possible sources of systematic error. The ap- 
paratus was plated with gold at frequent intervals to reduce surface charge; he found 
some evidence ofa change of about 140 ppm in the excitation potential during the 
average time delay between plating and taking a run; however, the error in this cor- 
rection was such that he did not feel justified in including it. On the other hand, there 
was a small, but statistically significant, variation of ther esult with V,; DUNNINGTON 
made a correction of .0002 V to extrapolate from V, = .0047 V to V, = 0. Since there 
is no good explanation of this variation or even for the average value of V,, it does 
not seem obvious that one obtains the true result by correcting to V, = 0 as deter- 
mined by the potentiometer. A change of .01V in the value of V, used would change 
the answer by 20 ppm. In the light of all these considerations it seems safer to make 
at least a small allowance for additional systematic errors and treat the uncertainty 
given by DUNNINGTON as a probable error rather than standard deviation. Thus his 
result is 

h/e = (1.3790 + .0002)-10-1? erg s/esu (145 ppm). 


3-4. Excitation potentials for discrete X-ray lines. — A modified form of the method 
described in the last section can be used with X-ray spectra. When a metallic target 
is bombarded by high velocity electrons both a continuous and a discrete X-ray spec- 
trum are observed. The latter is, of course, characteristic of the atoms of the target. 
‘The so-called K lines are produced when one of the two electrons in the innermost 
level is removed and another electron falls into its place with the emission of a photon. 
Clearly it is impossible to excite this radiation unless the bombarding electron has at 
least enough energy to raise a K electron to the lowest unoccupied level of the target. 
This energy is equal to 108Ve/c, where V is the minimum target voltage necessary to 
excite K lines. 

The wave-length corresponding to this voltage may be determined from X-ray 
absorption measurements. When the frequency v is raised to a value such that an 
X-ray photon can eject a K electron to the lowest unoccupied level, a sudden rise in 
the absorption coefficient of the material is noted; this « absorption edge » represents 
the desired wave-length. (Since neither this absorption edge nor the excitation voltage 
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is defined with perfect sharpness there is some uncertainty in deciding exactly what 
points on the two curves correspond to each other). The photon energy is then he/A;. 
by equating this to the corresponding electron energy it follows that 


(hje)(A,/A,) = 108A’V /c? , 


where 4’ denotes wave-length on the Siegbahn scale. 

Some preliminary results by this method were briefly reported by BEARDEN and 
Scuwakz [88]. It is interesting to note that these results, which they then felt were: 
less reliable than their measurements of the continuous X-ray spectrum, now appear 
to be their best value obtained at the time. 

More recently NILSson [89] completed an extended series of measurements of 
K excitation voltages for Fe, Co, Ni, and Cu. The absorption data needed was taken 
from BEEMAN and FRIEDMAN [90]. NIiLsson gives a limit of error of about 450 ppm. 
The systematic errors in voltage measurement are combined as the square root of the 
Squares, giving 230 ppm; other errors are added arithmetically to this. It seems safe: 
to halve the voltage uncertainty in computing a probable error. Among the other 
errors, that in c? may be considered negligible while that in (A,/A,) does not appear 
in the present statement of the result. The uncertainty in locating corresponding points. 
on excitation and absorption curves is given as 50 ppm; this estimate should probably 
be retained without change in the absence of a good theoretical explanation of the 
shape of the curves. When all remaining errors are combined by taking the square 
root of the sum of the squares, the probable error is about 138 ppm. The original 
data, as given by NiLsson in Table XXIII, has been recalculated with the present. 
value of c; the result may be stated as 


(h/e)(A,/A,) = (1.3766 + .0002)-10-17 erg s/esu (145 ppm), 


where the probable error can be considered rather conservative. 


3°5. High-frequency limit of continuous X-ray spectrum. — The X-ray 
« h/e experiments » are essentially the inverse of the photoelectric determinations. 
A beam of electrons is accelerated by a known voltage and allowed to strike 
a metal target (the anode), where X-rays are produced. The maximum energy 
(i.e., the highest frequency) X-ray photon will be emitted when all of the 
energy of an electron is converted to radiation in a single collision. Thus is 
i' is the wave-length of the high frequency limit on the Siegbahn scale, the 
resulting equation is 


(h/e)(As/Ag) = 1082'V /e? 


the same expression as was given in the previous section. 

The most common procedure consists in setting an X-ray monochrometer 
for a given wave-length and observing the intensity as the voltage is varied 
about a small range near cut-off frequency; the resulting curve is called an 
isochromat. 
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Fig. 4 shows a typical isochromat for tungsten. The intensity is usually 
measured by a Geiger counter and must be corrected for the background count. 
The experimental curve does not, of 
course, show a perfectly sharp cut-off, 
but approaches the axis assymptotically. 
The usual way of locating the high fre- 
quency limit has been to take the point 
where the second derivative of isochro- 
mat is a maximum; this method was first 
proposed by DuMonD and BoLLMAN [91]. 
There is some experimental error involv- Fig. 4. 
ed, which depends on the number of 
counts for each point on the curve, as well as a possible theoretical error 
which will be discussed later. 

Other possible sources of error are impurity films on the target, thermal 
emfs, ripple voltage, thermal velocities of electrons, and cathode work function. 
It has been shown by DuMonD and BoLLMAN [91] that an energy equal to the 
work function is supplied to each electron from the filament power supply; 
the work function voltage must therefore be added to the plate voltage. The 
correction is known quite accurately for tungsten, but varies in a rather un- 
predictable manner for an oxide cathode. Most of the errors just mentioned 
involve constant voltages and can be reduced by going to a high voltage. 

The first precision experiment of this type was that of DUANE, PALMER 
and YEH[92]. However, the principal weight must now be given to three 
recent determinations (the figures quoted are the values of 108A’ V in cm emu): 


W TARGET 


INTENSITY ——> 


Vv—> 


a) In 1950 BEARDEN and ScHWARZ [93] obtained 12370.76 + .50 as a 
mean of 9 determinations with Cu, Ni, Mo, Ta, Au, and W targets. (In com- 
puting the above error, the uncertainty in the work function has been lowered 
to 30 ppm.) Most of the work was done with an indirectly heated oxide cathode 
at 8 to 10 kV, although one run was taken with a tungsten filament and two 
at about 20kV. A double crystal spectrometer was used. 


b) In 1951 BEARDEN, JOHNSON, and WATTS [94] performed a rather 
similar experiment with one major improvement; they used a tungsten fila- 
ment which was heated during alternate half cycles of the 60 cycle supply 
and employed a « gating circuit » to record counts only during the half cycles 
when the filament power was off. Thus they obtained essentially a unipo- 
tential cathode with an accurately known work function. The target was also | 
tungsten and runs were taken at 6 and 10 kV with wave-length as the variable. 
A third result was obtained by considering the difference in the high frequency 
limit for these two voltages. The mean of their values was 12371.09 + .50. 
(The published result includes the BEARDEN and SCHWARZ value in the average 
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as stated in Fig. 3 of the BJW paper. This seems to have caused some con- 
fusion. DuMoNnD and COHEN [1] used this final average and then introduced 
the B and S value as a separate item, thus effectively including the latter 
twice in their average.) 


e) In 1953 DuMonD, FELT and HARRIS[95] determined the high fre- 
quency limit at 24.5 kV for a tungsten target, by means of a curved crystal 
spectrometer. A tungsten filament was used with the same « gating circuit » 
arrangement employed by BEARDEN, JOHNSON and WATTS. The result was 
12370.02 + .42. 


If these three determinations are given equal weight, the mean value is 
(12370.63 + .28) cm emu. As will be shown later in this report, this average is 
significantly lower than that implied by other experimental data. The discre- 
pancy is also clearly shown in the DUMonD and CoHEN [1] evaluation of atomic 
constants though they do not consider the evidence strong enough to reject 
these X-ray experiments. 

There are at least two causes which may account for this disagreement. 
The first is the thermal energies of the electrons; the average value is 2 kT 
and the most probable value is 47 (see, for example, MILLMAN and SEELEY [96]). 
It can be shown that the correction should lie somewhere between these limits. 
For the oxide cathode used by BEARDEN and ScHwARZ kT should be of the 
order of .1 eV while for the tungsten filament used in the other experiments 
it is about .2eV. In any case it will tend to raise the experimental value. 

A more important correction may perhaps be required due to the method 
of determining the high frequency limit by the maximum second derivative [97]. 
In introducing this method, DuMoND and BOLLMAN [91] assumed that the 
ideal isochromat ccnsists of a series of straight 
lines and intersects the axis with a finite slope, 
as shown by the solid line in Fig. 5. 

However, they also pointed out that, if a 
discontinuity exists at the high frequency li- 
mit, as shown by the dotted line, then the cor- 
rect value is obtained by taking the maximum 
first derivative. They actually employed both 
methods and found a difference of about two 
volts; the observed isochromat could be re- 
constructed almost equally well on either hypothesis. For their final average 
a mean of the two values was used. While the difference between the two 
methods was not too important at the time, it becomes more critical with 
the accurate data now available. However, the second derivative method 
seems to have been accepted without much question in all recent work. 


INTENSITY —- 
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The assumption that the ideal isochromat consists of a series of straight 
lines with small discontinuities in the slope (the «knees» referred to by 
DuMonD and BOLLMAN [91]}) is clearly no longer tenable. The recent work. 
on tungsten by BEARDEN and ScHWARZ[93] and BEARDEN, JOHNSON, and 
WATTS [94] shows, at 3 to 4 V above what is 
taken as the cut-off voltage, a peak followed 
by a distinct minimum; this was first noted 
by OHLIN [98]. NITBOER [99] has suggested there 
may be a range of forbidden energies between 
the Fermi levels of the target and that the 
intensity of an ideal thin-target isochromat 
may be zero in this region, as shown in Fig. 6. I fe 
DuMonp [100] also considered the problem and dot 
reached substantially the same conclusions. 

Thus the changes in slope near the cut-off voltage may well be of major im- 
portance. 

The original DuMond and Bollman theory postulated that the ideal thin- 
target isochromat had a discontinuity at the cut-off voltage and was constant 

for higher values, as shown in Fig. 7(a); 

they cited SOMMERFELD’s [101] work as 

theoretical justification. The thick-target 

curve was then considered as formed by 

_ the superposition of many thin-target ones 

with the electrons losing some energy in 

each layer. The resultant curve was a 

Fig. 7. straight line intersecting the axis, which 

is indicated in Fig. 7(b). This analysis, 

as well as the later treatment by DuMonp [100], assumed that the electrons 

. lose energy continuously and that the energy loss at a given depth is the same 
for each electron. 

However, there is considerable evidence that, energy losses occur largely 
in discrete steps of the order of 10 eV. Furthermore, even at depths where 
the average loss is much larger than this, there appear to be an appreciable 
number of electrons which have not undergone any measurable loss in energy. 
(See, for example, the work of RuDBERG [102], RUTHEMANN [103], LANG [104], 
and MARTON and LEDER[105]). If these statements are correct, then a dis- 
«continuity should be expected at the cut-off voltage even for a thick target; 
the magnitude of this discontinuity is difficult to fix. GERTHSEN and Ar- 
BERT [106] have developed a theory from a simple model with discrete energy 
losses and obtained an isochromat with a series of sharp steps. However, 
their work gave no explanation of the peaks and valleys which are observed 
and are qualitatively explained by  Nijboer’s model. 
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To form some estimate of how much the experimental results may be misinter- 
preted, a simplified example will now be analyzed. Assume that for any constant vol- 
tage V, with corresponding cut-off frequency ») = 108eV/ch, the characteristic curve: 
of intensity vs. frequency is 


B+C(m—?), v <> 
Lv; 5) = 


0. Vv > Vg 


This is indicated in Fig. 8; the intercept of the project- 
ion of the characteristic is given by d=v —wm = 


V = CONSTANT : 

Fe = B/C. B and C are considered constant over a. 
2 he small range of V. Assume further that the ideal iso- 
» o}———=-—-—= x chromat is blurred only by the limited resolving 
= x 


power of the spectrometer and that this can be re- 
% presented by the Gaussian function 


Fig. 8; ie, ») La exp |- ( | 1 


where »; is the frequency setting of the spectrometer. The observed intensity for given’ 
values of V (i.e., v9) and », is now found by integrating f(»,, v)I(v, v9) between the limits. 
of 0 and »; there is, of course, no radiation above »,. The resulting function F(v,, v9) 
s thus 


Vo 


v—v;\? 
FO x) = fa exp|-( 2 J] + oro. 
0 


The equation for an isochromat is. obtained by holding v; constant and varying wp. 

The maximum of the second derivative of the isochromat is found by differentiating 
F(vs, v9) three times with respect to » (which is directly proportional to V) and setting 
the third derivative equal to zero. The result is 


Vo 
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The third derivative will vanish when 
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If B = 0, it follows that n —v = 0 and the maximum of the second derivative gives 
the true cut-off frequency. On the other hand, if C = 0, the above equation gives 
%—v; = + 0/2; the negative root corresponds to the maximum of the second deriv- 
ative. The first or second of the above cases will be approximated if o is very large 
or very small respectively compared with the quantity ò in Fig. 8. In general, the 


desired root is 
vv 1 fo0 o0\2 + 
= —_| +2|}. 
o 2|2B |\2B 


From this equation it follows that, for positive C 


1 LIT 


V2 eo 


Thus, under the conditions assumed, the second derivative criterion can misinterpret 
the cut-off frequency by as much as 0//2. It may easily be shown that the half width 
w of the Gaussian curve of the spectrometer is approximately .830; hence, roughly 
speaking, the second derivative method may misinterpret the result by w. Since », <,, 
it follows that the voltage obtained by the second derivative method (i.e., that cor- 
responding to ») will be less than or equal to the true voltage for the setting »,. The- 
refore, if any error exists in this method, the calculated value of h/e will be lower 
than the true value. 


Since precise structure of the ideal isochromat is not known, it seems unsafe 
to assume either that such a correction should be applied or that it should 
not: but it appears probable that the true value lies somewhere in this range. 
Thus it seems reasonable to apply a correction of 2k7'+-(hw/2) and increase 
the error appropriately. The details of these corrections are shown in Table V; 
the resulting unweighted mean of the three experiments is 12371.7. It will 
be noted that the BEARDEN, JOHNSON, and WATTS results are much more 
consistent after corrections. 

It is not obvious what error should be assigned to the result. The cor- 
rection factor is different for each experiment. On the other hand, the errors 
are not completely independent, since all involve the structure of the ideal iso- 
chromat (in most cases for a tungsten target). The correction to the DUMoND 
result is 139 ppm; the correction to BEARDEN and SCHWARZ and BEARDEN 
JOHNSON, and WATTS average 46 ppm and 59 ppm respectively. The average 
correction for the three experiments is thus 81 ppm; it seems reasonable to 
assign a probable error of 60 ppm due to the corrections. Combining this with 
the probable error of the mean due to other causes (23 ppm) gives a total error 
of 64 ppm. When the result stated in the preceding paragraph is divided 
by c?, the final value now becomes 


(h/e)(A,/4y) = (1.37653 + .00009)-10-1 erg s/esu (65 ppm). 
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If the first derivative criterion is applied, the change in the experimental 
values is in general considerably greater. In most cases the shift is too large 
to be plausible since the observed intensity should be essentially zero when 
the wave-length setting of the spectrometer is below cut-off by several half- 
widths. The shape of the curves near the maximum of the first derivative is 
probably influenced by the same mechanism producing the peaks and valleys. 
It seems unlikely that the present curves can be interpreted with high accuracy 
until a better theory is found to explain the form of the isochromats. 


TABLE V. — Corrected values of « h/e experiments ». 
3 he?/e)(As/A 
i SI Voltage n. kT w 2kT wr: és da i ra 
xperimen (kV) UBS I (ev) |- (eV) +hw/2 rection 
(ppm) | Original Corrected 
BEARDEN 8 5 1 | .4(*)| .4 | + 50] 12370.96 | 12371.58 
and 9.9 2 sì .5(*)| .45 | + 46|12370.06 | 12370.63 
SCHWARZ 19.6 2 a! 1.0 (*) “yl + 36! 12370.96 | 12371.40 
Wt. average 12370.76 | 12371.33 
BEARDEN, 10.2 1 2 LEA 65 | + 64 | 12371.20 | 12371.99 
JOHNSON 6.1 1 “2 -3(*) | .55 | +. 90 | 12370.75 | 1237186 
and WATTS ASTE) ia — VARI + 24| 12371.92 | 12372.22 
Average 12372.02 
DuMoND, 
FeLT and 24.5 .| All 2 | 6.0 3.4 +139 | 12370.02 | 12371.74 
HARRIS 
Unweighted average of experiments 12 371.70 


(*) The spectrometer used in the Bearden experiments had a resolving power 4/44 ~ 10000. 
Thus the half-width w is about 1 eV for 10 kV and is approximately proportional to the voltage. 
It is assumed that other factors contributing to effective half-width are negligible. 

(+) This value was obtained by taking the difference between the 6.1 kV and 10.2 kV 
results; it therefore involves no error due to thermal energies. In the original report of the 
experiments it was considered less accurate than the other two determinations, but in the light 
of the various uncertainties cited it is now given the same weight as the others. 


3:6. Electron diffraction. — The wave properties of the electron were first suggested 
by pe BroGLIE and soon afterwards confirmed in the experiment of DAvisson and 
GERMER [107]. At non-relativistic velocities, the wave-length is h/mv. This may be 
determined by diffraction experiments; Brage’s Law applies just as in the case of 
X-rays [108]. However, since the specimens used are usually not single crystals but 
rather thin films formed by evaporation, the technique is similar to X-ray powder 
photographs. For a given Bragg angle 6, there will be some microcrystals oriented 
so that the beam of electrons strikes them at this angle of incidence. The contributions 
of various layers will then reinforce one another and strong direct reflection will take 
place; since the various microcrystals have random azimuth angles, the result will be a 
cone of rays making an angle of 20 with the incident beam. For each Bragg plane there 
will be one cone corresponding to each order of the Bragg equation. Thus the wave- 
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length is obtained in terms of the measured angle and the grating space of the crystal, 
usually known in Siegbahn units. 

The velocity of the incident electrons may be determined in two ways. The first 
is by means of direct measurement with the technique used by KIRCHNER [54] in his 
e/me measurement. Such experiments have been performed by MEIBON and Rupr [109] 
and Gwnav [110]. They yield essentially a value of (h/m)(A,/A,); but the accuracy is not 
high, the error being of the order of .1%. The second method is to accelerate electrons 
by means of a known voltage V; in this case the non-relativistic energy is 108 Ve/e and 
the corresponding velocity v = (2-108 Ve/mc)?. Combining this expression with the 
wave-length equation gives (A,/A,)h/(em)* = 4'(2-108V/c)3, where A’ is measured on the 
Siegbahn scale. When relativistic corrections to momentum and kinetic energy are 


included, this becomes 
h 3 V 108 Ve\ |? 
SSL em a a a 1 ; 
temi | IAs) = a me ) 


The first precision determination of this type was made by von FRIESEN [111]. 

A more accurate version of this experiment was recently performed by RYMER 
and WRIGHT [112]. They used a modified form of an electron microscope at 50 kV. The 
voltage was determined by a precision voltage divider, which was calibrated under 
full load. Four different substances were used to produce diffraction, viz., Al, NaCl, 
and CsI. The grating constant for NaCl [113] is taken from Bergen’s [63] work. The 
constant for CsI (and possibly that for KCl which is based on unpublished work) is 
obtained from powder photographs and may be somewhat less accurate for this reason. 
The estimate of uncertainty in grating constants, 20 ppm, may perhaps be low; but 
the total error appears reasonable. It is stated in terms of standard deviation and cor- 
responds to a probable error of about 67 ppm. Their result then becomes 


cal (4/3) == (.999 54 + .00007)-10-8 erg s/(esu g)? (70 ppm). 

3:7. Compton effect. - When a high energy photon is scattered by a free electron 
or by one lightly bound to an atom, the frequency of the scattered photon is less than 
that of the incident one. In the case of a free electron the theoretical explanation can 
be obtained by treating the photon as a particle with energy hy and momentum hy/e 
and using the principles of conservation of energy and conservation of momentum 
(see, for example, RicutmMpyER and KENNARD [114] or Compron and ALLISON [115]). 
The resulting change of wave-length is found to be AA = (h/me)(1 — cos g) where p is 
the angle between incident and scattered radiation. 

Actually, while atomic electrons behave almost as free particles for high frequency 
photons, some correction is required to take account of the binding forces. Ross and 
KIRKPATRICK [116] developed a simplified theory indicating that this correction was pro- 
portional to the square of the wave-length for a given scattering material. They then per- 
formed an experiment using three different wave-lengths with carbon as scatterer. By 
extrapolating their results to zero wave-length they obtained h/me = (.02418 +.00004) A. 
As far as can be determined from their paper and a recalculation of some of their data, 
wave-lengths were measured on the Siegbahn scale and the conversion factor 4,/4s 
was not applied. If this is correct, the result should be expressed as 


I 
2 _(1,/2;) = (2.418 + .004)-10-:° em (1650 ppm). 
me 
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3:8. Compton wave-length by annihilation radiation. — When a radioactive sub- 
stance emits a positron, the latter will eventually encounter an electron; both will 
then be annihilated with the emission of two photons. The probability of this process 
occurring is very small until the positron has slowed down to thermal velocity; hence 
the energy of the photons, i.e., 2hy must come almost entirely from the rest mass 
of the two particles, which is 2mc?. (The potential energy cancels out since the po- 
tential energies of the positive and negative particles must be equal and opposite 
when they collide). If 2, denotes the wave-length of the photons on the Siegbahn 
scale, the frequency v is equal to (c/A,)(A,/A,). Substituting this value for v and 
equating the two energies yields 


h ' 
ANA 
me 


This is known as the Compton wave-length, which was discussed in the previous section. 

DuMonp, Linp and Watson ([117] measured this wave-length by means of a curved 
crystal spectrometer and positrons emitted from a ®Cu source. The value was slightly 
higher than the best indirect value; DuMonp interpreted this as possible evidence 
for a difference between the mass of the electron and that of the positron, the latter 
being about 80 ppm lighter. Another measurement by Linp and HEDGRAN [118], com- 
paring annihilation radiation with y-rays emitted by !*8Au, seemed to confirm this 
view. However, DuMonp further improved the accuracy of his spectrometer and 
discovered small mechanical nonlinearities which had affected his previous results. 
A newer measurement by MuLLER, Hoyt, KLEIN and DuMonp [119] gave the value 


h 
— (A,/A,) = (2.4213 + .0003)-10-1° cm (124 ppm). 
me 


They also recalculated the Lind and Hedgran result, based on their own measurements 
for !98Au y-rays, and found almost perfect agreement. 


3:9. Radiation constants. — The spectral distribution of energy radiated by a black 
body was theoretically determined by PLANCK [120]. (See RicatmeYveR and KEN- 
NARD [121] for a good explanation of the theory, but note that most of the results are 
stated in terms of energy density rather than energy radiated by the body). The 
resulting equation is 

Gq 
br 45(exp [¢,/AT] —1)’ 


C, = 2nc*h, Gsi=Gh/h 


where 4 is the wave-length, k is the Boltzmann constant, J,d/ is the energy radiated 
per unit area per unit time in the wave-length interval dA, and the other symbols 
have their usual meanings. €, and c, are known as the first and second radiation con- 
stants. By integration it is possible to derive the Stefan-Boltzman Law, which states 
that the total energy emitted by unit area per unit time is 074. Differentiation of J, 
gives Wien’s displacement law, i.e., Z,,7 is a constant (where 7, is the wave-length 
for which J; is a maximum). The values of these constants are thus shown to be 


o = (27°/15)(k4/e2h8), 
Aml = 6,/4.96511 = he/4.96511%k. 
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Various methods of measuring these constants have been described in a review 
by WENSEL [122] and more recently in a shorter paper by RuTGERS [123]. The con- 
stant o may be measured by constructing a cavity to approximate a black body and 
measuring the total radiation at a known temperature with a thermocouple or bolo- 
meter. WENSEL lists 20 measurements and takes the unweighted mean; RuTGERS 
favors omitting the first four, which involve work previous to 1910. If this is done 
a more consistent set is obtained and the result is o = (5.74 + .02)+10-5 erg/cm? s °C4, 
with the probable error computed by external consistency. Obviously the accuracy 
is extremely low as compared with other measurements of atomic constants. 

The constant e, is determined principally by two methods: 1) The constant 7,7 
may be found and ¢, computed by the exact numerical relationship quoted above. 
WENSEL gives a table quoting seven values; the last five (those obtained since 1910) 
will be used in the present calculation, as in the Rutgers paper. 2) An optical pyrometer 
may be used to find the ratio of radiant energy flux at two different temperatures for 
a constant wave-length or at two different wave-lengths for a constant temperature. 
In this way the constant ¢, cancels out of the resulting equation and e, is obtained. 
A total of six such experiments by various workers are quoted by WENSEL and RUTGERS. 
The unweighted mean of these eleven values gives c, = (1.434 + .001) cm °C. Un- 
doubtedly a careful analysis and proper weighting would give a more accurate result, 
but the accuracy of the data does not seem to justify it. The value of e, has also been 
determined by RurGERS [123] in photometric measurements giving the energy equivalent 
of luminous flux; the results are slightly lower than the mean quoted above. 

In discussing these various results, both HOFFMAN [124] and pr Groot [125] point 
out that they all depend on the accurate knowledge of the absolute temperature of 
the melting point of gold, which is used as a reference for high temperature work. This 
is now being redetermined carefully for the first time since 1911; Moser [126] has reported 
his preliminary value is in essential agreement with the accepted figure. However, 
it seems best to allow for a possible error of at least 1 °C common to the above measu- 
rements. When the resulting error is then calculated, the above results become 


o = (5.74 + .03)-10-5 erg/cm? 8 °C4 (5250 ppm), 
Cy = (1.434 + .002) em °C (1400 ppm). 
3-10. Index of refraction for X-rays. — According to the theory of X-ray dispersion 
(see Compton and ALLISON [127] and WHEELER and BEARDEN [128]), the index of re- 


fraction x of a material (or the decrement 6), defined as «1 can be computed by 
the equation 


e212No 
> (1+ As) . 


ae 2ame*?M *; 
Here o is the density of the material, M its molecular weight, A the wave-length of 
the X-rays, n, the effective number of electrons in the st* level and A, a factor which 
takes into account the electronic binding and the X-ray absorption for the s™ level. 
If the wave-length 7’ is measured on the Siegbahn scale, the equation may be rearranged 
to give 

Ne? 275M 


— (A,/A,)? = 
eat als) Ko 


[Dn +4). 


BEARDEN [129] employed this relation, together with the value of the faraday, to 
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obtain e/mc; he used a diamond prism since this crystal had low absorption and the: 
factors n, and A, could be evaluated quite accurately. 6 was measured for the Cu-Kg line 
and the density determined by placing the diamond in a special solution whieh could. 
be adjusted so that the diamond would neither float nor sink. With the use of the 
wave-length on the Siegbahn scale quoted by BEARDEN [58], the result is 1.7056 -+ .0003. 
BETHE and LONGMIRE [130] have pointed out that, when the scattering of X-rays by 
nucleii is included, the result is lowered by 2.7 parts in 104. The revised figure is thus- 


Ne? 
na (45/44)? = (1.7052 + .0003)-101 (esu)?/erg mole (175 ppm). 
to 


The principal source of error remaining may lie in the approximations used in the 
dispersion theory. 


3:11. Zeeman effect. — The only spectroscopic measurement which can now be 
considered a determination of e/me is that of the Zeeman effect. The energy levels of 
an atom in a magnetic field are shifted by an amount Mgu,H, where M is the magnetic: 
quantum number, yy = eh/4rme is the Bohr magneton, and H is the magnetic field 
strength and g is the so-called Lande g-factor which is unity (or very nearly so) for atoms: 
showing the normal Zeeman effect, i.e., having total spin zero. When the shift is ex- 
pressed in wave numbers and when g is unity, the expression beemes Av = MHe/42me®.. 
The selection rules for transitions require that AM = 0 or + 1, which means that in 
the normal Zeeman effect a single spectral line splits into three components in a mag- 
netic field (see HERZBERG [131] for further details). By measuring the splitting and. 
the field strength, it is thus possible to determine e/mc?. Since c is known to high 
accuracy, it is possible to convert this to e/m or e/me. 

KINSLER and Houston [132] determined such a value by measuring the Zeeman effect 
for Zn, Cd, He and Ne, All of these have singlet states and give the normal Zeeman effect... 
By observing the emitted light in the direction of the field, only two components were ob- 
served (AM= +1); these wave-lengths were measured by a Fabry-Perot interferometer. 
The mean of these measurements was given as (1.7570 + .0007)-107 emu/g. This error 
is the mean deviation, about twice the statistical probable error. However, it seems. 
best to make some allowance for possible systematic error and to let the original 
figure stand. It is not entirely clear just what auxiliary constants were used, but it 
would appear that the result should be raised about 100 ppm to correct for a change 
in the conversion factor from international to absolute electrical units. The result 
then becomes 


e/me = (1.7572 + .0007)-107 emu/g (400 ppm). 


3°12.* Rydberg constant and the atomic mass of electron. — The Bohr theory 
gave energy levels (in wave numbers) for hydrogen and hydrogen-like atoms. 
(e.g. Het) as Z?R,/n2(1+m/m,) where n is the principal quantum number, 
Z is the atomic number, R_ = 22?me4/ch® is the Rydberg constant for infinite 
nuclear mass, and 1+(m/m,), is a reduced mass correction factor depending 
on m/m,, the ratio of the mass of the electron to that of the nucleus. Z and n, 
being integers, are precisely known; the energy levels can be determined to: 
high accuracy by spectroscopic measurements of wave-length. If measurements. 
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are available from two or more such spectra for atoms of known mass, it is. 
possible to calculate both R, and m. R, can be determined to extremely 
high accuracy, while m will have a much higher error since it is determined 
by the difference between two quantities which are nearly equal. If the known 
masses are given in a.m.u., then the quantity directly determined is M,= Nm, 
i.e., the atomic mass of the electron. 

Actually, of course, the simple levels given by the Bohr theory do not 
exist. It has long been known that a given n level is split into levels cor- 
responding to different j values. This effect is due to relativistic and spin cor- 
rections, and is given by the Dirac theory (see, for example, RoJANSKY [133]). 
Recently a further splitting has been discovered for levels with the same j, 
but different 1. This is the so-called Lamb shift and consists principally of 
radiative corrections given by quantum electrodynamics. COHEN [134] has. 
recently analyzed spectroscopic data on He, D, and H in the light of all these 
corrections and obtained for the Rydberg constant and the atomic mass of 
the electron 


R_ = (109737.31 + .01) em .(.1 ppm) 
M, = Nm = (54.895 + .008)-10- amu (145 ppm). 


(One figure has ben rounded off in R,). The error in R, could be ten 
times as great as that stated and still not appreciably effect any of the 
constants to be derived from it. Thus it may be considered ag precisely 
known when combined with other data. For this reason also, it turns out 
that the Lamb shift corrections to R, have no appreciable effect on the 
derived constants; hence they will not be discussed in any detail in this 
section. 

It should be noted that the measurement of M, = Nm mentioned above 
has frequently been referred to as a spectroscopic value of e/me. It becomes 
so only when combined with the faraday Y= Ne/c. Since there is consider- 
able uncertainty as to the best value of F and since it is desirable to avoid 
correlations between various experiments so far as possible, it is here consi- 
dered as a value of Nm. 


3°13.* Gyromagnetic ratio of the proton. — The gyromagnetic ratio of the 
proton is defined as the ratio of its magnetic moment to its angular momentum. 
This is of interest in itself and is also valuable in furnishing a means of ca- 
librating magnetic fields. It can be measured by determining the nuclear 
resonance frequency of the proton in a known field. The proton has spin #/2, 
i.e., its projection in the direction of an external magnetic field is either #/2 
or —h/2. The corresponding values of energy are + ,H where mw, is the 
magnetic moment of the proton. The frequency »v, of a transition between 
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these two states is easily found from the equation hy, = 2u,H (+, is just the 
classical Larmor precession frequency). It follows that the gyromagnetic ratio 
Y, is given by 


A good general survey of the theory and techniques of nuclear resonance has 
been given by PAKE [135]. Actually in experimental work the protons are not 
free particles, but are usually contained in samples of water or mineral oil. 
A small diamagnetic correction must then be made to account for shielding 
by the electrons. Since other experiments use a similar proton resonance 
device, it is convenient to make use of the directly measured quantities without 
diamagnetic correction, i.e., the apparent magnetic moment TÀ and the ap- 
parent gyromagnetie ratio Vr 

This ratio has been measured very accurately by THoMAS, DRISCOLL, and 
HIPPLE [136] at the National Bureau of Standards. The proton sample was 
placed in a small inductance coil about lin. in diameter, which formed part 
of a resonant circuit fed by a 20 MHz oscillator. When the magnetic field was 
adjusted to produce proton resonance at 20 MHz, the nuclear absorption of 
energy lowered the Q and thus changed the voltage across the coil. The field 
strength was approximately 4700 gauss; it was produced by a large electro- 
magnet with a 2in. air gap about a foot in diameter. A pair of small Helmholtz 
coils immediately adjacent to the proton sample provided the fine adjustment 
to obtain resonance. The field of the electromagnet was measured by deter- 
mining the force on a 10-turn rectangular coil carrying a known current. 
Almost all the force was on a 10 em length on the lower end of the coil; but 
«dimensions were measured very accurately and the necessary corrections made. 
The actual field was not quite uniform; it was carefully mapped by means 
of the proton resonance probe, which gave the ratio of the field at any point 
to that at a chosen reference position. 

The error given is 22 ppm, which is intended as a limit of error. While 
numerous sources of error were considered, there were four principal ones, 
each estimated at 9 or 10 ppm, viz., precision of the balance, electrical standards, 
field distribution, and width of coil. The precision of the balance produces 
a statistical error; the last two errors may also be in part statistical. The 
statistical probable error is only 2 ppm. The error in the ampere has been 
previously discussed [137] and estimated at 6 ppm. Field distribution may be 
the most serious problem; relative field strengths are stated to at most five 
significant figures; therefore it would seem that the probable error due to this 
cause cannot be much less than 10 ppm. On the whole, it would thus seem 
that the figure of 22 ppm is fairly conservative, but that the probable error 
should be appreciably greater than 10 ppm. It will therefore be taken as 
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15 ppm and the result, without diamagnetic correction, stated as 


y, = 2y,/fi = (2.67523 + .00004)-104 s-! oersted-! (15 ppm). 


In regard to the diamagnetic correction, there are two points to be considered. 
‘The first question is how much the correction factor varies between different proton 
samples. GutTowsky and McCLURE [138] found that the correction factor shift com- 
pared to H, gas as standard was (.3 + .3) ppm for H,O and (3.7 + .3) ppm for mineral 
oil. THomas [139] obtained the values —.6 ppm and + 1.6 ppm for water and oil 
respectively. A small amount of ferric salt was added to the water in the Thomas, 
Driscoll and Hipple experiment to reduce the relaxation time; this was compared with 
:a standard oil sample used as a reference. Thus it is clear that the differences between 
various samples are small, a few ppm at most. The second problem is the absolute 
diamagnetic correction in H, gas; this is needed only in calculating the absolute values 
‘of the magnetic moment and gyromagnetic ratio. This correction has been calculated 
as + 26.8 ppm by Ramsey [140], + 26.6 ppm by NEweELL [141], and + 29.5 ppm 
(stated as an upper limit) by HyLLERAAS and SKAVLEM [142]. It seems reasonable to 
adopt the value of + 27 ppm. With the inclusion of the Thomas correction between 
the standard sample and H,, the final correction becomes approximately + 28 ppm. 
This appears to be correct within about 1 ppm, which is completely negligible compared 
with the original experimental error. The corrected value now becomes 


Yo = (2.67530 + .00004)-104s-1 oersted- (15 ppm). 


3°14.* Cyclotron resonance of the proton (proton magnetic moment in nuclear 
magnetons). — A proton in a uniform magnetic field will revolve in a circle 
with angular velocity ©, = He/m,c; this is known as the cyclotron resonance 
frequency of the proton. If the proton resonance frequency described in the 
last section (w, = 27, = 4au,H/h) is measured in the same field, the ratio is 


n 
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This is the magnetic moment of the proton expressed in nuclear magnetons. 
The same diamagnetic correction is involved as before so that the directly 
measured frequency ratio involves TÀ rather than y,. 

SOMMER, THOMAS, and HIPPLE [143] at the National Bureau of Standards 
devised an ingenious experiment for measuring this ratio. The electromagnet 
was the same as that used in the THOMAS, DRISCOLL, and HIPPLE experiment, 
and the technique of measuring proton resonance was similar. The frequency », 
was measured by means of a device called the « omegatron », which was some- 
what like a small cyclotron and is shown in Fig. 9. It employed a uniform 
magnetic field, perpendicular to which was an electric field between two con- 
denser plates. A gas pressure of 10-* mm or less was maintained within the 
omegatron; protons or other positive ions were produced near the center by 
electron beam bombardment. A number of guard rings surrounded the ome- 
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gatron on four sides; these were given a positive d.c. potential (trapping vol- 
tage) of about .1 V to prevent protons from escaping in the axial. direction, 
i.e., the direction of the magnetic field. The electric field was supplied by 
an oscillator with a frequency of approximately 7 MHz. When the frequency 
was the same as the cyclotron frequency y,, then any protons which started 
in proper phase were continuously accelerated since the direction of the field 
changed each time they completed half a revolution. (The principle is the 
same as that of the cyclotron except that the particles gain energy continuously 
rather than in discrete steps when crossing a gap.) 


PATH OF ION AT R.F. OSCILLATOR 
RE SONANCE- 


ION ION COLLECTOR 
R.F PLATE COLLECTOR TRAPPING 
È =, VOLTAGE 
CUT-AWAY VIEW (SCHEMATIC) SIOE VIEW(MAGNETIC FIELO OUT OF PAPER) 
(Fig. 9. 


RE If the magnetic field was varied slightly, a proton might still gain energy 
from the field for a considerable number of revolutions, but would eventually 
become out of phase and decelerate. Thus, for any field strength such that 
the applied frequency was not quite equal to v,, there was a certain maximum 
radius attained by the protons. An electrode was located near one of the 
plates to collect the ions which reached this radius, about 1em. As the mag- 
netic field was varied through the resonance value, a-curve was obtained with 
half-width in the order of one part in 35000. A small correction to the 
measured v, was needed for the effect of space charge and trapping voltage. 
This was made by carrying out measurements on H*, H#, and D; and extra- 
polating to zero M, as indicated by a theoretical analysis. The theory was. 
well supported by the consistency of the final results. The final value without 
diamagnetic correction was v//v, = 2.792685. 

The error stated is 21 ppm. This includes an allowance of 5 ppm for sta- 
tistical error (5 times the probable error) and the rest for possible systematic 
errors. The latter are believed to be due to three causes: (1) magnetic effect 
of the filament current, (2) magnetic shielding of the electrodes, and (3) asym- 
metry of the resonance peak. The first effect was checked by reversing filament 
current and the second investigated by placing the proton sample inside the 
omegatron. No definite change was observed in either case; such errors must 
have been in the order of 1 ppm or less. The unexplained asymmetry in the 
resonance peak amounted to a few ppm; it may have been in part statistical. 
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In light of these considerations it seems justified to divide the stated error 
by three and use 7 ppm. Thus the result becomes 


poli, = 2.79268, + .00002 (7 ppm) 
Ml, = 2.79276 + .00002 (7 ppm). 


The same quantity has been measured by JEFFRIES [144], using an inverse or dece- 
lerating cyclotron, which is shown schematically in Fig. 10. A narrow beam of protons 
was accelerated to an energy of about 20 kV 
and injected into the cyclotron, which was 8.5 em ORBIT OF RESONANT 
in diameter and was placed in a magnetic field of PROTONS 
approximately 5300 G. The cyclotron resonance 
frequency », (vz in Jeffries’s notation) was about 


OSCILLATOR 


8 MHz for this field. When it was operated at 
precisely », protons starting in the proper phase i Sire COLLECTOR 
were decelerated each time they crossed the gap INJECTOR PLATE 


between the dees. For frequencies close to this, (MAGNETIC FIELD INTO PAPER) 
they could still be decelerated for many revolu- Fig. 10. 
tions, but eventually became out of phase with 
the voltage and were accelerated again; this was 
the principal factor determining the band width of the observed resonance curve. If 
the cyclotron was operated at some odd multiple of »,, say nv,, it was equally effec- 
tive in decelerating the protons. However, for frequencies slightly off ny, the timing 
was more critical and the relative band-width reduced by 1/n. In practice, most of 
the data were taken at the ninth harmonic. The frequency v, was obtained by inserting 
a proton resonance probe in the same field, just as in the Sommer, Thomas and. Hipple 
experiment. 

Experimental curves were obtained by varying the dee frequency v through a small 
range around the resonant frequency ny,; an approximate theory developed by Jeffries 
indicated that the curves thus obtained were not centered about the resonant frequency. 
The method of interpreting the results was essentially as follows. There existed a cri- 
tical dee voltage V, (about 200 V), below which no collector current could be obtained, 
since the protons had to lose a minimum amount of energy in the first revolution in 
order to clear the injector plate. At the voltage V, only those protons whose phase 
was such that the decelerating voltage was a maximum as they crossed the gap could 
reach the collector. Jeffries’s approximate theory showed that such protons would 
be defocused and strike the dees (i.e., move into or out of the paper in Fig. 10) if 
the dee frequency v was slightly above the nv,; for frequencies slightly below ny, focusing 
would occur. Hence the desired frequency nv, could be determined in principle by 
observing the high frequency cut-off of the curve when the dee voltage was equal to V,. 
This was difficult to observe. Therefore curves were taken at a number of higher 
voltages and the results extrapolated to V, by means of an approximate equation 
developed in the theory. The average obtained from many such determinations was 
a value of RITA equal to 2.79237, which was rounded off to 2.7924 + .0002. 

JEFFRIES considered various sources of error; the above result includes estimated 
corrections for known systematic errors (except diamagnetic correction). However, 
the one error which overshadowed all others was that due to the approximate theory 
referred to above. He estimated that this error was not greater than one-third of the 
width of the resonance curve at the highest harmonic used, i.e., the 9th harmonic. 
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This gave an error of about 70 ppm, apparently intended as a limit of error. However,. 
it seems safer to regard this figure as a probable error. The theory given by Jeffries- 
was in excellent agreement with experiment for lower harmonics; however, for the 
9th harmonic the band width predicted was in error by as much as a factor of three.. 
He attributed this discrepancy to weak focusing action at higher frequencies. In 
extrapolating to V, he used the observed band width rather than the calculated value 
in his theoretical equation. It is by no means certain that this one change is sufficient- 
to make his equation a valid basis for extrapolation when the original theory has clearly 
broken down. Furthermore it seems very probable that, while an exact theory might. 
place », inside the resonance curve for V= V,, it. would not place », totally outside 
this curve. Thus, if the extrapolation to cut-off voltage is carried out correctly, the 
residual error is more likely to be in one direction than the other. 

Comparison of the JEFFRIES experiment with that of SomMER, THOMAS and HIPPLE: 
shows that (1) the theoretical error in the Jeffries experiment is much more difficult 
to assess quantitatively than the experimental ones in the STH work, (2) there is some 
reason to believe the Jeffries error is more likely to make the apparent value of », too: 
high than too low (which might explain the fact that his value is lower than the STH [143] 
result), and (3), if weighted inversely as the square of the errors as assigned above, 
the JEFFRIES result carries only about 1% of the weight of the STH one. Consequently, 
while the JEFFRIES result is very helpful in essentially confirming the STH value by a 
different method, it seems best to omit the former entirely in calculation of atomic 
constants. 


3:15. Cyclotron resonance of the electron (proton magnetic moment in Bohr mag- 
netons). — If the proton nuclear resonance frequency is compared to the cyclotron 
frequency of the electron v,, the ratio is the same as that derived in the preceding section 
with m, replaced by m, i.e., 
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Hence the result is the magnetic moment of the proton expressed in Bohr magnetons.. 

This quantity was measured by GARDNER and PuRCELL [145] in an experiment 

shown schematically in Fig. 11. A narrow rib- 

bon of electrons was produced by an indirectly 

heated cathode and accelerated by a small vol- 

tage V,. A uniform field of about 3300 G was 

applied parallel to the beam. Due to small ther- 

mal velocity components perpendicular to the 

beam direction, the electrons moved in helical. 

orbits with radii of the order of 3um. The rota- 

va: ANALE tional frequency, i.e., the cyclotron frequency,. 

was approximately 9400 MHz, the electrons com- 

Fig. 11. pleting roughly 1000 revolutions between cathode 

and target. The magnetic field thus provided a 

focusing action and kept the beam from spreading 

appreciably between cathode and target under normal conditions. (It is interesting to 
compare this experiment with the Buscu [87] method of determining e/me). 

However, the section between cathode and target (about 2.3 cm long) was a part. 

of a wave guide, which was excited so as to produce a high frequency electric field. 
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perpendicular to the beam. When this frequency was equal to the cyclotron frequency 
of the electron, there was a gradual increase in the radius of the electron orbit, just 
as in the case of the omegatron described in the preceding section. Hence it was ex- 
pected that a current minimum would be obtained in this case. Such a minimum was 
indeed observed when the accelerating voltage V, was 5 V or greater. However, for 
values of V, less than 3 V the current was found to be limited chiefly by space charge. 
In this case the high frequency field tended primarily to remove the slowest electrons- 
from the center of the beam, thus decreasing space charge and actually increasing ob- 
served current. Since the width of this maximum was about one-tenth that of the 
minimum obtained for high values of V,, the former was adopted as the more sensitive 
criterion. It was also found that the band width was substantially independent of 
high frequency power input provided this was kept below a certain maximum. 

The proton resonance probe was placed in the field near the wave guide and the 
resonant frequency »,, determined by substantially the same method described in previous 
sections. The absolute frequency was not needed; instead the 657-th harmonic of vp. 
was compared with », and the difference was measured. A small error in the difference 
thus had a negligible effect on the result. The final value without diamagnetic cor- 
rection was given as velvi, = 657.475 + .008. 

Sources of error were carefully investigated.. The filament current was reversed 
and the positions of proton probe and electron beam interchanged. These results were 
averaged for each run; the deviation from the average under various interchanges 
was about 50 ppm. The maximum variation of any run from the overall average was 
less than 10 ppm and the statistical probable error of the mean was about 2 ppm. 
Various sources of systematic error were investigated and found to be a few ppm at 
most. The error stated is 12 ppm and seems to he closer to a limit of error than a 
probable error. It thus appears safe to reduce this figure to 6 ppm and to state the 
final result, without diamagnetic correction, as 


bol = 657.475 + .004 (6 ppm) . 


4. — Quantumelectrodynamic experiments. 


4°1.* Theoretical considerations. — The remaining experiments to be dis- 
cussed are concerned with the microwave spectra of hydrogen and deuterium. 
Due to the well-known techniques of frequency demultiplication and com- 
parison with a standard by means of beats, these measurements have a high 
inherent accuracy, which is of the order of 1 ppm. Consequently the necessary 
theoretical relations must be known to a corresponding accuracy in order to 
interpret the results properly. This requires the use of the recently developed 
theory of quantum electrodynamics to calculate a number of small but sig- 
nificant corrections. The experiments may then be regarded as fulfilling two 
different purposes: a) to provide a confirmation of the theory of quantum 
electrodynamics, and b) to furnish data for calculation of atomic constants. 
Fortunately there is enough data available to do both fairly satisfactorily. 


Quantum electrodynamics has been treated in a comprehensive work by HEITLER [146]; 
CorBEN [147] and Dyson [148] have briefly described some of the basic concepts in 
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excellent qualitative discussions. It is impossible to do more than briefly outline a 
few points in the present report. 

The theory began with Drrac’s [149] relativistic treatment of the electron. This equa- 
tion succeeded in giving the spin and magnetic moment as results rather than requiring 
them as additional assumptions; it had the apparent drawback of requiring an infinite 
number of negative energy states. Dirac postulated that these levels were normally 
filled with an infinite sea of negative electrons, which usually had no observable con- 
sequences. However if sufficient energy were supplied, e.g., by a photon, an electron 
could be transferred from the negative levels with energy E<— me? to a positive level 
with E>me?; the hole in the negative energy sea would then be recognized as a 
positively charged particle with the same mass as the electron. A few years later this 
prediction was verified with the discovery of the positron. 

Positrons are observed experimentally only in processes involving energies of the 
order of 1 MeV or more. However in processes of much lower energy there may be 
so-called « virtual »-electron-positron pairs. In a crude way these may be explained 
by saying that the uncertainty principle of quantum mechanics permits violations of 
the conservation of energy by an amount AF for a time At provided that AHAt ~ #. 
These virtual pairs are not observed directly, but may have significant indirect effects 
on interpretation of measured quantities. Thus Dirac’s theory of the electron is 
automatically a many particle theory. In some processes the possibility of virtual 
photons must also be taken into account; these processes result in what is termed 
vacuum polarization and the fluctuation of the vacuum field. 

Vacuum polarization may perhaps be clarified by another approach. If an electron 
is surrounded by other ordinary electrons, they will of course be repelled. However, 
electrons in negative energy states may be considered as having negative mass; hence 
the same force which would repel an ordinary electron attracts those in the infinite 
negative sea. This effectively leaves a deficiency of electrons in the negative sea (i.e. a 
positive charge) at infinity. Hence the effective charge of the electron is increased by 
the vacuum polarization. In this process the total charge remains constant, which 
means that an additional electron near the original charge is always accompanied by 
a hole (i.e., a positron) at infinity. Such processes are known as second order processes; , 
the corresponding terms in the various equations, which are of the order of «, are called 
second order corrections. ; 

Unfortunately, it turns out that the effective charge thus computed for an electron 
together with its induced charges is always infinite, even in vacuum. This has been 
a major stumbling block in the theory. However, the experimentally observed charge 
must always include the effects of the vacuum. Hence the value of e originally inserted 
must be unobservable; the infinite expression calculated from it must be identified 
with the experimental value of e. This process is called charge renormalization. Similar 
infinities arise in the calculation of the effective mass and self energy of the electron. 
However, when both charge and mass are renormalized as described above, all these 
infinities may be removed, as shown by Dyson [150]. 

Significant effects occur when the same calculations are carried out in the presence | 
of an electric or magnetic field. When these results are compared with the corresponding | 
ones for vacuum, finite differences are obtained; these differences should be observable. 
For example, it is found that the field of a nucleus departs significantly from a Coulomb | 
field due to vacuum polarization when the distance is of the order of the Compton 
wave-length (i.e. about 10-1° em). Near the nucleus the potential energy of an electron 
is somewhat higher (less negative) than one would expect from Coulomb’s law. The 
hydrogen-like wave functions [151] show that S electrons, have a high probability of 
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being found near the nucleus, while for other states the probability is almost zero. 
Fig. 12 shows the energy levels for n = 2 in H, computed by Lamp and others on the 
basis of quantum electrodynamics; they are compared with the Dirac levels and with 
the still simpler prediction of the Bohr theory. The separation is called the Lamb 
shift; this will be discussed further in a subsequent 
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spin along the field. It may be interpreted as 
a correction to the magnetic moment, which 
had previously been assumed to be exactly wo, Fig. 12. 

as given by the Dirac treatment. This is now 

known as the anomalous magnetic moment of the electron. It is of course 
really a second anomaly; the first lies in the fact an that electron with spin #/2 
has a magnetic moment 4, instead of po/2, as would be expected on classical 
grounds. KARPLUS and KROLL [152] have calculated the necessary corrections 
through fourth order terms (i.e., terms in «) and obtained: 


HYDROGEN FINE STRUCTURE FOR n=2 
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where «x denotes the fine structure constant (a = 27e?/he). Additional terms 
are of order x? and higher; thus, if the theory is exact, the error should be 
less than 1 part in 108. 


4°2.* Anomalous moment of the electron. — The anomalous moment of the 
electron has been experimentally measured to high accuracy; the theory is 
summarized by RAMSEY [153] and is basically similar to that of the Zeeman 
effect measurements [154]. If the proton had no magnetic moment, then the 
only magnetic moment in the ground state of H would be due to the electron 
spin. Thus there would be two degenerate levels in the absence of a field. 
A magnetic field would split these into two states with additional energies 
2u,Hm,, where m, = +} depending on the orientation of the spin with 
respect to the field. This hypothetical situation is shown by the two dotted 
lines in Fig. 13. 

In this case the Zeeman splitting of the two levels would correspond to a 
frequency f = 2u,H/h. The field H might be measured in terms of the nuclear 
resonance frequency of the proton[155], as previously described, giving 
H=h,/2u,. It would then follow that w,/u, =f/,- 

Actually the situation is complicated by the magnetic moment of the 
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proton, which is responsible for the hyperfine structure. This introduces two: 
additional energy terms. One is due to the interaction of the proton moment. 
with the magnetic field of the 
DR electron, while the other is due to 
‘its interaction with the external 
magnetic field. In -a field of a 
thousand gauss the former is of 
the same order of magnitude as. 
the interaction of the electron 
moment with the external field, 
hepo| 2000 while the latter is smaller by a 
MAGNETIC. FEL factor of about 10%. Thus the se- 
cond term, while not negligible, 
is of minor importance for fields. 
normally used in experimental 
work. 
At zero field the only signi- 
ficant factor is the interaction 
Fig. 13. of the proton moment with the 
electron field. When there is an 
orbital angular momentum of the electron, it is easy to see that when the 
total electronic angular momentum J is directed upward, the field at the 
nucleus is down; hence the position of highest energy occurs when I, the 
spin angular momentum of the nucleus, is 
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parallel to J. When J consists only of spin A 
angular momentum, i.e., in an S state, the Ts 
simple classical treatment breaks down and Bre 
the resulting expression becomes indetermi- CL 


nate (see RAMSEY [153]). However, a more 
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if F = I+ J is the total angular momentum, Fig. 14. 

the state Y = 1 (spins parallel) has a higher 

energy than Y = 0 (spins antiparallel), as indicated in Fig. 13. This dif- 
ference AW corresponds to a frequency »,, of about 1240 Hz. For low fields, 
in the order of 100 gauss or less, I and J may still be considered as coupled 
together to form a resultant F, which precesses about the field H. Hence a Zee- 
man splitting is obtained, the F = 0 level remains unchanged, while the P=1 
states split into three levels, corresponding to different values of M,. For 
high fields. i.e., a few thousand gauss, the coupling between I and J breaks 
down and they precess independently about H. These two cases are indicated 
ni Fig. 14. The relevant quantum numbers for high fields are m, and m,. 


Prepara = 
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In this case the primary splitting is due to m,; Fig. 13 shows that each pair 
of levels is split almost symmetrically and centered near the dotted line which 
would be obtained for zero proton moment; the level for which I and J are 
parallel is always slightly higher for fields under 10° gauss. The transition 
from the weak field to strong field case is analogous to that between the 
Zeeman and Paschen-Back regions in optical spectra. 

The region of experimental interest is near this transition. This involves 
the solution of a secular determinant for almost degenerate states (the two 
levels for which M, = 0) and has been carried out by BREIT and RABI [158]. 
Their result is 


W 
W(F,Mx.}=- — — 2u,HM, A +2M,0 + a, 


where the plus sign before the radical correspond to F=1 and the minus 
sign to F=0. The quantity x is defined by 


9 


SAW. 


(Lz =e Ly) ? 


in which both 4, and 4, are taken as positive quantities. It will be noted that 
the levels (1, 1) and (1, — 1) are linear functions of field throughout, just as 
in the low field approximation. 

If the field H is measured in terms of the proton nuclear resonance fre- 
quency [155] », (i.e., H = hy,/2u,) then u,H becomes hy,/2 and u,H becomes 
(hy,/2)(u,/u,). Since AW= hy,, the energy equation in terms of frequency 
units may be expressed as 


WF, My) = -5 —»,M,+2//1+2M, (2) (È a i) ag (2) (È fs i) 


H/ \Mp Yu} \Mp 


The observed frequencies can now be computed by taking differences between 
the appropriate levels; the transitions involved in the present experiments 
are those denoted by o, 2,, and x, (in analogy with optical Zeeman effect 
terminology) in Fig. 13. If the frequency of one transition is measured along 
with »v, (in the same field) and vy (without field), then clearly yw,/u, may be 
calculated. Since proton resonance is observed in oil the directly determined 
quantity is actually u alta,» A small relativistic correction [159], giving a slight 
difference between the magnetic moment of the free electron and that bound 
in the H atom converts the result to HH, 

KOENIG, PRODELL, and KuscH [160] have measured the anomalous magnetic 
moment using the z, transition (1, 0<>1, —1) indicated in Fig. 13. (The 
determination of v, will be discussed in the following section.) Their appa- 
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ratus is shown schematically in Fig. 15. The main electromagnet C produced 
a uniform magnetic field over a large area. The deflecting magnets A and B 
were of the two wire type (without 
iron core) and produced a strong in- 
homogeneous field over a small area. A 


a) TIBET beam of hydrogen atoms was obtained 
y (J DETECTOR from a Wood’s discharge tube at the 
SOURCE ne nay eae, left and collimated by a narrow slit. 
4 MAGNET M CAUCE) ‘he atoms were deflected in one direc- 
tion by the inhomogeneous field of A; 
rice this effect was normally just compen- 
ATOMIC BEAM APPARATUS sated by the field of B. Thus the beam 
a : was brought to a focus at the detector 
Fig. 15. 


on the right; the position of the focus 
was essentially independent of velocity. 
The detector consisted of a Pirani gauge; this measures pressure by the elec- 
trical resistance of a wire, which changes due to the effect of convection on 
its temperature. 

Transitions were obtained by applying a rf. field of 3655 MHz between the 
A and B magnets. This was furnished by a shielded, open-end rectangular 
box with small opening to admit the beam; it was fed by a coaxial line. When 
the field was adjusted to produce z, transitions at the applied radio frequency, 
this changed the direction of the magnetic moment of many atoms between 
the A and B magnets and a minimum was obtained at the detector. The 
proton resonance sample (not shown in the figure) was placed as near the 
beam as possible and allowance made for the gradient of the field. This cor- 
rection as well as small relativistic and diamagnetic corrections appear to be 
the chief sources of systematic error. Their final result was TA 653.2283+ 
+.000 6, where TA is stated for a spherical sample of mineral oil and the error 
is about twice the statistical probable error. 

BERINGER and HEALD[161] measured the same quantity by microwave 
absorption technique, using both the z, and 7, transitions (see Fig. 14). Atomic 
hydrogen was obtained from a Wood’s discharge tube and passed axially in 
a fused silica tube through a microwave resonant cavity. The cavity resonated 
in the TE,,, mode and was fed by a waveguide transverse to its axis. The 
waveguide continued on to a bolometer which measured the power transmitted 
through the cavity; this was, of course, greatly reduced when the radio frequency 
coincided with that of one of the hyperfine transitions. The radio frequency 
was fixed while the d.c. value of the magnetic field was varied. The field and 
the bolometer were both modulated at 30 Hz and the final signal obtained 
from a lock-in mixer; this technique greatly reduced background effects. One 
proton resonance unit was employed to maintain the d.c. level of the field 
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constant; a second was inserted in place of the resonant cavity and used to 
measure v,. Thus there was no problem about any spatial gradient of the 
field, but it was necessary to avoid or to average out time variations in the 
interval needed to replace the cavity by the proton probe. The final result 
(after corrections similar to those in the KoENIG, PRopELL and KuscH 
experiment) gave uu, = 658.2298 + .0003. However, the error represents 
only the statistical standard deviation and is obtained after rejecting 10 out of 
14 runs (giving high values) as statistically improbable. 

In view of the more conservative basis on which the KPK [160]error is stated, 
it seems reasonable to give equal weights to both measurements. The error 
of the mean will be taken as that given with the KPK result; this makes some 
allowance for any error in small correction factors common to both exper- 
iments. Thus the mean value is 


pu |, = 658.2293 + .0006 (1 ppm). 


This value may be combined with that of the GARDNER and PURCELL 
experiment [162] and the result compared with the theoretical value [163] pre- 
viously discussed. The experimental value thus obtained is u,/u, = 1.001147 
+ .000 006, while the theoretical value is u,/u, = 1.0011454. This remarkably 
close agreement is strong evidence in favor of the essential correctness of the 
quantum electrodynamie calculations. Consequently the theoretical ratio of 
ult, Will henceforth be treated as an accurately known theoretical relationship. 
(This admittedly assumes:a higher degree of accuracy than has been exper- 
imentally verified as yet.) Then, in view of the high accuracy of the anomalous 
moment experiment as compared with that of GARDNER and PURCELL, the 
latter will carry little weight in the final calculation of the constants, although 
it is of major importance in lending confidence to the recent theoretical de- 
velopments. 


4:3. Hyperfine structure of hydrogen. — The hyperfine splitting of the ground state 
of atomic hydrogen is the energy difference denoted by AW in Fig. 13. In frequency 
units this is equal to the quantity vg used in the last section. This is very nearly equal 
to the frequency of the o transition (1, 0<> 0, 0) in a low field; they differ only by qua- 
dratic (and higher) terms in the field strength [164]. The frequency vg has been expe- 
rimentally measured with high accuracy, but its theoretical interpretation is in some 
doubt. 

Narr and NELSON [165] determined vg by means of an atomic beam experiment; 
PRODELL and Kuscu [166] refined the method considerably, using essentially the same 
apparatus as in the KoENIG, PRODELL and Kuscu experiment described above and shown 
in Fig. 15. The chief difference was in the removal of the C magnet, only the residual 
laboratory field being used. An open loop of vertical copper plates shorted below the 
beam was employed and oriented to produce both x and o transitions. The 7, and zg 
lines were measured with lesser precision and used to make a small quadratic field 
correction to the o line; this amounted to about 1 ppm. The result of this experiment 
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then gave vg = (1420.4051 + .0002) s-1, where the probable error (14 ppm) is purely 
a statistical one. It should be noted that the full width of the resonance curve at half 
maximum was about 33 kHz; the ratio of this to yy is 23 ppm. PRopELL and KuscH 
point out that it would be difficult to prove the absence of asymmetry in the resonance 
curve sufficiently to justify confidence in the statistical error. However, they did 
succeed in eliminating certain experimental causes of asymmetry by taking data 
over.a wide range of rf. power input. 

Wirke and Dicks [167] determined vg by means of microwave absorption technique 
similar to that used by BERINGER and H®ALD in measuring Me| My They used a eylin- 
drical resonant cavity excited in the TE,,, mode with the sample in a glass container 
at the center. The atomic hydrogen was produced by a Wood’s discharge tube and 
mixed with molecular hydrogen to prevent it from diffusing to the walls too rapidly 
and then recombining. A series of runs were taken at three pressures and then extra- 
polated to zero pressure. The magnetic field strength was only .06 G; this was oriented 
parallel to the magnetie field of the rf. so as to produce only o transitions. The expe- 
riment gave vg = (1420.40580 + .000 05) -106 s-1, where the error is again only a sta- 
tistical one. 

The small but significant disagreement between these results may be due to the 
possible asymmetry discussed above. It therefore appears best to take the mean of 
these values and estimate the probable error conservatively on the basis of external 
consistency; the result is 


vg = (1420.4055 + .0003)-1068-1 (.2 ppm). 


The theoretical interpretation begins with the equation derived by FERMI [153, 157]. 
This may be qualitatively understood as follows. The hyperfine separation is due to 
the interaction of the nuclear moment with the field of the electron, which in an S state 
is due only to its spin. The field of a dipole falls off as 7°: hence the interaction energy 
should be proportional to p,u,/(7*),,- AS previously mentioned, the actual expression 
is indeterminate unless treated relativistically; however, (7°),, may be expected to 
be proportional to the cube of the Bohr radius a,. It follows that 


AW = My My _ Hy he \? (4x?me? * _eRo Hp ne 
to} ai uo) \4sme h? 2 PENI 


The detailed Fermi calculation gave a numerical factor of 16/3 rather than 4. 

This expression is in turn modified by several small corrections. Three of these 
are relatively straightforward: (1) Breit [168] obtained a correction factor of 3x°/2 by 
refining the mathematical approximations used in Fermi’s relativistic treatment. 
(2) Breir and MeyEROTT [169] showed an additional correction was necessary because 
of the motion of the nucleus. This is approximately equivalent to replacing m by 
the reduced mass in the expression for the Bohr radius; thus the correction is very 
nearly (—3m/m,). (3) When the anomalous moment of the electron is considered, the 
product mu, is equal (4,/1,)(Ms/Mo)"“o. thus involving the square of the theoretical ex- 
pression [163] for y,/u). BETHE and LONGMIRE [130] pointed out that two additional 
corrections are needed for a more exact treatment; i.e., allowance for distributed nature 
of the proton and electron magnetic moments rather than treating them as point 
dipoles. The electron moment extends over the order of a Compton wave-length; this 
effect can be calculated by a relativistic quantum electrodynamic treatment of the 
two particle problem; a small additional reduced mass correction is also obtained. How- 
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‘ever, the effect of the proton moment poses a more serious problem. While the 
electron moment is very nearly one Bohr magneton, the proton moment is almost 2.8 
nuclear magnetons. This anomalously high value is presumably due to its interaction 
with the meson field; unfortunately there is no generally accepted meson theory by 
which the necessary calculations may be made. 

KaARrPLUS and KLEIN [170] and KroLL and PoLLocxk [171] computed the second- 
‘order radiative corrections, which include the distributed nature of the electron dipole. 
Their calculations were in close agreement. The Kroll and Pollock result (with the 
correction of a typographical error in sign) was 


1602 Bolt 3a2\ /m,\3 ; = 125973018 
rg = cRo DE (pee so DE 3 
È 3 2/\m T 27 a? 


2.97302 5 
(Yo Tigri. re o2|—-—-In 2)| P,P, , 
Ms Qa a 2 


where m, = m(1 — m/m,) is the reduced mass. The first and second terms in paren- 
theses are Breit’s relativistic correction and the approximate reduced mass correction. 
The next two terms together give the proton moment in terms of the Bohr magneton. 
‘The bracket gives the electron moment in terms of the Bohr magneton plus additional 
radiative corrections. P, and P, are addition terms, then uncalculated, due to proton 
recoil and to the distributed nature of the proton moment. Using this equation, KROLL 
and PorLock obtained «1 = 137.03651 + .000 28, if P, and P, are set equal to unity. 

Further work has been carried out by ARrnowITT [172] and by NEWwcomB and SAL- 
PETER [173]. Their procedure essentially consists of dividing the proton magnetic 
moment into two parts — a Dirac moment which arises naturally from a relativistic 
treatment of the proton and a Pauli type moment which is added to bring the total 
up to the experimentally observed value. Their results contain a logarithmic diver- 
gence, which can be removed only by introducing a cut-off corresponding to the 
Compton wave-length of the proton or meson. ARNowiTt states that the result is 
insensitive to which one is used. NewcomB and SALPETER indicate that the dif- 
ference is significant; they conclude that this crude non-relativistic spreading of the 
nucleon bears little relation to reality. Consequently the theoretical situation still 
appears unsatisfactory. 

An indication cf this difficulty is found in the discrepancy between the theoretical 
and experimental ratios of the hyperfine splitting in hydrogen and deuterium, i.e. vg/vp. 
PRODELL and KuscH [166] point out that this difference is 170 ppm. Of course, this 
may be largely due to the more complicated structure of the deuteron. 

If the Kroll and Pollock result is corrected for the present values of the auxiliary 
constants (including a diamagnetic correction [155] of 29 ppm) without adding any 
new correction terms, the result is x-1 = 137.0369. It is difflcult to estimate the prob- 
able error properly. The principal experimental errors are in €, (,/“,), and the dia- 
magnetic correction, each being of the order of 1 ppm. However, the uncertain terms 
related to proton structure were originally estimated by Berne and LONGMIRE [130] 
to be of the order of 16 ppm; SALPETER [174] predicted that they might be more than 
20 ppm. (There is now some indication [173] that certain terms may be fortuitously 
small). Furthermore, terms in « have not been calculated for any corrections. Con- 
sequently it seems unsafe to put the probable error at less than 10 ppm in vg or 
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Hence. 


he 


= — = 137.0369 + .0007 
2rne? 


(5 ppm), 


where the error is given on a rather questionable basis. 


It should be remarked that the above uncertainties do not affect the use of vg in 


the anomalous moment experiments. These require only the experimental quantity vg, 
not its interpretation in terms of «. 


ti ENERGY (MHz) > 


4°4.* Fine Structure Separation. — The fine structure constant x was orig- 
inally determined by interferometric measurement of the fine structure sepa- 
ration of optical spectra [175]. BIRGE[176] pointed out that x could be obtained 
with higher accuracy from the L"L™ separation of X-ray lines; he calculated 
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Fig. 16. 


a value of x-1= 136.95 + .05 with 
the aid of a theoretical analysis by 
CHRISTY and KELLER [177]. BETHE. 
and LONGMIRE [130] later noted that 
this value should be raised to 137.11 
due to the anomalous moment of 
the electron. 

A great increase in accuracy can 
be attained by using microwave 
techniques. This makes it possible 
to measure directly the energy diffe- 
rence between two states of almost. 
equal energy rather than compute it 
from the small difference between 
two frequencies almost equal in mag- 
nitude. An extensive and painstak- 
ing investigation of the fine structure: 
of the first excited state of H and D 
has recently been completed at Co- 
lumbia University. The earlier work, 
including most of the theoretical 
calculations, was carried out by LAMB. 
and RETHERFORD [178], while the 
final experimental results needed in 
the present report were obtained 
by TRIEBWASSER, DAYHOFF, and 
LAMB [179, 180]. 


The quantity of primary interest is the separation AH between the 2Py 
and 2P, levels (see Fig. 16). A direct transition between these states is pro- 
hibited by the selection rule Al = +1. Hence AF can best be determined 
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by two measurements, the separation S between the 28, and 2P, levels and 
the separation AE — 5 between the 2P, and 28, levels. It should be em- 
phasized that the 28, level is only a convenient stepping stone to go from 
2P, to 2P,. Any other level to which transitions from these states are per- 
mitted would do as well in principle, although the precision might be much 
less. Thus any uncertainties in the theoretical interpretation of the Lamb 
shift S have no bearing on the experimental accuracy of measuring AE. 


These separations were actually computed from measurements in moderately high 
magnetic field. One reason was the large half-widths of the microwave lines. It was 
convenient to maintain a constant rf. power input of fixed frequency and to observe 
the resonance curve by changing the magnetic field. For this reason it is necessary 
to discuss the Zeeman splitting of the fine structure of hydrogen. 

Fig. 16 shows the energies as functions of field for the n = 2 state of H or D 
(which are almost identical). The notation designating different levels is the same as 
in the TRIEBWASSER, DAvHorr and Lamp work. Lamb’s [181] results for the various 
energy levels, in terms of frequency units, may be rewritten as: 


- For 2P3, m = + $# (a and d levels), 


AE Is MoH 
a aie dt : 
Yaa b= (1, a o) ( n 


For 2P3, m= +4 (b and e levels), 


For 253, m = +4 (« and f levels), 


S IsttoH 
a 


> 


where the plus or minus sign refers to the first or second subscript respectively. The 
exact form of the equations for y, , and», was given and used by Lamb. However, 
they are expressed here in Taylor series expansions through quadratic powers of H; 
this form is more instructive in showing the dependence of the levels on the different 
variables. The factor g,= 2/9 is given theoretically [163]; the effective g,=1—(m/m,), 
which includes a reduced mass correction, is known to very high accuracy. 

The a, d, x and f levels are linear functions of (“,H/h) and have an obvious physical 
interpretation; the slopes are approximately + 2 and + 1 respectively. The bd, c, e 
and f levels are almost linear functions of (“,H/h) with respective slopes of +. % and 
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+ 4; however, they also contain quadratic terms in H involving the ratio (u,H/AZ). 
For the fields employed in the experiment the quadratic term was always less than 
10% of the total field correction. Hence it was easily possible to assume a value of AF 
in the quadratic terms, calculate AH from the measured frequencies (e.g., ax and «f) 
and fields, and refine the result by successive approximations if necessary. 

It should pe noted that the hyperfine structure has been ignored in the above dis- 
cussion. Hence the lines in Fig. 16 really correspond to the dotted lines of Fig. 13. 
However, the hyperfine splitting is much smaller than in the ground state due to the 
increase in average electron radius; the S state is split by 88 MHz and the P state 
by only 29 MHz. The line widths of the observed transitions are such that the hyper- 
fine splitting is barely resolved, if at all. 


The experimental apparatus is shown schematically in Fig. 17. An oven, 
which dissociated about two-third of the H, molecules, was used as a source 
of atomic hydrogen (or deuterium.) A collimated beam of atoms passed into 
a magnetic field; here it was bombarded by a transverse beam of electrons 
accelerated by a grid (not shown) at roughly 13.5 V. As a result a small fraction 

(of the order of 1 in 10’) of the 


| . atoms was excited by collision 
SE TAPE to the 28, state, which is me- 
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of the P states. Since under 
operating conditions the ma- 
gnetic field was fairly close to 
the fe cross over point (574 G), 
most of the f component was 
removed through transitions 
to the e level; a small d.c. electric field in the rf. cavity furnished suffi- 
cient perturbation to make this practically complete. The beam then passed 
on to a tungsten plate, which acted as a metastable detector. In losing their 
excitation energy by collision, many of the metastable atoms were able to 
eject electrons from this surface. These were collected by a second electrode 
and the resulting current measured. 

If transitions from a 2S to 2P level were induced in the beam by an rf 
field, the atoms would then go quickly to the 1S ground state by spontaneous 
emission. This resulted in a lower current output of the metastable detector. 
At resonance peaks, about 30% of the metastable x component was « quenched » 
in this way. Two different rf cavities of appropriate design were used, de- 
pending on whether the desired transition required Am = 0 or Am = +1 
(o or x transition). Both were fed by a coaxial line from an oscillator with 
carefully controlled frequency and power level. Resonance curves were ob- 
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tained by scanning the magnetic field through the region of interest. In some 
cases the resulting curves were as wide as 100 MHz at half maximum. Since 
it was desired to obtain a result good to .1 MHz, all known causes of asymmetry 
were carefully investigated and corrections made wherever significant. 

The magnetie field was measured by means of a search coil and bridge 
circuit. For all runs used in final calculations, this was calibrated by means 
of a proton resonance probe. The calculation of the field corrections [180, 181] is 
rather complicated to follow; it appears to be equivalent to taking wo/h from 
BW [17] constants and measuring the field in terms of the THOMAS, DRISCOLL, 
and HIPPLE [136, 155] experiment. Since the field corrections total almost 20% 
of the value of AH, the error due to this correction and the correlation 
with other experiments should not be overlooked. This point was noted in 
the recent BEMT [2] evaluation; it was tentatively concluded that the ad- 
ditional error and the correlation were not too serious. 

A closer study shows that this argument can now be put on a much firmer 
basis. Calibrating the field in terms of the TDH [136] experiment meant 
essentially computing H from 


h il 
H = (5 = 2(—| My 
ZU] TDH Yp/TDH 


where v, was the measured quantity and y, was the TDH value. The field cor- 
rection in frequency units thus employed 


E 
h h Jaw Yo/TDH 


‘On the other hand, this quantity could be now computed as follows: 


(ol © He (Ta) E 5 (4) ; 
h h \2uy 7A 7779 ae 
The ratio (wo/#,) is known to about 1ppm from the experimental measure- 
ments [164] of (u,/u,) and the theoretical relation [163] giving (u,/“,). This value 
of (uo/4,) should be accurate to about 1 ppm, while the TDH experiment and 
the BEARDEN and WATTS [17] constants involve errors many times as great. 
However, the two methods will give the same results if it happens that 
the factor 27(Mo/h) py (1/y,)enn Used in the actual experiment equals the factor 
(uo/2u,) appearing in the second equation. These values are 328.7375 and 
328.7381 respectively, i.e., the difference is less than 2 ppm. Since the entire cor- 
rection is about 20% of AH, this difference amounts to about .4 ppm in AF, 
which is negligible. Thus any correlation between this experiment ando ther 
data may safely be ignored. Furthermore, even if the TDH value should be 
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substantially altered, it will not be necessary to revise the value of AH as 
long as (4/4) does not change appreciably. 

The Lamb shift S for H and D has been measured by TRIEBWASSER, 
DAYHOFF, and LAMB [179]. They use xe and of transitions at 1160 gauss and 
700 gauss respectively (2195 and 2395 MHz), as shown in Fig. 17. The mean 
of these two values for deuterium is 1059.00 MHz with a statistical probable 
error of .03 MHz based on internal consistency of the two measurements. Actually 
each resvlt differed from the mean by about .05 MHz, which gives a very reason- 
able statistical picture. In the case of H, where asymmetry corrections for 
the two transitions differed by 2 MHz, the close agreement between the re- 
sulting values of S gave strong evidence for the validity of the corrections. 
(The corrections for D were smaller and almost equal.) The field measurements. 
have already been discussed; there is no other obvious source of systematic 
error. The experimenters multiplied the average deviation of the mean by 
three to obtain a limit of error of .1 MHz. They gave no specific estimate of 
systematic error here, although in the experiment to be discussed below, they 
allowed .05 MHz, intended as a limit of error. On the basis of all this it seems. 
reasonable to adopt a probable error of .05 MHz for the S measurement. 

It may be noted in passing that the experimental value of S differed from 
the theoretical result given by SALPETER [182] by about .5 MHz. While this dis- 
crepancy poses an important problem, the agreement may still be regarded 
aS a major achievement of quantum electrodynamics. It should be empha- 
sized that the entire Lamb shift, not merely some small corrections to it, is. 
due to quantum electrodynamic effects. As already noted, this discrepancy 
is of no real concern in the fine structure measurement unless it is taken to 
mean that the experimental value of S is in error. 

DAHYOFF, TRIEBWASSER, and LAMB [180] measured AH — 5 using the ax 
and ca transitions at 630 and 1190 gauss respectively (10795 and 7195 MHz). 
They rejected the «e measurements from the final result; these were based 
on a method of field calibration using an «f transition, which was proved to 
be less dependable than the proton resonance technique. (The difference 
between the ax and ca results was .2 MHz.) They obtained a limit of error by 
adding 3 times the probable error (.017 MHz) to an estimated possible system- 
atic error of .05 MHz. <A probable error of .05 MHz will be used here. Adding 
the results of the two experiments and combining probable errors in the usual 
way gives 


E 
Van = = = (10971.58 + .07) MHz (6 ppm). 


It is now necessary to use this value in a theoretical relation to calculate «.. 
The order of magnitude may be quickly obtained by recalling that this is an 
interaction between the spin moment of the electron and the magnetic field 
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created by the motion of the proton as seen from the electron; the separation 
will be of the order of the Bohr radius a,. Hence the energy should be of the 
order of u3/43; this is roughly the same as the estimate for the hyperfine splitting 
with the omission of the factor (u,/u,). The necessary calculations are sum- 
marized by LAMB [183] and the result stated in the fine structure experi- 
ment [180]. It may be rewritten as 


: Fa | 5 a 5.946 I 
alt= 1+-a% + —-— OB ae 
8 I 


16v4z ora 


where »,, is the frequency corresponding to AH. The first two terms come 
from the Dirac treatment [149], while the last two are based on the second and 
fourth order anomalous moment corrections. The equation may, of course, 
be solved by successive approximations. Correcting Dayhoff’s figure for the 
values of ¢ and R, adopted in the present report gives. 


ot = = 137.0371 + .0004 (3 ppm). 


2re* 
‘The above equation is exact,.so far as is known, except for higher order terms 
in x. Hence, unless these are unexpectedly large, there should be no theo- 
retical uncertainty in the value. 


C) Analysis 
1. — General principles. 


The experiments just discussed represent a rather formidable mass of data. 
It is fairly obvious that experiments with errors much greater then 100 ppm, 
will carry negligible weight in the final adjustment. However, it seems de- 
sirable to adopt some criterion for accepting or rejecting borderline data. 

The relative error by itself is not an infallible standard. For example, the 
RyMER and WrIGHT determination of h/(em)? has been assigned a probable 
error of 70 ppm. It could equally well have been listed as a determination 
of h?/em. In this case, the error would be doubled, but the significance of 
the data would be unchanged. 

Consequently the following procedure will be adopted: The constants will 
first be calculated by a direct solution (i.e., one which is not overdetermined) 
of the most precise experimental data. (A unique choice of such data will 
not always be possible, of course.) This experimental data will be referred to 
as the «trial data». From the trial data it is also possible to calculate any 
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desired combination of constants with the corresponding probable error. The 
error of any other experiment may thus be compared with that obtained from 
the trial data. This will be illustrated by an example later on. 

At this point it should be noted that the probable error of an experiment 
is never known exactly. For example, it was shown by GAuss [184] that on 
a purely statistical basis the relative probable error of a probable error is 
.48/Vn, where n is the number of observations. Thus with 25 observations 
the probable error itself may have an uncertainty of about 10%. In the 
experiments considered here, the statistical probable errors are subject to this 
order of uncertainty; the systematic errors are generally much more difficult, 
to evaluate. Since data will be weighted as the square of the probable error 
in the final adjustment, the weighting is subject to an uncertainty of at 
least 20%. 

This discussion shows the fallacy of including data of low weight. If one 
experiment has three times the probable error of another, it will carry only 
about 11% as much weight. However, the weight of the second experiment is 
uncertain by more than this amount. Hence, as a rule of thumb, it seems 
pointless to include any experiment when its error is more than three times 
that of any other combination of experiments yielding the same result. There 
might be an exception if a great quantity of additional data where available, 
with each item having an error barely too high to meet the above criterion. 

By the same token when some items used in the adjustment have errors 
much smaller than the others, it is pointless to treat such quantities as variables. 
for adjustment. These quantities may be then termed auxiliary constants; 
the value and error of each one is fixed by a single type of experiment. This 
procedure assumes that such constants do not overdetermine any variable or 
appear to too high a power in any calculation. 


2. — Selection of input data. 


Table VI summarizes all of the experimental results with errors of 100 ppm 
or less (as well as the theoretical relation for w,/u,); the left hand column 
indicates the experiment involved and the section of Part B in which it is 
described. The first group comprising five relations all involve errors of 1 ppm 
or less; they do not overdetermine any of the unknowns. Clearly they may 
be treated as auxiliary constants; in fact they have already been so used in 
many of the other equations. 

A direct solution is possible with the first two groups of equations. With 
the help of the auxiliary constants, Eqs. (6), (7), and (8) can be rewritten in 
terms of N, e and h. Obviously e may be eliminated wherever it appears by 
Eq. (4); m may be similarly treated by Eq. (2). Through the use of Eqs. (1) 


Yepreeo ne er e-4 


Rydberg constant (3°12) 


Theoretical anomalous 


of electron (4°2) 


Velocity of light (2°1) 


Mass of proton (2°8) 


of D (44) 


Cyclotron resonance of 


proton (3°14) 


Gyromagnetic ratio of 


proton (3°13) 


Electrochemical 


faraday (2°3) 


Avogadro’s number 


Gas constant (2°9) 


Cyclotron resonance of 


electron (3°15) 


H yperfine structure 


of H (4:3) 
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TaBLe VI. — Summary of experimental data. 
Note. - All values are in c.g.s. units. e is given in e.s.u. and magnetic moments in erg/Oe. 

Experiment Quantity Beatlt Error Eq. 
(Section Ref. Part B) measured N (ppm) Ne. 
Mello - 1.001 1454 0 (1) 

moment of electron (4°1) 
Ro = 2n?met/ch? 109 737.31 DI (2) 
Exp. anomalous moment Ms| My 658.229 3 1 (3) 
e 2.997928-1010 1 (4) 
MU; 1.007 5957 1 (5) 
Fine structure separation |. «1 = he/2re? 137.037 1 3 (6) 
RITA 2.792 685 7 (7) 
ITAL. 2.675 23-104 15 (8) 
iH = Neje = 9651.4 30 (9) 
N'= (4,/4,)*N 6.0614-1023 33 (10) 

(X-ray Cryst. Dens.) (2°6) 
Ruled X-ray gratings (2°7) AglAs 1.002 02 20 (11) 
R - 8.38170-107 12 (12) 
Molt, 657.475 6 (13) 
at — he/2ne* - 137.0369 5 (14) 
X-ray high freq. limit (3°5) (h/e)(Ag/ Ag) 1.37653: 10-17 65 (15) 
.999 54-10-8 70 (16) 


Electron diffraction (3°6) | [h/(em)*](A,/A,) 
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and (3), TE is expressible in terms of 4, which in turn is defined as eh/4ame. 
Since, by definition, the nuclear magneton pw, =,(m/m,) and M, = Nm,, 
Eq. (7) may be treated as follows: 


A Bi (2 bo = 
M, Ms My Mn Nm, Ho(M/Mmp) Nm 


The first three quantities on the left are given by Egs. (5), (1) and (3), 
while on the right m may be eliminated by Eq. (2). 
Thus Eqs. (6), (7) and (8) become 


(17) hje® = 2.8720818:10-8= B, (3 ppm), 
(18) e4/Nh? = 3.0371164-10!7 = B, (7 ppm), 
(19) e°/h*+*=;817883001033=—B, (15 ppm). 


None of the errors is significantly changed by the auxiliary constants; in the 
case of Eq. (17) reference to the original calculation [185] of a! will show that 
ec appears only as c-? in the final result. 

The solution of these equations is 


e = 4.802876-10-10 
(20) h = 6.625208-10-27 
N=6.024837-1023. 

The remaining items in Table VI will now be considered: it will first be 
shown that the next group, Eqs. (9), (10) and (11), should be retained. In 
terms of the trial data with the B's defined by Eq. (17), (18) and (19), it is 
seen that Y = Ne/c = B:/(B,c), with a corresponding error of (15*+7?)* =17 ppm 
Since the error in the electrochemical value of F given by Eq. (9) is less than 


twice as great, it will carry appreciable weight and should be included. With 
elimination of c it becomes 


(21) Ne = 2.893420-1014 (30 ppm). 


Eqs. (10) and (11) may be combined to give N with an error of 68 ppm. The 


trial data give 
Noa BoB Bi 


P.E. =[(3x3) + 72+ (2x15)?] = 32 ppm. 


Since this error is only about half that obtained from Eqs. (10) and (11), they 
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too will be included. The resulting equation is 
(22) N = 6.024817-1023 (68 ppm). 


The remainder of the data, Eqs. (13) to (16), will not be used in the adjust- 
ment. As previously discussed, the GARDNER and PURCELL result Eq. (13) 
can also be obtained from Eqs. (1) and (3); therefore, once the theoretical 
relation for u,/u, has been accepted unreservedly, Eq. (13) carries little weight. 

The hyperfine structure measurement Eq. (14), determines the same quan- 
tity as the fine structure result, Eq. (6). The hfs. result carries a somewhat 
lower weight and differs from Eq. (6) by less than half the probable error of 
the latter. Therefore the practical effect of including or rejecting Eq. (13) 
will be small; in view of the theoretical uncertainty involved in the hfs. inter- 
pretation it seems safer to omit this item entirely. 

Combining Eq. (11) with Eq. (15) gives h/e with an error of 68 ppm; 
combining it with Eq. (16) gives h/(em)* with an error of 73 ppm, The cor- 
responding errors calculated from the direct solution data are 16 ppm and 
10 ppm respectively. Hence both Eqs. (15) and (16) carry too little weight 
to be of value. 

It should noted that the ratio /,/A, need not be treated as a variable at 
this stage of the adjustment since it appears in only one equation which is 
to be included, (Eq. (22)). Thus the data can not be correlated through this 
factor. After the primary constants, e, h, N, and m, have been determined, 
the best values of A,/A, will be considered separately. 


3. — Linearized equations. 


Five equations have now been selected, viz., Eqs. (17), (18), (19), (21) and (22). 
Mathematically these equations represent five surfaces in the three dimensional space 
of the variables e, h and N. If the experiments were all ideally precise, these 
surfaces would intersect in a common point which would then give the true values 
of the variables. Of course, in practice they do not quite fulfill this ideal; hence a 
least squares treatment becomes necessary to find the best compromise solution. 
However, since the surfaces all come very close to meeting in a common point, 
they may certainly be approximated as planes in the region of interest; in other words, 
the equations may be linerized. 

Consider for simplicity an equation involving two variables, e and N; let f(e, N) = 
= e*N? — B, where B is the experimentally observed value. A pair of «origin values », 
e and N,, which are known to be close to the final solution, may now be selected; in 
this analysis, the values obtained from the trial data will be used. The function f(e, N) 
may now be expanded about this point 0 in a Taylor series, only linear terms being 
retained. Thus 


È of of <a 
fle, N) = fle, No) + (=) (€ --&) + (as) (N-MN) Hosp =Bs 
0 


de ON} 


22 — Supplemento al Nuovo Cimento. 
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Since f(e, N) = e*N®, the linear approximation gives 


B— Wt = (aes INo)(e — e0) + (be NS N — Mo). 
Multiplying by 16%/e} NG yields 


BEN; A 
(23) AX, a= bay = TEN ‘ 10 > 
where 
€—Cy N—WN, 
(24) E = “108, Dar -108. 
lo No 


The w’s represent differences in ppm between the values of the variables and the 
origin values; the term on the right gives the difference in ppm between the experi- 
mentally observed quantity and the value which would have been expected on the 
basis of the origin values. It will be noted that the coefficients of the x°s are just the 
exponents of the original variables. The extension to three variables is straightforward 
and yields equations of the same form. 


The linearized forms of the five equations to be used are given below, with 
the original number of the equation shown on the left. Since the direct so- 
lution values of Eq. (20) are used as origin values, the constant terms of the 
first three equations all vanish. Experiments are weighted inversely as the 
square of the probable errors; the weighting factor is 900/(P.E.)?, where the 
probable error is expressed in ppm. Some of the weights have been rounded 
off slightly, but this is not significant in view of the uncertainty of the errors. 


Original Eq. Linearized Equations Weight Eq. Number 
! 
(17) [Son Berea =~ 0 100 (25) 
(18) | da; = dt, ay =i. 0 ' 18 (26) 
(19) Sari nariri = 0 4 (27) 
(21) | xe + xy = — 81.0 1 (28) 
(22) | ty = — 3.3 i (29) 


4. — Least squares solution. 


The most generally accepted means of obtaining the best compromise so- 
lution for a set of overdetermined equations such as given above is the method 
of least squares. This is a well-known technique, which is discussed in numerous. 
mathematical texts. (See, for example, CRAMER [186] or WHITTAKER and Ro- 
BINSON [187]). The application of this technique to adjustments of atomic con- 
stants has been well treated by DUMOND and COHEN [1, 29, 188]; COHEN has 
extended the method to the case when the experimental errors are corre- 
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lated. Some of the principal results needed here are summarized in Appendix 
B; for further details the above references should be consulted. 

It should first be noted that there are no serious correlations between the 
above equations. Several of the equations are correlated through e, R,, and 
u,lu,; but each of these factors involves errors of about 1 ppm or less. In 
the case of the equation with highest weight, Eq. (25), c actually appears 
only as c? as previously noted, while u.Ju, does not appear at all. Eqs. (27) 
and (28) are correlated by the use of the same electrical standards to deter- 
mine the ampere. The error in the ampere [137] has been estimated at 6 ppm; 
this is about one-third the uncertainty of Eq. (27) and only one-fifth that of 
Eq. (28). Thus, in view of the uncertainties of the errors, it seems an um- 
warranted refinement to allow for the effects of these correlations. 


The procedure for obtaining the normal equations is derived in Appendix B. For 
each variable an equation is found corresponding to minimizing the sum of the squares 
of the residuals with respect to that variable. This is obtained by taking the coefficient 
of the variable in each linear equation times the weight of the equation and multiplying 
through the linear equation by this factor; all of the resulting equations are then added. 
For example the equation corresponding to minimizing the sum with respect to xy 
is obtained as follows. 


O—1-18(4r, = 3x,--xy—-0) + 0 + 1-1:(x. + xy + 81.0) + 1-.2(xy + 3.3) = 0, 
—Tlax, + 54x, + 19.2xy + 81.66 =0. 
Equations corresponding to minimizing with respect to x, and x, are similarly derived. 


The result is 
normal equations 


71892, — 476%, —7lzy = —81.0, 
(30) —476x, +2980, +54¢, = 0, 
— Tlx, + 54a, + 19.22, = — 81.66 


The solution of the normal equations is 


pes 1341, 
(31) Vaia 261 
Uy = ae OT . 


In terms of the original variables the solution becomes 


e = 4.802938-10-10, 


I 


(32) h 6.625381-10-27, 


N = 6.024 658-103 , 
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All of the figures given above are not, of course, significant in these particular variables. 
However, they will be retained in calculating derived quantities; the errors will then 
be calculated and the results appropriately rounded off. 


5. — Errors. 


Before the errors are computed the y? test [189] will be applied to determine the 
consistency of the data. Clearly some such criterion is necessary to decide whether 
the errors have been assessed on a reasonable basis, i.e., whether the discrepancies 
between various experiments are of the same order as the experimental errors. The 
quantity x? is defined by 


d;\} 
330) 


where o; is the error associated with the i-th linear equation and d; is the difference 
between the right and left hand sides of the i-th equation when the least squares 
solution values are substituted for the w’s. It can be proved that the expectation value 
of x? for N equations and n variables is N —n. Roughly speaking this means that if 
there were one variable and many experiments (equations) each term in the sum 
would have an average value of unity and that y? would be approximately equal to N. 
Since n is of the same order as N, the discrepancy between the observed values and 
the solution values of the «’s is less than that between the observed values and the 
true values; in fact, if N= », all d; must vanish and y® = 0. Hence it seems reason- 
able that the expected value of y? is equal to the degree of overdetermination of the 
equations. 

When the solution values of Eq. (31) are substituted back in Eq. (30), the result is 


I 

bo 
a 
bo 
ot 


(33) ie 


The expected value is5 — 3 = 2. The value actually obtained seems quite reasonable, 
since large statistical fluctuations are to be expected when the degree of overdeter- 
mination is so small. 

The ratio of error by external consistency (based on agreement between experi- 
ments) to that by internal consistency (based on estimated errors of various experi- 
ments) is obtained by dividing z? by its expected value (N — »), and taking the square 
root, since y? involves the square of the differences. Thus 


(P.E.)ext/(P-E-)int = V(2.625/2) = 1.146. 


Hence the errors are about 15% larger by external consistency. Since this basis is 
thus more conservative, all errors in the final results will be computed by external con- 
sistency. 

As shown in Appendix B, the error matrix (based on internal consistency) is simply 
Kg", where K is the weighting constant (the factor 900 used in computing the weights 
for Eqs. (25) to (29)) and g is the inverse of the matrix of coefficients in Eq. (30). 
These values must be multiplied by 1.1462 to express the errors on the basis of external 


Lai. 
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consistency. The matrix then becomes 


e h N 
298.2 563.9 — 483.2 
eA (Kee) = 563.9 1074.4 — 936.4 
— 483.2 — 936.4 908.3 


These terms are in (ppm)?. The diagonal terms represent ¢e> i.e. the squares of the 
errors of e, h and N. Thus, for example, the probable error of h is (1074.4)? — 33 ppm. 
The off-diagonal terms are needed in calculating derived quantities, as explained below. 

The matrix of correlation coefficients is defined by 7,;; = Kee] Ke;)Key). This is 
easily computed to be 


e h N 
1.000 .996 — .928 
eS) (ra) = 996 1.000 — .948 
— .928 — .948 1.000 


It will be noted that there is a high degree of correlation among these variables. For 
example, if it were found by an entirely new and more accurate experiment that the 
value of e should be raised, then it is almost certain that the value of h should also 
be raised, while the value of N should be lowered. 


6. — Errors in derived quantities. 


The high correlations just obtained show clearly that the errors in e, h and N can 
not be treated as independent in determining the errors of derived quantities. How- 
ever, the correct procedure is quite simple when the error matrix has been calculated. 
As shown in a previous example [190], a function f(e, N) = e*N® may be linearized in 
the region of interest and expressed as 


(f —fo)/fo = a2, + bay. 


If the relative errors in e and N are e, and ey respectively, then it follows that e,, the 
the resulting relative error in f, is given by 


& = de, + bey. 
Hence the square of the probable error in f is given by 
<> = <(ae, + bey)®> = a2e$> + 2ab<esen) + b*ey) - 


The generalization to three variables is straightforward. 
As an illustration, the error in the fine structure constant x = 2me?/he will be com- 
puted. The factor 27 is known exactly; since c is an auxiliary constant, its error may 


a Pa Ay > sia? (ey er pe Oe, eS eee eS A 
È t K 
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be neglected. Thus the probable error of e?/h must be computed. This is given by 


Ke) = 482) + 2+2+(—1)<e,e,> + <> 2 
<e2) = 4:298.2 — 4-563.9 + 1074.4 = 11.6. 


The probable error in « is then 4/11.6 ~ 3.4 ppm, which is much less than the probable 
error in either h or e. This result was to be expected since « is determined principally 
by the fine structure experiment, Eq. (25), which gave « with an error of 3 ppm. 

It may be noted that the correlati:;n coefficients in Eq. (35) can never be greater 
than unity. Consequently the error of a product of two factors can never be greater 
than the sum of the errors of the individual factors nor less than the difference. Hence 
correlation coefficients may have an important role only when the errors of the two 
factors are comparable in magnitude. 


7. Calculation of 4,/4.. 


It still remains to calculate /,/A,, one of the primary objects of this inves- 
tigation. There are now available two ways of determining this quantity. 
The direct ruled grating measurement, Eq. (11), gave 


(36) 7,1: = 1.002 02 (20 ppm). 


On the other hand, it may also be determined from Eq. (10) together with the 
least squares value of N from Eq. (32). This yields 


(37) AglAs = 1.002 029 (11 ppm for weighting) 


where the error quoted represents only the error of Eq. (10) after taking the 
cube root. If these are weighted inversely as the square of the errors quoted, 
the mean is A,/A, = 1.002027. It may be shown that this result is exactly 
the same as would have been obtained by including 4,/4; as one of the variables 
in the adjustment. Since the algebra is somewhat lengthy, the proof of this 
statement will be omitted. 

The error can not be calculated rigorously by a simple combination of the 
errors in Eqs. (36) and (37). However, it is certainly less than the 20 ppm 
of Eq. (36), on the basis of internal consistency. By assuming for the moment 
that Eqs. (36) and (37) are completely independent and including the least 
squares error in N in Eq. (37), it follows that the error in the mean value of 
As/s is about 12 ppm. Thus the error would appear to lie betwen 10 and 
20 ppm. Since most errors are to be rounded off to one significant figure, a 
further study seems superfluous. The result will be stated as 


(38) 2] = 1.00203 + .00002 . 


boe 
> 
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In the judgment of the writers the above value is the best choice on the basis 
of the data available. However, it differs from the accepted value [3] now 
in general use by about half a probable error. It does not therefore appear 
possible to make a strong case for ahy change in the accepted value. 


8. — Transformation of variables. 


For the purpose of plotting a consistency chart, it is convenient to make a 
transformation from the quantities e, h and N to new variables «’, F’ and 6. The 
transformation is indicated below, both in terms of the original variables and in the 
linearized forms 


a — Fl hy Lar = 2%, — > 
(39) F'= Ne Cp =U, + ay, 
Scie Le = 5x, — 3%. 


Eqs. (25) to (29) now transform to 


(40) — ty SIRO (3 ppm), 
(41) —tp+xg= 0 (7 ppm), 
(42) vg = (0) (15 ppm), 
(43) Cp = — 81.0 (30 ppm), 
(44) — ty + 2p + %g= — 3.3 (68 ppm) . 


Combining Eqs. (40) and (44) yields 
(45) tp +xg= — 3.3 (68 ppm) . 


Actually this transformation might have been employed earlier. «’ could have 
been treated as an auxiliary constant since it is essentially determined by one equation, 
and the adjustment carried out with two variables. The final result would not have 
been changed by more than small fractions of the probable errors. 


D) Conelusions 


1. — Recommended values (with tables). 


Table VII gives the values obtained for the principal constants; Table VIII 
lists the conversion factors for various energy units in common use. These 
values are recommended by the writers for general adoption. 
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It should be noted that nuclear masses have not been studied as critically 
as other data, although they are believed to be of ample accuracy for use as 
auxiliary constants in deriving most other results. Also the value of pV at 
standard conditions has been taken from chemical sources [191] and not studied 
in detail; this affects the value of R, k, V,, and constants derived from them. 


TABLE VII. — Values of atomic constants as of August 1955. 


Note. — These constants are arranged alphabetically according to the usual symbols 
employed for them; Greek letters are included after the English equivalents. Atomic 
masses are all denoted by M with the proper subscripts. All atomic masses are on 


the physical scale. 


Ay First Bohr radius 

x Fine structure constant 

at 

c Velocity of light 

cr = 8ahe First radiation constant 

6, = helk Second radiation constant 

e Electronic charge 

elec 

ejm Specific electronic charge 

e/me 

F Faraday constant 

Vp Gyromagnetic ratio of proton 

Vo Apparent value of y, in oil 

h Planck’s constant 

h 

h/e 

hm 

k Boltzmanu’s constant 

ie Compton wave-length of 
electron 

dig Compton wave-length of 
proton 

Ag! As Conversion factor from Sieg- 
bahn X-ray scale to ab- 
solute units 

Jed Wien’s displacement law 


constant 


(5.29173 + .00002)-10-? cm 
(7.29729 + .00003) -10-% 

137.037 1 + .0005 

(2.997928 + .000004)-101° cm s* 
(4.9920 + .000 2) -10-15 erg cm 
(1.43879 + .00003) em °K 
(4.80294 + .00008)-10-1° esu 
(1.60209 + .00003)-10-2° emu 
(5.27299 + .00008)-1017 esu g-t 
(1.75888 + .00003)-107 g-1 s cm? 
(9652.0 + .2) emu eq! 

(2.67527 + .00006)-104 Oe-1 s-t 
(2.67519 + .00098)+104 Oe s-1 
(6.6254 + .0002)-10-2" erg s 
(1.05446 + .00004 
(1.37944 + .00002 


)- 10-7 ergs 

) 
(7.27378 + .00005) 

) 

) 


em? s-1 
(1.38049 + .00005)-10-16 erg °K-1 


-10-17 erg s esu-? 
*.10-=1%cma 


(2.426 27 + .00002 
(1.321 40 + .00002)-10-13 em 


1.00203 + .00002 


(.289 780 + .000005) em °K 
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(TaBLe VII continued.) 


Ao Wave-length associated with 
lev 

My Atomic mass of deuterium 

Ma Atomic mass of deuteron 


M.= Nm Atomic mass of electron 


My Atomic mass of hydrogen 

Ma Atomic mass of neutron 

M, Atomic mass of proton 

M,/M, Ratio of proton to electron 
mass 

m Electron mass 

Mo, Neutron mass 

My Proton mass 

bu Reduced mass of electron in 
H atom 

Un Nuclear magneton 

Ho Bohr magneton 

Up Magnetic moment of proton 

Ms Magnetic moment of electron 

Hs/ Ho 

Ms/ My 

N Avogadro’s number 


n= N/V, Loschmidt’s number 


Po Standard atmosphere 

R Ideal gas constant 

Rp Rydberg constant for D 

Ru Rydberg constant for H 

Bio Rydberg constant for infinite 
mass 

r Ratio of physical to chemical 
scale of atomic weights 

o Stefan-Boltzmann constant 

Pa Molar volume 


*) Conventionally adopted values. 


(1.23978 + .00002)-10-4 em 


2.014740 + .000002) amu 
2.014191 + .000002) amu 
5.487 60 + .00004)-10-4 amu 
1.008 1445 + .0000010) amu 
(1.008 984 + .000002) amu 
(1.0075957 + .0000010) amu 
1836.13 + .02 


( 
( 
( 
( 


(9.1086 + .0003)-10-28 g 
(1.67476 + .00005)-10-24 g 
(1.67245 + .00005)-10-24 g 
(9.1036 + .0003)-10-28 g 


(.50505 + .00002)-10-7* erg 0e-1 
(.92733, + .00002)-10-2° erg Oe? 
(1.41048 + .000 03)-10-?% erg Oe-* 
(.92840 + .00002)-10-2° erg Oe-1 
1.001 1454 (theoretical) 

658.211 + .001 

(6.0247 + .000 2)-10?8 mole 
(2.68709 + .00009)-1019 em-3 
1.013 250-106 dyn em-? (*) 


(8.3170 + .0001)-107 erg °K-1 mole 


109 707.42 + .01 cm“! 
109 677.57, + .01 cm-t 
109 737.31 --’.01 cm 


1.000275 (*) 
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(5.6696 + .0004)-10-5 erg em-2 s-1 °K-* 


(2.24208 + .00003)-104 em? mole! 
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TaBLe VIII. — Energy conversion factors. 
erg V sai emo SK: 
1 erg 6.2419 -10!! |1.50935 -10265.0346 -10!5 |7.2438 -1015 
1.602 09-10-12 1V 2.41810 -10!48.0659 -103 |1.16052-104 
6.6254 -10-274.13547-10-15 1s-1 3.335 637 -10-11/4.799 29-10-11 
pee oe 1.239 78-10-4 |2.997 928-101° 1 em 1.43879 
1.38049-10-158.6168 -10-5 2.083 64 -10!9 .69503 ek 
_ |8.987 57+ 10°° 5.6099 -1082 1.35654 -10474.5249 1036 /6.5104 -1086 lg 6.0247 
|1.49180-10-° 9.3116 -108 |2.25164 -10?3/7.51066 -10!2 |1.08063-10! |1.65985 -10-24 lamu 


Note. - The above table is to be read horizontally. For example, 1 erg is equal to 6.7033°10? amu, i.e., 
the conversion factor is 6.7033-10* amu/erg. 
The probable error indicates an uncertainty in the last digit of each factor. 


The value of A,/A, quoted here is slightly higher than the presently accepted 
value of 1.00202. However, since the difference is only about half the prob- 
able error, there is no compelling reason for discarding the present widely 
accepted value. 

As elsewhere in this report, the errors given are probable errors; it is there- > 
fore to be expected that the true values lie outside the stated limits in some 
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2. — Consistency 
of data. 


In many stu- 
dies of the con- 
stants an isometric 
consistency chart 
introduced by Du- 
Monp [192] has 
been used to pre- 
sent the adjust- 
ment graphically. 
However, in the 
present instance a 
still simpler picture 
may be obtained by 
using the transfor- 
mation of variables 


previously given 
[193]. This trans- 
formation shows 
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that the fine structure equation is nearly independent of the other equations and that 
2' may be used, without serious error, as an auxiliary constant in the one other equation 
where it appears. Thus the essential features of the adjustment may be shown in a 
two-dimensional plot using the variables 7'= Ne and f = e5/h3. This is given in Fig. 18. 

Since F’, is proportional to the faraday constant F = Ne/c, the effect of changing 
the electrochemical value of F may easily be seen by moving the vertical line repre- 
senting this experiment to the right or left by the desired number of ppm. By a for- 
tunate coincidence the other three equations all contain “vg to the first power. Thus 
the effect of changing any one of these experimental values by a given number of 
ppm may be seen by moving the line parallel to itself and changing its 6 intercept by 
the same number of ppm. 
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Fig. 19. 


The adjusted least squares values and the ellipsoid of error are also indicated 
in Fig. 18. It is clear that the adjusted values are in excellent agreement with all 
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experiments except the electrochemical faraday and that the discrepancy is not too 
implausible even in this case. 

The other experiments, which were not used in the adjustment, will now be con- 
sidered briefly. The various graphs in Fig. 19 compare the adjusted least squares values. 
with the less accurate data. The least squares values are denoted by « LS»; in most: 
cases, where the corresponding errors are too small to be seen on the scale used, the 
LS values are simply indicated by crosses. The origin is taken as the LS result, which 
is the numerical value given; the deviation of other experimental values and the 
spread of the probable errors are then shown in terms of ppm. The oil drop deter- 
mination of the electronic charge is omitted for reasons previously stated [194]. 
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Fig. 20. 


Almost half of the probable errors overlap the least squares values; in no case 
is the discrepancy greater than three probable errors. This may be considered as fairly 
satisfactory agreement. One of the worst discrepancies is in the experimental value 
of 0; this may come from some uncertainty in the temperature scale, as mentioned 
earlier [195]. 

Three of the X-ray results are presented in a slightly different manner in Fig. 20. 
Each of the experiments shown is tentatively assumed correct and used to calculate 
a value of /,//, with the aid of the least squares data for the other constants. The 
ruled grating value as well as the least squares value of 4,/4; are given for comparison. 
Nilsson’s measurements of excitation potentials are in excellent agreement with the 
least squares value. The results of the electron diffraction work may be considered 
in fair agreement in view of its probable error. The same holds for the measurement 
of the X-ray high frequency limit in the corrected version adopted in the present report. 

However, the corresponding uncorrected value with its original weight is seen to 
disagree drastically with the least squares result. It is believed that the reasons pre- 
viously cited [196] account for most of the discrepancy. Use of the uncorrected figure 
with its original weight obviously accounts for the high values of 4,/A, obtained in 
previous evaluations. If the reasoning of this report is valid, such a procedure must 
now be considered erroneous and the resulting high values of 4,/4 ignored. 
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3. — Recommendations for new experiments. 


a) The principal discrepancy in the present adjustment seems to be 
that between the electrochemical value of the faraday and that implied by 
the other experiments. The most likely causes seem to lie in the uncertainty 
as to the purity of the sodium oxalate used by CRAIG and HOFFMAN and the 
uncertainty in the isotopic abundances and therefore the atomic mass of silver. 
Further refinement of the CRAIG and HOFFMAN work, as well as the necessary 
measurements on silver would seem to be of first importance. 


b) New theoretical work leading to a better interpretation of the observed 
isochromats in the X-ray «h/e experiments » would be highly desirable. If 
properly interpreted, these experiments might again play an important part 
in fixing the atomic constants. 


e) There is still some evidence for a possible variation in the velocity 
of light with frequency, although the negative case appears much stronger. 
‘One more optical measurement confirming Bergstrand’s value would do much 
to settle any lingering doubts. In this connection, the possibility of refining 
Houston’s method deserves study. 


d) In view of the very limited amount of overdetermination in the 
present adjustment, any new precision experiments (with errors of the orders 
of 10 ppm) would be of great value. For example, SIEGBAHN’s [197] proposed 
method of measuring h/e by B-ray spectroscopy warrants careful consideration. 

In the final analysis, improvement in our knowledge of the atomic con- 
stants depends primarily on new and more accurate experiments. If this re- 
port serves in any way to clarify the present situation and to encourage new 
work along useful lines, it has fulfilled a major purpose. 


The present authors were associated with Dr. M. D. EARLE and Mr. J. M. 
MINKOWSKI in a previous study of the atomic constants. Dr. EARLE and Mr. 
MINKOWSKI carried out a great deal of the preliminary work in connection 
with the present investigation. The authors are also greatly indebted to them 
for many helpful discussions during the later stages of the work; in particular, 
they wish to thank Mr. Minkowski for a careful reading of the manuscript 
and for his valuable suggestion as to the cause of the discrepancy in the 
X-ray «h/e experiments ». 

It is probably impossible at the present time to compile a complete list 
of those who have contributed significantly to this work. However, the writers 
are particularly indebted to the following workers: Drs. J. I. HOFFMAN and 
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D. N. Crate for several valuable discussions on the electrochemical measure- 
ment of the Faraday; Dr. J. A. HIpPLE ifor detailed {discussion of the mea- 
surements of the gyromagnetic ratio of the proton and the cyclotron frequency 
of the proton, as well as the problem of the best value of the Faraday; 
Mr. R. DRIScoLL for informative conversations on the absolute determination 
of the ampere and for reviewing some of the work on the gyromagnetic ratio 
of the proton to answer questions raised by the authors; Dr. J. L. THOMAS. 
for recalibrating certain electrical standards in a search for the possible cause 
of the h/e discrepancy; Dr. D. H. RANK for a full discussion of his measure- 
ment of ¢ before publication; Dr. E. S. DAyHorr for many informative con- 
versations on the fine structure measurement and the determination of €; 
Dr. H. M. WATTS for his assistance in many ways during the early part of the 
study; Dr. E. E. SALPETER for valuable advice on the interpretation of the 
fine structure and hyperfine structure measurements; Dr. H. Lyon for helpful 
discussions of the velocity of light; and Dr. T. H. BERLIN for reviewing and 
criticizing the section on quantum electrodynamie experiments as well as his 
assistance in other matters. 

Since the original report was distributed they have had helpful discussions 
from Commander C. I. AsLAKSON, Dr. L. EssEN, Dr. R. A. HousToUN, and 
the Rev. J. F. MULLIGAN, on the velocity of light. Dr. J. H. SANDERS and 
Dr. K. R. TRIGGER have kindly furnished their respective results on the cy- 
clotron resonance frequency of the proton in advance of publication. Dr. R. D. 
HUNTOON and his associates at the Bureau of Standards have been most co- 
operative in describing the current status of various experiments at the Bureau. 
Dr. J. W. M. DuMonp has kindly furnished us with advance copies of a 
number of studies by himself and his coworkers. 

The writers also wish to thank Mr. O. RASK for help in the computations; 
Mr. W. E. TREXLER and Mr. M. LA PORTE for preparing the figures; Mrs. M. B. 
Rowe and Mrs. W. REICHELT for typing the manuscript; and Mrs. J. S. THOMSEN 
for assistance in proof reading. To others who have inadvertently been omitted 
from this list the authors extend their apologies and their thanks. 


APPENDIX A 


Recent developments 
1. — Velocity of light. 
SILSBEE and GARDNER [198] have recalculated the results of the Rosa and 


Dorsey experiment, correcting it for humidity of the air and for the best con- 
version factors between the old international electrical units and absolute 


' 
j 
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units. They obtained (299800 + 31) km/s; this would raise the old Birge 
average by about 1 km/s. RANK [199] has recently revised his value of c as a 
result of a more precise determination of the band origin used in his infra-red 
work; his new figure is (299 791.9 + 2.2) km/s. HousToUN [200] has recalculated 
his result, using a more precise correction for the refractive index of air, and 
obtained 299 782 km/s. All of these revisions bring the respective experiments 
in better agreement with the value recommended in the present report. 

HENDERSON and MUNGALL [201] have reported a resonant cavity measure- 
ment at 400 MHz, which gave c = 299780 km/s. However, they concluded 
that various systematic errors were sufficient to explain the discrepancy be- 
tween their result and the accepted value. 

Three new precision measurements [202] of ¢ have been carried out with 
the Bergstrand geodimeter, in Sweden, Great Britain, and Australia respec- 
tively. All of the results were about .8 km/s lower than Bergstrand’s. Although 
these results could conceivably all have some systematic errors in common 
with Bergstrand, it is difficult to imagine a significant one which would be 
essentially independent of length of base line. Furthermore, in the case of 
the British measurement the frequency of the crystal was independently 
calibrated. 

While in the opinion of the authors it would still be desirable to have 
another completely independent precision measurement at an optical frequency, 
preferably in vacuum, the need for such confirmation has certainly lessened. 
As for the recommended value, the lowering of the geodimeter value is par- 
tially compensated by the rise in Rank’s result; the weighted mean is not. 
changed significantly. 


2. — Avogadro’s number (X-ray crystal density method). 


SMAKULA, KALNAJS and SILS [203] have completed a series of measurements. 
of erystal densities and lattice spacings (by X-ray measurements). The best 
results were obtained with Si, Al, CaF, and Ge; these gave values of Avogadro’s. 
number mutually agreeing within 10 ppm, with the average value of N’ about 
80 ppm higher than the calcite figure used by the authors. (An error in the 
calculated calcite value might be explained if the impurity atoms occur inter- 
stitially rather than substitutionally, as assumed in the Birge and Brogren 
corrections). If the accuracy is as high as indicated by the excellent internal 
agreement for different crystals, this average should probably be adopted to 
the exclusion of calcite data. The only effect of major significance would be 
on A,/A,. The value implied by N’ would then be 1.00205 and the weighted 
average about 1.00204. 

However, there are possible sources of systematic errors common to all 
these measurements, both in the X-ray spectrometer and in the density mea- 
surements. (For instance, the density of quartz was determined as an inter- 
mediate step in the weighing; this density was in serious disagreement with 
some of the previous measurements, e.g., DuMond’s). Thus it seems premature 
to put the entire weight on these values. In any event, the change in 4,/4. 
would only be about half its assigned probable error. 
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3. — Mass of the proton. 


QUISENBERRY, SCOLMAN and NIER [204] have performed a new and highly 
precise measurement of the mass of the proton, obtaining My= (1.0081442 + 
+ .0000002) amu. A new evaluation of nuclear data by WAPSTRA [205] leads 
to a value My = 1.0081452 amu. Both results fall within the probable error 
of the value recommended here. 


4. — X-Ray «h/e experiments ». 


KERSCHER and KULENKAMPFF [206] and SCHEER and ZEITLER [207] have 
made experimental and theoretical studies respectively of thin-target X-ray 
spectra; however, they are concerned primarily with angular distribution. 
ALBERT [208] has carried out further work on the correlation between discrete 
electron losses and the form of the X-ray isochromat near cut-off frequency; 
his theory now includes a qualitative explanation of the observed minima. 
Recent developments in measurement of discrete electron losses and their 
correlation with X-ray fine structure has been well summarized by LEDER, 
MENDLOWITZ and MARTON [209]. While all of this new work seems to strength- 
en the analysis of the X-ray «h/e experiments » presented in this report, 
it would not appear that it is yet possible to interpret the results of these 
experiments with any higher accuracy. 


5. — Gyromagnetic ratio of the proton. 


The gyromagnetic ratio of the proton has been measured by KIRCHNER 
and WILHELMY [210] in a new experiment employing an iron-free field, who 
obtained a value y, = 2.67562 + .00016. This figure is about 150 ppm higher 
than that of THOMAS, DRIScOLL and HIPPLE, the difference being considerably 
greater than the sum of the errors. They suggest that the difference may be 
due to a discrepancy in electrical conversion factors and plan to continue their 
work. Meanwhile at the Bureau of Standards in Washington work is well 
under way on a new and more precise determination, also using an iron-free field. 


6. — Cyclotron resonance frequency of the proton. 


TRIGGER [211] has refined the Jeffries experiment with the inverted cyclotron. 
The effect of impurities and of stray electric and magnetic fields was carefully 
investigated. Furthermore he developed a more exact theory of the operation 
of the cyclotron, including certain non-linear terms in his equations and solving 
them with the aid of an IBM computer. The preliminary result indicates 
ssltn = 2.79267 + .00010, without diamagnetic correction. Still more refined 
calculations are in progress. His work essentially answers the criticism of the 
Jeffries value advanced in this report. 
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Another precise measurement of the same quantity has been reported by 
COLLINGTON, DELLIS, SANDERS and TUBERFIELD [212]. They employed a mo- 
dified form of inverted cyclotron with two grounded outer dees and a straight 
center conductor connected to a high frequency source. This was designed 
so that at resonance all protons approached asymptotically a particular cir- 
cular orbit, in which the transit time between the dees was exactly one cycle 
of the rf. source. They obtained us/u, = 2.79273 + .00004, without diamagnetic 
correction; the error is intended as a limit of error. 

Both of these results are essentally in agreement with the Sommer, Thomas 
and Hipple value used in the present study. A properly weighted average 
might now give about 2.79270; in view of the external consisteney of the 
data the present error could not be appreciably decreased. The effect on most 
of the other constants would be very small and, in consideration of the pro- 
isional character of the Trigger value, there seems as yet no reason for a 
change. 


‘7. — Hyperfine structure of hydrogen. 


The hyperfine splitting of H has been remeasured with improved apparatus 
by KuscH [213], who now reports Av(H)= (1420.405 73 + .000 05) MHz, in essen- 
tial agreement with the result of Dicke. However, the theoretical interpretation 
of the hfs. still appears unsatisfactory. More precise calculations of the 
«average excitation potential» have been carried out by HARRIMAN [214]. 
MOELLERING, ZEMACH, KLEIN and Low [215]have used the hfs. result to estimate 
the spatial extent of the proton magnetic moment; this i just the reverse of 
the procedure which would be desirable from the standpoint of evaluation of 
atomic constants. 

HEBERLE, REICH and KuscH [216] have determined the hfs. splitting in 
the metastable 2S state; the ratio Av(2S)/Av(18S) checks with present theoretical 
predictions within about 1 ppm. However, it would appear that this new 
datum has the same uncertainty in interpretation as for the 18 state, in so 
far as the calculation of the fine structure constant « is concerned. 

More precise theoretical knowledge would be of great value since without 
it the highly weighted fine structure determination of x stands virtually alone, 
without any other data of comparable accuracy to confirm or disprove it. 
Of course, new microwave fine-structure measurements might serve the same 
purpose. 


8. — 6-Ray spectroscopic measurements. 


SIEGBAHN and EpVARSON [217] have recently reported a value of (h/me) = 
=(2.4262-4-.0016)-10-1° em, obtained by f-ray spectroscopy on electrons ejected 
by internal conversion of f-rays. When a number of lines are measured and 
the energy differences of various initial electronic states known from X-ray 
data, the necessary relativistic equations can be solved for h/me. In this pro- 
cedure errors are compounded in such a way that the accuracy is low; while 
the individual lines are measured to about 10 ppm, the above error is given 
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as almost 700 ppm. Thus it can not as yet carry significant weight in evaluating 
atomic constants. However, it is in almost perfect agreement with the values 
recommended here. It may also be noted that SIEGBAHN [197] has proposed 
a rather similar method for determining h/e. 


9. — Other evaluations and reviews. 


STILLE [218] has carried out studies of various atomic constants, leading 
to recommended values for c, F, e, N and k. Within their errors, these values 
agree with those recommended here; but the errors are much higher. More 
recently he has published a book [219] on physical units and constants, but 
the writers have not yet had an opportunity to study this work. Mention 
should also be made of a review of the velocity of light by A. KAROLUS [220]. 

Almost simultaneously with the completion of the body of the present, 
report, in August 1955, CoHEN, DUMOND, LAYTON and ROLLETT [221] published 
a new adjustment of the atomic constants. This was primarily based on data 
which they used in their 1952 evaluation, but involved a systematic analysis 
of variance in an effort to eliminate experiments with unsuspected systematic 
errors. It was reassuring to find that their conclusions in regard to the X-ray 
« h/e experiments » were substantially the same as reached by the writers. They 
further pointed out that in order to make all three « h/e experiments » agree 
satisfactorily the form of the ideal isochromat must be assumed to be voltage 
dependent, a point not properly emphasized by the present authors. On the 
other hand, even after offering a tentative explanation of the systematic error, 
they retained the least accurate of these three experiments in the final eval- 
uation. Although the practical effect of this is small, it certainly appears. 


more logical to reject all three of the questionable items. (Thus the adjust- 


ment denoted as 1111120 in Table IV of their paper would appear to be 
the best choice. After allowing for the fact that it has one less degree of 
freedom than the one they adopt, the value of y? is still considerably more 
favorable). 

They also differed from the present analysis in the use of electrochemical 
data. They used only silver and iodine values, with the latter based on an 
older value of atomic mass; furthermore the errors quoted for these meas- 
urements appear very low. As a result, they retained the iodine value, which 
is in almost perfect agreement with electromagnetic experiments, and rejected 
the silver one. (The writers considered the same thing during the early stages 
of their work). However, with the iodine result revised for the latest value 
of atomic mass and the Craig and Hoffman figure included, there are three 
electrochemical values in fair agreement with one another, all lower than the 
electromagnetic value. Hence it now appears very questionable to completely 
reject the silver result. As a result, the CDLR [221] value of faraday differs from 
that recommended here by almost three times the probable error of the former; 
there are lesser effects on the values of other constants. 

Finally they employed Birge’s average for the measurement of Avogadro’s 
number by means of the X-ray crystal density method; this gives a figure 
for N’ about 60 ppm higher than the revised calcite value used by the writers. 


However, as noted earlier in this Appendix, the SMAKULA, KALNAJS and SILS. 
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data indicates that the Birge average should be raised rather than lowered. 
Thus the CDLR evaluation seems to occupy a safe middle ground in this 
respect. At any rate the practical effect of this difference is small, except 
perhaps on the recommended value of 4,/4;. 

CoHEN [222] recently developed a more precise method for evaluating the 
significance of differences between items of input data and least-squares output 
values. He concluded that none of the items retained in the CDLR study 
(except possibly the BEARDEN and ScHWARZ «h/e experiment ») showed signi- 
ficant differences. 


APPENDIX B 


Summary of the theory of least squares 


1. — Most probable values of the variables. 


Consider a set of N experiments to determine n variables, where N > n. 
Assume that all of the experimental equations can be linearized as described 
in Sect. C. 3. Then the i-th experiment is described by the equation 


(B.1) >, Carly. 76 « 1<i<N. 


The matrix a is thus a rectangular (N row x» column) matrix. Assume further 
that the experiments are completely independent (no correlation of errors) 
and that the errors have a Gaussian distribution. Let o; be the standard 
deviation of the i-th equation. 

Since the set of equations is overdetermined, it will in general be impos- 
sible to find a set of xw’s which satisfy all equations simultaneously. Let 7; be 
the difference between the right and left sides of Eq. (B.1) for a given set 
of ws. On the basis of the i-th experiment the probability of these values 
being correct is proportional to 


exp [— ri/20î] . 
Since the experimental errors are all independent, then on the basis of all 


the experiments the probability of a given set of «’s being the correct one 
is simply proportional to a product of N such factors, viz., 


To find the best compromise value (the most likely value) the above prob- 
ability must be maximized, i.e. the sum must be minimized; hence the name 
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«least squares » is used. In terms of the a’s this sum becomes 


The factor (1/07) may equally well be taken as K/o?, where K is any convenient 
constant; this will not affect the values of the 2s which give a minimum. 

To minimize the sum it must be differentiated with respect to each # in 
turn and each derivative set equal to zero. The p-th equation obtained in 


this manner is 


N K n i 
(B.3) > a) dip.(Da:2—)=0. 1<p<n. 
ped 


isl a 


There are nm such equations known as «normal equations »; these uniquely 
determine the best values of the #’s. The normal equations may be put in 
the form 


(B.4) > Int, = hy , 
r=1 
where 
Me K 
(B.5) Yor = > (È) AipAir , h, = Oy (È) Dipl; . 
al \LUf t=1 i 


It is clear that the matrix g is a square symmetric nxn matrix. The solution is 
(B.6) Ly = > Irony ’ 
p=1 


where g;; is the rp-th element of the reciprocal of the matrix g. 


2. — Probable errors of the variables. 


Let e, be the error associated with the calculated value of x;: Since e, is 
equally likely to be positive or negative, its average value is zero. However, 
<e?>, the average value of the square of ¢,, is non-zero and will now be com- 
puted. As the average of the product of errors, <¢,¢,> will also be needed, it 
will be calculated at the outset; <e?> may then be obtained as a special case 
by putting r =t. 

In the following analysis, è and j will be used for indices ranging from 
1 to N (number of experiments), while p, r, s and ¢ will denote indices ranging 
from 1 to n (number of variables). Combining Eqs. (B.6) and (B.5) gives 


n n N K 
(B.7) a =D gh = DI X (3) DiGi - 
p=1 p=1 t 


i=1 
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Each constant c; is, of course, equal to its true value plus some error 6;. The 
corresponding error in @, is 


(B.8) => a8 IE i) iad 
Hence 


(B.9) 6,6,= > di > (i 1) ud] Zan S (5) Dy Oil 


p=1 j=1 È) 


Since all factors on the right are known constants except 6; and ò;, the average 
value is given by 


Sodi © (E\(K 
(B.10) CEE x = 2 ZI Yro Its I Sla o o inh 350i Oj) « 

Dp j î 
As the experiments are assumed completely independent, their errors are 
uncorrelated. Therefore <6,6;> = 6,,0; where 6,; is the Kronecker 6. Thus 


(B.11) (EE: = x S Gro Gis > LE i) VipWis » 


p=1 s=1 


From the definition of g in Eq. (B.5) it follows that the last sum is equal 
to Kg,,. Therefore 


(B.12) Cele War IY “> Us, Gaps 
p=1 s=1 
Since > 9:9. = dm, the Kronecker 6, it follows that 
s=1 
(B.13) Kee = Kga' - 


Thus <> = Kg. Since the probable error is simply a constant times the 
standard deviation, this same method will give probable error if the experi- 
ments are weighted as A/(P.E.)? instead of K/o?. 

It must be emphasized that the above errors have been computed by 
internal consistency, i.e., on the basis of the errors estimated for each exper- 
iment. It may now be asked whether or not the agreement between different 
experiments is consistent with these estimates. This question is answered 
by the x? test. 


3. — The y? test. 


The y? test involves the expected value of the sum in Eq. (B.2) when the 
least square solution for the «’s is substituted in it; in other words, it is the 
value of this sum after it has been minimized. If N > n, the least squares 
solution should approximate the true values of the variables; each term would 
then have an average value of unity and the sum would be approximately N. 
However, the case when N and n are comparable is of great practical interest. 
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Let d; be the value of r; when the least squares solution is substituted for the as. 
Let 6; now denote the error in c;, which is, of course, unknown, while X; 
denotes the true value of «;, also unknown. Then, since the true value of 
the constant is the observed ¢; minus the error 6;, it follows that 


| Dd Girt, — Ci =d, 
(B.14) pai 
| DÀ Gir XK, — Ci e 0:5 
r=1 
(B.15) Sa, ee ee 


The quantity x, — X, is simply the error in 2,, i.e., €, 
(B.16) ieee d; = 0; “3 > Oe onus 
r=1 


The expected value of 7? now becomes 


N dî n 
(B.17) (x? =(39)- DEI 0; Zare), 


(B.18) 4) = i 4) -2($ > oS a) (Za Save, > ei) . 


i=1r= i Gig hres 


These terms will now be considered separately. Since the average value of 
<6;> is oî, it follows that the first term is 


The second term may be simplified by use of Eq. (B.8) and by recalling that 
<6,6;> = 6;;07 (where 6,; is the Kronecker 6) since the errots are independent. 


n Wein ae Oe nt N K 
(B.20) x Di dadi i :) — © 3 = 2 > dro (5) ad) Fe 


i=1r=1 Oj j 


K 


aoe £6; 6; = FE Sì Sg > (5) Uirlip » 


r=1p=1 tel i 


By Eq. (B.5) the last sum is simply g,,. Hence, 


eo n n 3 n 
(B.21) > x Sa ai > > Ore Yer 53 >? Ors TAM 
r=1 


i=l1lf=1 r=1p=1 
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The final term in the sum may be evaluated with the help of Eqs. (B.13) and (B.5) 


Eon 1 r=1t=1 is1 di r=1t=1 


(B.23) Se > UirEr > U5 Et rs S Orr =n. 
r=1 


i Caen 


Eqs. (B.19), (B.21) and (B.23) are now substituted in Eq. (B.18) to give 
(B.24) (x) = N—-2n+4n=N—-n. 


This is the desired result. If in a given case x > N — n, it indicates that 
some or all of the experimental errors have been seriously underestimated; 
conversely, if y:< N — n, some or all of the errors have been stated much 
too conservatively. 

It should be noted that standard deviations, not probable errors are used 
in the definition of y?. If probable errors are used elsewhere, the sum may 
be taken on the basis of probable errors and then multiplied by (.6745)?. 


4. — Transformation and elimination of variables. 


The sum in Eq. (B.2) to be minimized is a function defined in n-dimensional 
space. If the function is expressed in terms of different coordinates, say y’s, 
this can not affect the geometrical point at which the function reaches its mi- 
nimum. Hence it may be first expressed in terms of the y’s and then mini- 
mized; when the corresponding a’s are calculated from the solution in the y’s, 
the result must be the same as if the entire problem had been solved in terms 
of the x°s. Thus it is only a question of convenience as to what quantities 
are taken as the independent variables. 

When one equation has a weight very much higher than the others (e.g. the 
equation for R.), it may be considered as an exact relationship defining one 
variable in terms of the others, i.e. reducing the number of independent va- 
riables by one. Since the final result is independent of a transformation of 
variables, it is a matter of indifference which variable is eliminated by this 
relation. 


5. — Non-Gaussian error distributions. 


In the case when the error distributions are not Gaussian (a most important 
practical case), COHEN [188] has shown that the least squares method can be 
derived by the Axiom of Minimum Error, i.e., the condition that the error 
should be a minimum with respect to possible variation of free parameters. 
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